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Chapter 1 
Introduction 



These lectures were delivered at the ECT* 1 in Trento (Italy) in 2002 and 
2003. They are addressed to physicists who wish to acquire a minimal back- 
ground to understand present day attempts to model the confinement of 
QCD 2 in terms of dual superconductors. The lectures focus more on the 
models than on attempts to derive them from QCD. 

It is speculated that the QCD vacuum can be described in terms of a 
Landau-Ginzburg model of a dual superconductor. Particle physicists often 
refer to it as the Dual Abelian Higgs model. A dual superconductor is a 
superconductor in which the roles of the electric and magnetic fields are ex- 
changed. Whereas, in usual superconductors, electric charges are condensed 
(in the form of Cooper pairs, for example), in a dual superconductor, mag- 
netic charges are condensed. Whereas no QED 3 magnetic charges have as 
yet been observed, the occurrence of color-magnetic charges in QCD, and the 
contention that their condensation would lead to the confinement of quarks 
was speculated by various authors in the early seventies, namely in the pio- 
neering 1973 paper Nielsen and Olesen [I], the 1974 papers of Nambu [2] and 
Creutz jsl], the 1975 papers of 't Hooft |4\, Parisi [H], Jevicki and Senjanovic 

and the 1976 paper of Mandelstam Qualitatively, the confinement of 
quarks embedded in a dual superconductor can be understood as follows. The 
quarks carry color charge (see App|Dj). Consider a static quark- ant iquark 
(qq) configuration in which the particles are separated by a distance R. The 

1 The ECT* is the European Centre for Theoretical Studies in Nuclear Physics and 
Related Areas. 

2 QCD: quantum chromodynamics. 
3 QED: quantum electrodynamics. 
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quark and anti-quark have opposite color-charges so that they create a static 
color-electric field. The field lines stem from the positively charged particle 
and terminate on the negatively charged particle. If the qq pair were em- 
bedded in a normal (non-superconducting) medium, the color-electric field 
would be described by a Coulomb potential and the energy of the system 
would vary as —e 2 /R where e is the color-electric charge of the quarks. How- 
ever, if the qq pair is embedded in a dual superconductor, the Meissner effect 
will attempt to eliminate the color-electric field. (Recall that, in usual super- 
conductors, the Meissner effect expels the magnetic field.) In the presence of 
the color-electric charges of the quarks, Gauss' law prevents the color-electric 
field from disappearing completely because the flux of the electric field must 
carry the color-electric charge from the quark (antiquark) to the antiquark 
(quark). The best the Meissner effect can do is to compress the color-electric 
field lines into a minimal space, thereby creating a thin flux tube which joins 
the quark and the antiquark in a straight line. As the distance between the 
quark and antiquark increases, the flux tube becomes longer but it maintains 
its minimal thickness. The color-electric field runs parallel to the flux tube 
and maintains a constant profile in the perpendicular direction. The mere 
geometry of the flux tube ensures that the energy increases linearly with R 
thereby creating a linearly confining potential between the quark and the 
antiquark. This qualitative description hides, in fact, many problems, which 
relate to abelian projection (described in ChapEJ), Casimir scaling, etc. As a 
result, attempts to model quark confinement in terms of dual superconduc- 
tors are still speculative and somewhat ill defined. 

The dual superconductor is described by the Landau- Ginzburg (Dual 
Abelian Higgs) model ||. Because the roles of the electric and magnetic 
fields are exchanged in a dual superconductor, it is natural to express the 
lagrangian of the model in terms of a gauge potential B^ associated to the 
dual field tensor F^ v = d^B v — d u B^. Indeed, when the field tensor F is 
expressed in terms of the electric and magnetic fields , as in Eq. (|2.ip . the 
corresponding expression ()2.2j) of the dual field tensor F^ v = e^^F^ is 
obtained by the exchange E — * H and H — > —E of the electric and magnetic 
fields. (Recall that electrodynamics is usually expressed in terms of the gauge 
potential N 1 associated to the field tensor F^ = d»A u - d u A fM .) When a 
system is described by the gauge potential B^ , associated to the dual field 
tensor F, the coupling of electric charges (such as quarks) to the gauge field 
B^ is analogous to the problem of coupling of magnetic charges in QED to 
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the gauge potential A^. Such a coupling was formulated by Dirac in 1931 
and 1948 |Hj and it requires the use of a Dirac string. The Dirac theory of 
magnetic monopoles is reviewed in ChapJU In Sect l2.ll! it is applied to the 
coupling of electric charges to the gauge field associated to the dual field 
tensor F. 

For a system consisting of a qq pair, the Dirac string stems from the quark 
(or antiquark) and terminates on the antiquark (or quark). The string should 
not, however, be confused with the flux tube which joins the two particles 
in a straight line. Indeed, as explained in Sect 12. 91 the Dirac string can 
be deformed at will by a gauge transformation. The latter does not modify 
the flux tube, because it is formed by the electric field and the magnetic 
current, both of which are gauge invariant. We refer here to the residual 
U (1) symmetry which remains after the abelian gauge fixing (or projection), 
discussed in ChaplU However, as explained in ChaptJ31 there is one gauge, 
the so-called unitary gauge, in which the flux tube forms around the Dirac 
string. Calculations of flux tubes have all been performed in this gauge. 

One attempt [TUj to apply the dual superconductor model to a system of 
three quarks is discussed in ChapEJ Ultimately, this is the goal aimed at by 
these lectures. We would like to formulate a workable model of baryons and 
mesons, which would incorporate both confinement and spontaneous chiral 
symmetry breaking and which could be confronted to bona fide experimental 
data and not only to lattice data. Presently available models of hadrons 
incorporate either confinement or chiral symmetry, but not both. It is likely 
that models, such as the one described in ChaplHl will have to be implemented 
by an interaction between quarks and a scalar chiral field, for which there is 
also lattice evidence [TT] . 

The model is inspired by (but not derived from) several observations made 
in lattice calculations. The first is the so-called abelian dominance, which 
is the observation that, in lattice calculations performed in the maximal 
abelian gauge, the confining string tension a, which defines the asymptotic 
confining potential o~R, can be extracted from the Abelian link variables alone 
H2lCI3|,|n|,IISl,[E|,in|,[ISl- Abelian gauge fixing is discussed in Chapg] 
both in the continuum and on the lattice. 

The second observation, made in lattice calculations, is that the confining 
phase of the SU (N) theory is related to the condensation of monopoles 
12DI EI] > [22H2HI- Such a statement can only be expressed in terms of an 
abelian gauge projection. The condensation of monopoles and confinement 
are found to disappear at the same temperature and it does not depend on the 
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chosen abelian projection (213 122]- However, confinement may well depend on 
the choice of the abelian gauge. In the abelian Polyakov gauge, for example, 
monopole condensation is observed but not confinement [2E]. in ChapEJ we 
show how monopoles can be formed in the process of abelian projection. It 
is often difficult to assess the reliability and the relevance of lattice data. For 
example, on the lattice, even the free U (1) gauge theory displays a confining 
phase in which magnetic monopoles are condensed [23 |2H] • This confining 
phase disappears in the continuum limit [2H] as it should, since a U (1) gauge 
theory describes a system of free photons. However, non-abelian gauge theory 
is better behaved than U (1) gauge theory (it is free of Landau poles) and 
lattice calculations point to the fact that, in the non-abelian theory, the 
confining phase, detected by the area law of a Wilson loop, survives even in 
the continuum limit. 

The third observation, which favors, although perhaps not exclusively, 
the dual superconductor model, is the lattice measurement of the electric 
field and the magnetic current, which form the flux tube joining two equal 
and opposite static color-charges, in the maximal abelian gauge |30].|35]. 
jHIl E21 E3] , |M| • They are nicely fitted by the flux tube calculated with the 
Landau-Ginzburg (Abelian Higgs) model, as discussed in Sect 13.41 

The model is, however, easily criticized and it has obvious failures. For 
example, it confines color charges, in particular quarks, which form the funda- 
mental representation of the SU (N) group and therefore carry non-vanishing 
color-charge. However, it does not confine every color source in the ad- 
joint representation: for example, it would not confine abelian gluons. (Color 
charges of quarks and gluons are listed in AppjDj) Because it is expressed 
in an abelian gauge, the model also predicts the existence of particles, with 
masses the order of 1 — 2 GeV, which are not color singlets. 

In addition, there is lattice evidence for competing scenarios of color 
confinement, which involve the use of the maximal center gauge and center 
projection, described in Sect J4. 31 They are usefully reviewed in the 1998 and 
2003 papers of Greensite jHHlEZl- They account for the full asymptotic string 
tension as well as Casimir scaling. In fact, both the monopole and center 
vortex mechanisms of the confinement are supported by the results of lattice 
simulations. They are related in the sense that the main part of the monopole 
trajectories lie on center projected vortices [SB], [HHl- We do not describe the 
center- vortex model of confinement in these lectures because it does not, 
as yet, lead to a classical model, such as the Landau-Ginzburg (Abelian 
Higgs) model. Instead, it describes confinement in terms of (quasi) randomly 
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distributed magnetic fluxes in the vacuum. It is however, numerically simpler 
on the lattice and flux tubes formed by both static qq and 2q2q have been 
computed [HI] ■ Further scenarios, such as the Gribov coulomb gauge scenario 
developed by Zwanziger, Cucchieri [HJI12] anc ^ Swanson [13], and the gluon 
chain model of Greensite and Thorn are not covered by these lectures. 

The relevant mathematical identities are listed in the appendices. 
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Chapter 2 



The symmetry of 
electromagnetism with respect 
to electric and magnetic 
charges 

The possible existence of magnetic charges and the corresponding electro- 
magnetic theory was investigated by Dirac in 1931 and 1948 [JHl E] The 
reading of his 1948 paper is certainly recommended. A useful introduction 
to the electromagnetic properties of magnetic monopoles can be found in 
Sect. 6. 12 and 6.13 of Jackson's Classical Electrodynamics [47J. The Dirac 
theory of magnetic monopoles, which is briefly sketched in this chapter, will 
be incorporated into the Landau- Ginzburg model of a dual superconductor, 
in order to couple electric charges, which ultimately become confined. This 
will be done in ChaptEl 

2.1 The symmetry between electric and mag- 
netic charges at the level of the Maxwell 
equations 

"The field equations of electrodynamics are symmetrical between electric 
and magnetic forces. The symmetry between electricity and magnetism is, 
however, disturbed by the fact that a single electric charge may occur on a 
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particle, while a single magnetic pole has not been observed on a particle. 
In the present paper a theory will be developed in which a single magnetic 
pole can occur on a particle, and the dissymmetry between electricity and 
magnetism will consist only in the smallest pole which can occur, being much 
greater than the smallest charge." This is how Dirac begins his 1948 paper 

— * — * 

The electric and magnetic fields E and H can be expressed as components 
of the field tensor F^ v : 



( 


—E x 


-Ey 


-E z \ 


E x 





-H z 


Ey 


Ey 


H z 





—H x 


\E Z 


~Hy 


H x 





(2.1) 



They may equally well be expressed as the components of the dual field 
tensor F^ v : 

( -H x 



p» v 



1 



a/3 







E, —E„ 



-H z \ 



Hy —E Z 







V H z 



Ey —E X 



E x 
I 



(2.2) 



where the antisymmetric tensor with 5 0123 = 1. Thus, the cartesian 

components of the electric and magnetic fields can be expressed as compo- 
nents of either the field tensor F or its dual F: 



E' 



1 -e 0ijk F ik H l = -F 0i = --e 0ijk F, 



3k 



(2.3) 



_ F oi = _ e «ifc_p ip = _ F oi z 

2 3 2 

The appendix 1X1 summarizes the properties of vectors, tensors and their dual 
forms. In the duality transformation F — > F, the electric and magnetic fields 
are interchanged as follows: 



E -> H H 



-E 



(2.4) 



The electric charge p and the electric current j are components of the 
4- vector j M : 

f=(pj) (2.5) 



Similarly, the magnetic charge p mag and the magnetic current j mag are com- 



ponents of the 4- vector j^ as : 



a-' 



Jmag ~ \ Pmag, Jmag 



(2.6) 
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At the level the Maxwell equations, there is a complete symmetry between 
electric and magnetic currents and the coexistence of electric and magnetic 
charges does not raise problems. The equations of motion for the electric 
and magnetic fields E and H are the Maxwell equations which may be cast 
into the symmetric form: 

d v F^ = f d u T" = f mag (2.7) 

It is this symmetry which impressed Dirac, who probably found it upsetting 
that the usual Maxwell equations are obtained by setting the magnetic cur- 
rent jf^ to zero. The Maxwell equations can also be expressed in terms of 

the electric and magnetic fields E and H. Indeed, if we use the definitions 
(EHl) and d22l, the Maxwell equations (j22|) read: 

q vF uh = f _^ V-E = p -d t E + V x H = j 

d v F vl "=f mag -> V-H = Pmag -d t H-V xE = j mag (2.8) 

2.2 Electromagnet ism expressed in terms of 
the gauge field associated to the field 
tensor F^ u 

So far so good. Problems however begin to appear when we attempt to ex- 
press the theory in terms of vector potentials, alias gauge potentials. Why 
should we? In the very words of Dirac 0: "To get a theory which can be 
transferred to quantum mechanics, we need to put the equations of motion 
into a form of an action principle, and for this purpose we require the elec- 
tromagnetic potentials." 

This is usually done by expressing the field tensor F^ v in terms of a vector 

potential = A) : 

F nv = d n A u _ qv A h g = _ dt £ H = f x A (2.9) 

However, this expression leads to the identity 1 d„F u = which contradicts 
the second Maxwell equation d v F = j^ ag . The expression F = d A A 

1 The identity d ■ d A A — is often referred to in the literature as a Bianchi identity. 
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therefore precludes the existence of magnetic currents and charges. In elec- 
tromagnetic theory, this is a bonus which comes for free since no magnetic 
charges have ever been observed. Dirac, however, was apparently more se- 
duced by symmetry than by this experimental observation. 

Let us begin to use the compact notation, defined in AppEJ and in which 
d A A represents the antisymmetric tensor (d A A)' 11 ' = d^A u — d v N 1 . The 
reader is earnestly urged to familiarize himself with this notation by checking 
the formulas given in the appendix El lest he become irretrievably entangled 
in endless and treacherous strings of indices. 

Dirac proposed to modify the expression F = d A A by adding a term 

—G: 

F = dAA-G F = dAA + G (2.10) 

where G^ v = —G Uf1, is an antisymmetric tensor field 2 . The latter satisfies the 
equation: 

d-G = Jmag (2.11) 

The field tensor F^ v then satisfies both Maxwell equations, namely: d-F = j 
and d ■ F — j mag . In the expressions above, the bar above a tensor denotes 
the dual tensor. For example, G^ u = \sn Ua pG a P (see AppEJ). For reasons 
which will become apparent in Sect. 12.51 we shall refer to the antisymmetric 
tensor G,„, Dirac string term. 

The string term G^ v is not a dynamical variable. It simply serves to 
couple the magnetic current j!^ ag to the system. It acts as a source term. 
Note that both G and the equation d-G = j mag are independent of the gauge 
potential A^. 

An equation for A^ 1 is provided by the Maxwell equation d-F — j. When 
the field tensor F has the form ()2.10|) . the equation reads: 

d- (dAA) -d-G = j (2.12) 

The Maxwell equation ()2.12|) may be obtained from an action principle. In- 
deed, since the string term G does not depend on the gauge field A, the 
variation of the action: 

W (A) = j d*x (- V -j-A)=J d'x {(d A A) - G) 2 - j ■ A^j 
(2-13) 

2 Remember that the dual of F is — F ! 
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with respect to the gauge field A^, leads to the equation (|2.12|) . The action 
(12.13)) is invariant with respect to the gauge transformation A — > A + (da) 
provided that d ■ j — 0. 

The source term G has to satisfy two conditions. The first is the equation 
9 ■ G — j mag . The second is that d ■ G ^ 0. If the second condition is not 
satisfied, the magnetic current decouples from the system. This is the reason 
why G cannot simply be expressed as G = d A B, in terms of another gauge 
potential B^. 3 String solutions (see Sect J2.4|) of the equation d ■ G — j mag 
are constructed in order to satisfy the condition 9-G*^0. 

The string term G^ u can be expressed in terms of two vectors, which we 
call E st and H st : 

Hl t = -G° l = ^e 0ljk G jk = E\ t = -G 0i = -^e 0ijk G jk (2.14) 

The equation d ■ F = d ■ G = j mag then translates to: 

V • H st = Pmag - d t H st + V x E st = J mag (2.15) 

Let us express the electric and magnetic fields E and H in terms of the 
vector potential and the string term. We define: 

A»=^A) (2.16) 

When the field tensor F has the form (j2.1(J|) . the electric and magnetic fields 
can be obtained from (|2.3|) . with the result: 

E = -d t A-V4> + E st H = VxA + H st (2.17) 

and we have: 

- l -F 2 = - l - (d A A - G) 2 = ~ {-d t A-V X ~ E st ) 2 -^(vxI+4) 2 

(2.18) 

• Exercise: Consider the following expression of the field tensor F : 

F = d AA-d7^B (2.19) 

3 The Zwanziger formalism, discussed in Sect 13. Ill does in fact make use of two gauge 
potentials. 
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in terms of two potentials N 1 and B**. Show that F will satisfy the 
Maxwell equations d ■ F = j and d ■ F = j mag provided that the two 
potentials A and B satisfy the equations: 

d-(dAA)=j d ■ (d A B) = j mag (2.20) 

Check that the variation of the action: 

I jJmag (A, B) = J d A x (-±F 2 -J-A + j mag ■ Bj (2.21) 

with respect to A and B leads to the correct Maxwell equations. What 
is wrong with this suggestion? A possible expression of the field tensor 
in terms of two potentials is given in a beautiful 1971 paper of Zwanziger 
@H1 (see Sect ED}- 

2.3 The current and world line 
of a charged particle. 

When we describe the trajectory of a point particle in terms of a time- 
dependent position R (t) , the Lorentz covariance is not explicit because t 
and R are different components of a Lorentz 4- vector. The function R (t) 
describes the trajectory in 3-dimensional euclidean space. Lorentz covari- 
ance can be made explicit if we embed the trajectory in a 4-dimensional 
Minkowski space, where it is described by a world line (r), which is a 
4-vector parametrized by a scalar parameter r. The parameter r may, but 
needs not, be chosen to be the proper-time of the particle. This is how Dirac 
describes trajectories of magnetic monopoles in his 1948 paper and much of 
the subsequent work is cast in this language, which we briefly sketch below. 
Let (r) be the world line of a particle in Minkowski space. A point 

r on the world line Z^ (r) = [t (r) , R{t) \ indicates the position R(r) of 

the particle at the time T (r) ,as illustrated in Fi d2.ll The current (x) 
produced by a point particle with a magnetic charge g can be written in the 
form of a line integral: 

f (x)=g J dZ>*5 4 (x - Z) (2.22) 
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along the world line of the particle. A more explicit form of the current is: 

9 I {x - Z {t)) (2.23) 



f (x) 



dr 



where To and t\ denote the extremities of the world line, which can, but need 
not, extend to infinity. 



z*Yc) = Cr(T),R(T)) 



Figure 2.1: The world line of a particle. For any value of r, the 4- vector 
Z^ (r) indicates the position R (r) of the particle at the time T (r). 

In order to exhibit the content of the current ()2.23j) , we express it in terms 
of a density p and a current j: 



(2.24) 



Let x M = (t,r). The current ()2.23|) . at the position r and at the time t, has 
the more explicit form: 



f Tl dZ^ 

f{t,r)=g / dr- J -5(t-T(r))5 3 (r- J R(r) 

J TO 



dr 



^2.25) 



The expression (I2.23|) of the current is independent of the parametrization 
Z M (r) which is chosen to describe the world line. We can choose r = T. The 
density p(t,r) is then: 

p(t,r)=3°(t,r)=g / dr—8 (t - T (r)) 5 (r - R (r 
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= g f 1 dr5(t-r)5 (f- R(t) \ = gS(f-M(t)\ (2.26) 
and the current j (t, r) is: 

j(t, r) = g jT dr^5 (t-r)5 (f-R (r)) = g™6 (V- R (*)) (2.27) 

The expressions ()2.26|) and ()2.27|) are the familiar expressions of the density 
and current produced by a point particle with magnetic charge g. 

2.4 The world sheet swept out by a Dirac 
string in Minkowski space 

The Dirac string, which is added to the field tensor F^ v in the expression 
(I2.1UJI . is an antisymmetric tensor G^ v (x) which satisfies the equation: 

d-G=j (2.28) 

As stated above, not any solution of this equation will do. For example, 
if we attempted to express the string term in terms of a potential B^ by 
writing, for example, G = d A B, we would have d ■ G = and the string 
term would decouple from the action ([2.130 . For this reason, string solutions 
of the equation d ■ G = j mag have been proposed. 

The string solution can be expressed as a surface integral over a world- 
sheet Z"(r,s): 

(x)=g [ drds d [ Z ^ Z ^ S 4 (x - Z) (2.29) 
J d{s,T) 

The world sheet Z^ (r, s) is parametrized by two scalar parameters r and s 
and: 



d{z^z v ) = dz^dz^ _ dz^ dz v 

d (s, t) ds dr dr ds 



(2.30) 



is the Jacobian of the parametrization. A point (r, s) on the world sheet 
Z^ (r, s) indicates the position R (r, s) at the time T (r, s) of a particle on 
the world sheet. The expression (J2.29[) for the string G is independent of 
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the parametrization of the world sheet Z M (r, s) and it can be written in a 
compact form as a surface integral over the world sheet Z: 

Gfiv (x) = g / da^ 8 (x - Z) (2.31) 
Js 

The surface element is: 




(2.32) 



Figure 2.2: The world sheet Z^ (s, r) swept out by a Dirac string, which 
stems from a particle with magnetic charge g and terminates on a particle 
with magnetic charge —g. The world line of the positively charged particle 
is the segment joining the points (sq, r ) and (so,Ti). The world line of the 
negatively charged particle is the joins the points (si, r ) and (s\, Ti). 

Figure 12.21 is an illustration of the world sheet which is swept out by a 
Dirac string which stems from a particle with magnetic charge g and termi- 
nates on a particle with magnetic charge —g. The word line of the positively 
charged particle is the border of the world sheet extending from the point 
(so,t ) to the point (so,Ti). The world line of the negatively charged parti- 
cle is the border of the world sheet extending from the point (si,r ) to the 
point (si,Ti). For any value of r, we can view the string as a line on the 
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world sheet, which stems from the point (sq, t) on the world line of the posi- 
tively charged particle to the point (s 1; r) of the world line of the negatively 
charged particle (see Fig |2.1|) . Often authors choose to work with a world 
sheet with a minimal surface. This is equivalent to the use of straight line 
Dirac strings. An observable, which is independent of the shape of the Dirac 
string, is independent of the shape of the surface which defines the world 
sheet. If the system is composed of a single magnetic monopole, that is, of 
a single particle with magnetic charge g, then the attached string extends 
to infinity. The corresponding world sheet has s\ — > oo and it becomes an 
infinite surface. 

Let us check that the string form (|2.29|) satisfies the equation d ■ G = j: 

dTds (^_^_a^a^\a s{x _ z) 



ds dr dr ds J dZ 
We have: 

f) f)7 a f) 

SF^-^TFaz^*-*) (2 ' 34) 

and a similar expression holds for (x — Z). We obtain thus: 

f / d7 fl d d7 fI d \ 

d a (F» (x) = -g J drds (x-Z)-^-jj-6(x- Z)\ (2.35) 

We can use Stoke's theorem which states that, for any two functions U (r, s) 
and V (r, s), defined on the world sheet, we have: 

fdUdV 3UdV\ I TT fdV 1 dV , \ , n „„. 

o — n ^r^r- ) = 4> u \ ^ ds + ^ dr ( 2 - 36 ) 



\ dr ds ds dr J Jc \ 9s dr 

where the line integral is taken along the closed line C which borders the 
surface S. A more compact form of Stoke's theorem is: 

d(U,V) 



d(r,s) 



UdV (2.37) 



c 



We apply the theorem to the functions U = 5 4 (x — Z) and V = Z^ 1 so as to 
obtain: 



'dZ» dZ^ 

d a G a » (x)=g<f> 5{x-Z)[ —dr + —ds ) = g <f> dZ^b (x - Z) 



c \ Ul uo J Jc 



^2.38) 
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We can choose, for example, the world sheet to be such that the path C 
begins at the point (s , r ) and passes successively through the points (s , Ti), 
(si,ri), (si,To) before returning to the point (so,r ). Then if the world line 
of the charged particle begins at (s , To) and ends at (s , 7~i), the other points 
being at infinity, the expression ()2.38|) reduces to: 

d a G afl (x)=g / dr^—5 4 (x-Z (r)) (2.39) 

Jro dT 

which, in view of (|2.23)) . is the current j M (x) produced by the magnetically 
charged particle. 



2.5 The Dirac string joining equal and oppo- 
site magnetic charges 

The string term G^ u (x) can be expressed in terms of the two vectors H st and 
E st defined in (12.14)1 . We can use ()2.29)) to obtain an explicit expression for 
these vectors. Thus: 



Hl t (*> r) = Gi ° (*> r)=9 J drd ^ (t -T)6 3 (r-R 



dR, dT dR< dT 



ds dr dr ds 



For a given value of s, we can choose r = T (r, s) in which case we have 
— 1 and I 1 = 0. The string term reduces to: 

Or as ° st 

H st (t,r) = g J ds^5 (f- = g J dR5 (r - R(t)} (2.40) 

The expression for i*^ is: 

^ (* , = = jf rfrrfs 6(t — T) 5 3 (f - i?) 

(2.41) 

so that: 

E st (t,?)=g J ds^ x ^ 3 (V-H(t,<r)) = g J dRx^d (r-R(t) 

(2.42) 
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— * — *■ 

The string terms H st and E st satisfy the equations (|2.15|) . Let us calculate 

V r • H st {t, r)=g J dR- V r 5 (r-R}=-g J dR- V R 5 (V- it) (2.43) 

Now, for any function / (i?) we have dR-V R f (r\ = f (R+SR\-f (it) . 

Let Ri (t) and R2 (t) be the points where the string L originates and termi- 
nates. We see that the expression ()2.43|) is equal to: 

V r • H st (t,r) = g5 (f- R x (*)) - g5 (r - R 2 (*)) (2.44) 

The right hand side is equal to the magnetic density of a magnetic charge g 
located at R\ and a magnetic charge —g located at R 2 - For such a system, 
we can choose a string which stems from the monopole g and terminates at 
the monopole —g. 



2.6 Dirac strings with a constant orientation 

Many calculations are made with the following solution to the equation d-G = 
j mag , namely: 

G = A j mag CT = ^ (/t 9 - nXj (2-45) 

where n M is a given fixed vector and n ■ d = n tl d fl . We can check that this 
form also satisfies the equation d • G — j mag - 

d a G ^ — — g^a [n j mag ~ n ^3mag) ~ jmag (2.46) 

where we assumed that the current j mag is conserved: d^j^ ag = 0. The 
solution ()2.45|) is used in many applications because it is simple and we shall 
call it a straight line string. 

Let us choose n M to be space-like: 

r^ = (0,n) n-d = n-V (2.47) 

The string terms E st and H st , defined in (I2.14J) are then: 
-* 1 -* -* n 

E s t — ~— —n x jmag-, H s t = — ~Pmag (2.48) 

n ■ V n ■ V 
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• Exercise: Use flA.80|) to check explicitly that the form (|2.48|) of the 

string terms satisfies the form (|2.15|) of the equation d • G — j mag - 

Consider first the case of a single magnetic monopole sitting at the point 
Ri. The monopole is described by the following magnetic current j mag : 

jmag = \PmagjjmagJ Pmag = 9$ yj" -Rlj jmag = (2.49) 

The equations (|2.48|) show that: 

V ■ H st = pma 9 = gS (r - Rl) E st = (2.50) 

The Fourier transform of H st is: 

H st (k)=g I d 3 r e* f -^5 (r- R x ) = ig^e illil (2.51) 
^ ' J n • V ^ ' n ■ k 

Let us choose the z-axis parallel to n so that tr=> = G( z )y where e*( 2 ) is a unit 

vector pointing in the z direction. We then have H st (k^j = ige^j-e 1 ^'^ 1 . 
The inverse Fourier transform is: 

H st (f) = 4orA*) I d 3 ke- fk (^)±- = ^3e (2) / £k e -*(^) j°° dz'e^' 

(2.52) 



(2tt) j W J k z (2tt) 

Let us define a vector R (z') = (X\, Yi,Zi + z'). We have: 

k- (r-Ri^j -k z z' = k- (r-R(z'f) dR(z') = e {z) dz (2.53) 



so that: 

"OO 



H st (r) = e (z) g dz 5 (f - R (z')) =g J dRS(f-Rj (2.54) 

where the path L starts at the point Ri and runs to infinity parallel to the 
positive ^-axis. 

Thus, when the density p (r) represents a single monopole at the point 
R\, the straight line solution (|2.48|) is identical to a Dirac string which stems 
from the monopole and continues to infinity in a straight line parallel to the 
vector n. If the system consists of two equal and opposite magnetic charges, 
located respectively at positions Ri and R2, then the straight line solution 
f!2.48|) represents two strings emanating from the charges and running to 
infinity parallel to the z-axis. Straight line strings, such as ()2.48|) with a 
fixed vector used in the Zwanziger formalism discussed in Sect 13. Ill 
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2.7 The vector potential A in the vicinity of 
a magnetic monopole 

Let us calculate the vector potential in the presence of the magnetic monopole. 
Since j = 0, the equation (j2.12|) for the vector potential A 11 is: 

d-(dAA)-d-G = (2.55) 

Let us write: 

A"=U,A\ (2.56) 
The equation ()2.55j) can then be broken down to: 

V • (-V0 - d t A} = 

d t (d t A + V</>) + Vx(vxl)=VxF sl (2.57) 

For a static monopole, it is natural to seek a static (time-independent) solu- 
tion. We can choose = 0. The equation for A reduces to: 

V x (V x Aj = V x H st (2.58) 

Let us distinguish the longitudinal and transverse parts of the vector poten- 
tial, respectively A L and A T : 

A L = ^V(V -A) A T = A-^W (V • A) Vxi L = V-A T = 

(2.59) 

The equation ()2.58|) determines only the transverse part At of the vec- 
tor potential A because V x A = V x At- It leaves the longitudinal 
part undetermined. Since the transverse part has V • At = 0, we have 

Vx^Vx At^ = —V 2 A T . Substituting for H st the string (|2.54|h the expres- 
sion for At becomes: 

-V 2 1 T = #V x J dRS (V- R) (2.60) 

At this point, a useful trick consists in using the identity ()2.69|) to rewrite 
the 5-function. We obtain thus: 



V 2 Ar = —V x / dRV 
4tt 



9 - f 1 



f- R 



(2.61) 
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so that we can take: 

9 



At f = —V x / dR 

An 



L 



r- R 



(2.62) 



In a gauge such that V • A = 0, the expression above becomes an expression 
for A, but no matter. The expression ()2.62|) gives the vector potential At 
for a Dirac string defined by the path L. 

An analytic expression for At may be obtained when the path L is a 
straight line, as, for example the straight line string ()2.54j) which runs along 
the positive z-axis. In this case the expression (|2.62|) reads: 

a C°° 1 

A T (f) = fVx e {z) / dz' - (2.63) 

47r J ° + (z - z'f 

In cylindrical coordinates f Appendix IA. 6. 2j ). At can be expressed in the form: 

A T {f) =e {e) A{p,z) (2.64) 

and ()A.105|) shows that this form is consistent with V ■ At = 0. Using again 
flA.106j) we obtain: 

After performing the derivative with respect to p, the integral over z' becomes 
analytic and we obtain the vector potential in the form: 

i T („,M)- ( „Al( 1 + _^_) (,66) 

In spherical coordinates ( Appendix I A . 6 . 3|) . the vector potential has the form: 

A T (p,9,z)=e {ip ° 1 -±^- (2.67) 
47r r sin 6 

We shall see in Sect. 14. II that the abelian gluon field acquires such a form in 
the vicinity of points where gauge fixing becomes undetermined. 
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Exercise: Use (JA.105|) to check that (|2.66|) is transverse: V • A T = 0. 
Use f|A.106|) to calculate the magnetic field from (|2.66|) . Check that, 
at all points which are not on the positive z-axis (where H st = 0), the 
magnetic field is equal to: 



H 



g r 



4wr 2 r 



(2.68) 



Exercise: Calculate the Coulomb potential produced by a charge sit- 
uated at the point f=R and deduce the identity: 

1 



5 r-R 



-V 2 



47r r- R 



[2.Q9) 



2.8 The irrelevance of the shape of the Dirac 
string 

Let us calculate the electric and magnetic fields E and H generated by the 
monopole. They are given by the field tensor (|2.1|) . There are two ways 
to calculate them. The complicated, although instructive, way consists in 
starting from the expression (j2.1()j) of the field tensor in terms of the vector 
potential and the string term. The dual string term G is then: 



El 



-G 



Qi 



HI 



t oijk G 



If we use (|2.1|) , the electric and magnetic fields become: 

E = H = VxA + H st 
Now V x A can be calculated from ()2.62j) : 



(2.70) 



(2.71) 



V r x A = — — 
Air 



V r x V r x dR 



f- R 



We use V x ( V x a 



V [ V • a I — V 2 a to calculate: 



(2.72) 



V r X V r X dR- 



V r • dR- 



-dR Vi 



R 



R 



R 
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1 



+ 4icdR6 [f-R 



-V I V R ■ dR 

f- R 

Substituting back into the expression for V x A, we obtain: 

9 



[2.73) 



V r xi = - j-V- ~ 9 I ,IR S ( F ~ R 



Anr 2 



: e( r ) — H st 



; (r) 



r 

(2.74) 

where we used (12.54)1 . Substituting these results into (|2.71|) we find that the 
electric and magnetic fields are: 



E = 



H = V x A + H . 



.sf 



47rr 2 



; (r) 



(2.75) 



The fields (|2.75|) could, of course, also have been obtained by simply 
solving the Maxwell equations (|2.8J1 with the magnetic current (|2.49jl , without 
appealing to the Dirac string. We have calculated them the hard way in order 
to show that the string term, which breaks rotational invariance, does not 
contribute to the electric and magnetic fields. As a result, the trajectory 
of an electrically or magnetically charged particle, flying by, will not feel be 
Dirac string. However, we shall see that the string term can modify the 
phase of the wavefunction of, say, an electron flying by, and this effect leads 
to Dirac's charge quantization. 

• Exercise: Start from the action (J2.13J1 and derive an expression for 
the energy of the system. Show that is does not depend on the Dirac 
string term. 

• Exercise: Equation ()2.58|) states that Vx(vxl)=-Vx H st . 
What would we have missed if we had concluded that V x A = —H st ? 



2.9 Deformations of Dirac strings and charge 
quantization 

Although the electric and magnetic fields are independent of the string term, 
the vector potential is not. What happens to the vector potential if we 
deform the Dirac string? Let us show that a deformation of the Dirac string 
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is equivalent to a gauge transformation. This is, of course, why the magnetic 
field is not affected by the string term. 

The vector potential is given by (|2.62|) . Let us deform only a segment of 
the path, situated between two points A and B on the string. The difference 
5 A (r) between the vector potentials, calculated with the two different paths, 
is the contour integral: 

5 Air) = -— V x / dR ' 



c 



f-R 



[2.7Q) 



where the contour C follows the initial path from A to B and then continues 
back from B to A along the deformed path, as shown on Fig )2.3l Using the 
identity (|A.94[) . we can transform the contour integral into an integral over 
a surface S whose boundary is the path C: 



5A(r) = x / ds x V R ' 

47r Js 



f-R 





1 \ 


( Vx / S d? - 




■ 




f- R 





^2.77) 



Using the identity Vx Vxa = V ( V • a I — V 2 a, we obtain: 



5 A (r) 



~V I I ds-V 

47T 






1 \ 










f-R 


J 



g I ds 5 [f — R 



(2.78) 



The second term vanishes at any point f not on the surface and can be 
dropped. The first term is the gradient of the solid angle Q, subtended by 
the surface S, when viewed from the point f. 



^« = ^v[/ s d s --v H 1 



f-R 



47T 



ds-if—R 



1 



f- R 



-- ^-Vfi (f) 



To see why, consider first a very small surface ds, such that 
essentially constant. Then 5s • if — R) , = i — ^-777 ( 5s 



f-R 

r~R 
\r-R\ 



[2.79) 
remains 

= 50 s . 
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Figure 2.3: The effect on the vector potential of deforming a Dirac string. 
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A finite surface S can be decomposed into small surfaces bounded by small 
contours which overlap (and therefore cancel each other) everywhere except 
on the boundary of the surface, that is, on the path C. The result (12.79)) 
follows. 

The expression (|2.79j) shows that 5 A is a gradient. The deformation 
of the string therefore adds a gradient to the vector potential A and this 
corresponds to a gauge transformation. A deformation of the Dirac string 
can therefore be compensated by a gauge transformation. 

This is, however, only true at points which do not lie on the surface S. 
Indeed, the solid angle Q (f) is a discontinuous function of r. The vector 
r — R changes sign as the point f crosses the surface. If the point r lies 
close to and on one side of the surface S (the shaded area in Fig. I2.3)) . the 
solid angle Q (f) is equal to 2tt (half a sphere). As soon as point r crosses 
the surface, the solid angle switches to — 27r. Thus the solid angle Q (f) 
undergoes a discontinuous variation of Am as the point crosses the surface 
S. As a result, the gauge transformation which compensates a deformation 
of the Dirac string is a singular gauge transformation. This point will be 
further discussed in Sect. 13.3.21 

Consider the wavefunction ip (r) of an electron. (We consider an electron 
because it is a particle with the smallest observed electric charge.) When 
the vector potential undergoes a gauge transformation A — > A — J^Vfi, the 
electron wavefunction undergoes the gauge transformation ip — > e ie ^ip. This 
means that, on either side of the surface S, the electron wavefunction differs 
by a phase e ie ^ = e ieg . This would make the Dirac string observable, unless 
we impose the condition: 

eg = 2mr (2.80) 

where n is an integer. The expression ()2.8()|) is the charge quantization con- 
dition proposed by Dirac. In his own words, "the mere existence of one 
[magnetic] pole of strength g would require all electric charges to be quan- 
tized in units of 2im/g and similarly, the existence of one [electric] charge 
would require all [magnetic] poles to quantized. The quantization of electric- 
ity is one of the most fundamental and striking features of atomic physics, 
and there seems to be no explanation of it apart from the theory of poles. 
This provides some grounds for believing in the existence of these poles" 4 
j^j. This was written in 1948. Dirac used a different argument to prove 

4 In the 1948 paper of Dirac, the quantization condition is stated as eg — hnHc. Ex- 
cepting the use of units h = c = 1, the charge e defined by Dirac is equal to Att times the 
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the quantization rule (J2.8U|) . The proof given above is taken from Chap. 6 
of Jackson's Classical Electrodynamics jUj. There have been many other 
derivations [1E1 EH]- For a late 2002 reflection of Jackiw on the subject, see 
reference j50j . 

The discontinuity of the solid angle fl (r) implies that, on the surface S, 
we have V x ^Vfij 7^ 0, so that, for example, (d x d y — d y d x ) Vt 7^ 0. To 

see this, consider the line integral § L (x/Q^j ■ dl taken along a path L which 

crosses the surface S (the shaded area on Fig. 12. 3J) . The integral is, of course, 
equal to the discontinuity of Q across the surface S. It is therefore equal to 
47r. However, in view of Stoke's theorem (jA.93l we have: 

j (vfi) ■ dl = J ds- (V x (VO)) = 4tt (2.81) 

— > 

where S (L) is the surface bounded by the path L. It follows that V x 
^Vfij 7^ on the surface S. 



2.10 The way Dirac originally argued for the 
string 

In his 1948 paper Dirac had an elegant way of conceiving the string 
term in order to accommodate magnetic monopoles. He first noted that the 
relation F^ u = d^A u — d u A^ implied the absence of magnetic charges. In an 
attempt to preserve this relation as far as possible, he argued as follows. If 
the field tensor had the form F^ u = d^A u — d u A^, then the magnetic field 
would be H — V x A. The flux of the magnetic field through any closed 
surface S would then vanish because of the divergence theorem (jA.90|) : 

J H.ds = J [Vxi) .ds = J d 3 r V. (vxl)=0 (2.82) 

where V is the volume enclosed by the surface 5* and ds a surface element 
directed outward normal to the surface. However, if the surface S encloses 

charge we use in this paper. For example, Dirac writes the Maxwell equation d ■ F = 4irj 
whereas we write it as d ■ F = j. The quotation of Dirac's paper is modified so as to take 
this difference into account. 
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— * — * 

a magnetic monopole of charge g, then the Maxwell equation V .H = g5 (r) 
states that the total magnetic flux crossing the surface S should equal the 
magnetic charge g of the monopole: 

/ H.ds= [ d 3 rV.H= [ d 3 rgS(r)=g (2.83) 
Js Jv Jv 

Dirac concluded that "the equation = d^A u — d u A^ must then fail some- 
where on the surface S" and he assumed that it fails at only one point on the 
surface S. "The equation will then fail at one point on every closed surface 
surrounding the magnetic monopole, so that it will fail on a line of points" , 
which he called a string. "The string may be any curved line, extending from 
the pole to infinity or ending at another monopole of equal and opposite 
strength. Every magnetic monopole must be at the end of such a string." 
Dirac went on to show that the strings are unphysical variables which do not 
influence physical phenomena and that they must not pass through electric 
charges. He therefore replaced the expression H = V x A, by the modified 
expression: 

H = V x A + H st (2.84) 

• Exercise: Try the formulate the theory in terms of two strings, each 
one carrying a fraction of the flux of the magnetic field. Under what 
conditions can a monopole be attached to two strings? 

2.11 Electromagnetism expressed in terms of 
the gauge field associated to the dual 
field tensor F^ v 

The scheme developed in Sect. 12.21 is useful when a model action is expressed 
in terms of the vector potential A^ related to the field tensor F = d A A — G 
and when magnetic charges and currents are present. We shall however be 
interested in the Landau- Ginzburg model of a dual superconductor, which 
is expressed in terms of the potential associated to the dual tensor F^ v 
and in which (color) electric charges are present. We cannot write the dual 
field tensor in the form F = d A B because that would imply that d ■ F = 
— <9 • d A B — 0, which would preclude the existence of electric charges. The 
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way out, of course, is to modify the expression of F^ v by adding a string 
term G^ u and writing the dual field tensor in the form: 

F = dAB + G F = -~dA~B + G (2.85) 

We require the string term G to be independent of B and to satisfy the 
equation: 

d ■ G = j (2.86) 

Note that, in the expression ()2.85|) . the string term G is added to the field 
tensor F, whereas, in the Dirac formulation (j2.10j) . the string term is added 
to the field tensor F. The roles of F and F are indeed interchanged when 
we express electromagnetism in terms of the gauge field 

This way, the first Maxwell equation d ■ F = d ■ G = j is satisfied inde- 
pendently of the field The latter is determined by the second Maxwell 
equation d • F — j ma g, namely: 

d ■ (d AB) + d ■ G = j mag (2.87) 

where j^ ag is a magnetic current, which in the dual Landau-Ginzburg model, 
is provided by a gauged complex scalar field. The equation for B^ may be 
obtained from the variation of the action: 

hu a9 (B) = J d A x ^-F 2 - j mag ■ B^j = J d 4 x (~ (d A B + G) 2 - Jmag ■ B 

(2.88) 

with respect to the gauge field B^ . The action ()2.88j) is invariant under the 
gauge transformation B — >• B + {d(3) provided that d ■ j mag = 0. 

The source term G has to satisfy two conditions. The first is the equation 
d-G = j. The second is: d- G ^ 0. Otherwise, the electric current decouples 
from the system. String solutions satisfy the second condition, whereas a 
form, such as G = d A A does not. 

We define: 

B» = (B, x) d "=Q^= ( d *> -^) f = (a J) (2-89) 

We can express the string term G^ and its dual G^ v in terms of two euclidean 
3- vectors E st and H st in the same way as the field tensor F^ v is expressed in 
terms of the electric and magnetic fields. In analogy with (J2.3)) we define: 

E\ t = -G 0i = ^e 0ijk G jk H { st = ~G 0i = ~\e 0ijk G jk (2.90) 
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Be careful not to confuse these definitions with the definitions ([2.14)1 ! 

With a field tensor of the form (|2.85|) , the electric and magnetic fields can 
be obtained from ()2.3j) with the result: 

E = -VxB + E st H — -d t B - Vx + H st (2.91) 

The equation d ■ G = j translates to: 

V ■ E 8t = P - d t E st + VxH st =j (2.92) 

and we have: 

- l -F 2 = -~ (d A B + G) 2 = ~ (-dtB - V X + H s tf- 1 - (-V x B + £ st ) 2 

(2.93) 

The source term E st has to satisfy two conditions. The first is the equation 
V • E st = p. The second is that V x E st ^ 0. 

The magnetic charge density and current are given by (|2.8J) : 

Pmag = V • H = V ■ (-5 t S - Vx + # st ) 

J ma9 = -d t H -VxE = -dt {-dtB - Vx + H*) - V x (-V x5 + E^j 

(2.94) 

Consider the case where the system consists of a static point electric 
charge e at the position Ri and a static electric charge — e at the position 
i?2- The charge density is then: 

p(f) =eS{f-R^ -e8(f-R 2 ) (2.95) 

and the electric current j vanishes. In that case, we can use a string term of 
the form: 

E st {f f ) = ejdR5{r- R) H st (r) = (2.96) 

where the line integral follows a path L (which is the string), which stems 
from the charge e at the point R\ and terminates at the charge — e at the 
point i?2- A point on the path L can be parametrized by a function R (s) 
such that: 

dR 

R 1 = R(s 1 ) R 2 = R(s 2 ) dR = ds— (2.97) 

ds 
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in which case the line integral (|2.96|) acquires the more explicit form: 

r S2 dR / \ 
E st (r) = e J ds— 5 [r- R(s)j (2.98) 

The argument which follows equation (J2.43)) can be repeated here with the 
result: 

V r -E st = e J dR-V r 5(r-R)=-e J dR-V R S (r - R) = eS (r- R^-eS 

' L ' L (2.99) 

If we had a single electric charge e at the point Ri the string would extend 
out to infinity. 

Many calculations are performed with straight line strings, discussed in 
Sect. \2.ti\ and which is the following solution of the equation d ■ G = j: 

G = -^—n A j G» u = «f - n v f) (2.100) 
n • a n ■ a 

where is a given fixed vector and n ■ d = n^d^. This form solves the 
equation d ■ G — j if <9 M j M = 0, that is, if the electric current is conserved. If 
we choose n 11 to be space-like, the string terms E st and H st are given by: 

-» -» n -* I -* 

n" = (0, n) , n ■ d = n ■ V, E st = -^^rp, H st = -^^n x j 

n ■ V n ■ V 

(2.101) 



Consider the Fourier transform of the source term E xi : 

E st (k) = I d 3 r e^-^p (f) = — ^ p[k) (2.102) 



n n 
_—^p{r) = — — 
n • V in - k 

The inverse Fourier transform is: 

Est (r) = -7^-3 / d 3 fc e^A / dh'e~ ik f p {?) (2.103) 



(27r) y m • A; 

Let us choose the z-axis to be parallel to n, so that — = k z and - = e 



*(*) 



where et z \ is the unit vector pointing in the ^-direction. We obtain: 

E st (r) = -e (z) -^3 J d 3 ke~^^- J d z r' e~ itf p{¥) (2.104) 
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We obtain: 



E st (r) = e (z) j d 3 r'5 (x - x') 6 (y - y') 8(z- z') p (f) = e (z) j dz' p (x, y,z- z') 

(2.105) 

Let 

R ( z ') = (0, 0, z) dR = e (z) dz f - R= (x,y,z - z) (2.106) 
We obtain: 

E s t (r) = dz'^^- p(r-R (/)) = J^dRpfr- if) (2.107) 

where the path L is a straight line, starting at the origin and running parallel 
to the z-axis. The expression ()2.107|) provides for a determination of the 
operator used in ()2.101|) . For example, if the system has a single charge 

e at the point Ri, the density is p (r) = e5 (r ' — Rij and the expression yields 

a string term E st (r) which stems from the point R\ and extends to infinity 
parallel to the z-axis. If there is an additional charge — e located at the point 
i?2, then the expression will yield an additional parallel string, stemming 
from the charge — e and extending to infinity parallel to the z-axis. It will 
not be a string joining the two charges, unless the two charges happen to 
be located along the z-axis. Of course, if the strings can be deformed so as 
to merge at some point, then the two strings become equivalent to a single 
string joining the equal and opposite charges. 
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Chapter 3 



The Landau- Ginzburg model of 
a dual superconductor 

We shall describe color confinement in the QCD ground state in terms of 
a dual superconductor, which differs from usual metallic superconductors in 
that the roles of the electric and magnetic fields are exchanged. The dual 
superconductor will be described in terms of a suitably adapted Landau- 
Ginzburg model of superconductivity. The original model was developed in 
1950 by Ginzburg and Landau [H]. Particle physicists refer to it today as the 
Dual Abelian Higgs model. The crucial property of the dual superconductor 
will be the Meissner effect |3T] , which expels the electric field (instead of the 
magnetic field, as in a usual superconductors). As a result, the color-electric 
field which is produced, for example, by a quark-antiquark pair embedded 
in the dual superconductor, acquires the shape of a color flux tube, thereby 
generating an asymptotically linear confining potential. 

It is easy to formulate a model in which the Meissner effect applies to the 
electric field. All we need to do is to formulate the Landau-Ginzburg theory 
in terms of a vector potential associated to the dual field tensor , as 
in Sect 12 .111 For an early review and a historical background, see the 1975 
paper of Jevicki and Senjanovic jHj- The presentation given below owes a lot 
to the illuminating account of superconductivity given in Sect. 21. 6 of vol.2 of 
Steven Weinberg's Quantum Theory of Fields jH2j. We first study the dual 
Landau-Ginzburg model with no reference to the color degrees of freedom. 
The way the latter are incorporated is discussed in ChapEl 
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3.1 The Landau- Ginzburg action of a dual 
superconductor 



The Landau-Ginzburg (Abelian Higgs) model is expressed in terms of a 
gauged complex scalar field xp, which, presumably, represents a magnetic 
charge condensate. The model action is: 

i s {b^) = J d A x (-If^+1(d^) (D^y-h^r-v 2 ) 2 ^ 

(3.1) 

where ip is the complex scalar field and F^ v the dual field tensor. The 
covariant derivative is _D M = <9 M + igB^, where B^ is a vector potential, and: 

(D^) = {d^ + igB^) (D^Y = (d,r - igB^*) (3.2) 

The dimensionless constant g can be viewed as a magnetic charge. As ex- 
plained in Sect 12. Ill the presence of electric charges and currents can be 
taken into account by adding a string term G to the dual field strength 
tensor, which is: 

W v = (8 A Bf v + GT = - (FA~B)^ + GT (3.3) 

The string term is related to the electric current by the equation: 

d«G°* = f (3.4) 

The Landau-Ginzburg action (|3.1j) is invariant under the abelian gauge trans- 
formation: 

n (x) = e- i9f3(x) -> QD^ 

B ft ^B li + (d^) i) -> e-^V (3-5) 

The last term of the action is a potential which drives the scalar field to a 
non- vanishing expectation value ipip* = v 2 in the ground state of the system. 
The superconducting phase occurs when ipip* = v 2 ^ 0, and it will model the 
color-confined phase of QCD. The normal phase occurs when ipip* = and it 
represents the perturbative phase of QCD. The model parameter v may be 
temperature and density dependent, and its variation can drive the system 
to the normal phase. Of course, other processes may also contribute to the 
phase transition. Note that, when ipip* ^ 0, the action ()3.8)1 is not invariant 
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under the gauge transformation B — > B + {0(3) of the field B^ alone because, 
loosely speaking, the gauge field B^ acquires a squared mass g 2 ipip* . In the 
compact notation described in App[XJ the action reads: 

Ij (B, if;, r) = Jd 4 x^(dAB + G) 2 + ± \dj> + igBif;] 2 - U - v 2 ) 2 

(3.6) 

The physical content of the model is often more transparent in a polar 
representation of the complex field if): 

if){x) = S (x) e i9 ^ x) if)* (x) = S (x) e~ i9 ^ x) (3.7) 

The Landau-Ginzburg action ()3.1|) can be expressed in terms of the real fields 
S and <p: 

7, (5, <p, S) = J d*x (^(9A5 + G) 2 + ^(B + d^f + \ (OS) 2 - \b (S 2 - v 2 ) 2 ^) 

(3.8) 

The action (|3.8|) is invariant under the gauge transformation: 

B^B + (df3) <p^<p-/3 S-^S (3.9) 



In the ground state of the system, S = v , and fluctuations of the scalar 
field S describe a scalar particle with a mass: 

m H = 2vVb (3.10) 

Particle physicists like to refer to S as a Higgs field and to mu as a Higgs mass. 
The field (p remains massless and is sometimes referred to as a Goldstone field. 
The gauge field develops a mass: 

my = gv (3-11) 

Properties of superconductors are often described in terms of a pen- 
etration depth A and a correlation length £, which are equal to the inverse 
vector and Higgs masses: 

A = — £=— (3.12) 
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In usual metallic superconductors, the penetration length is the distance 
within which an externally applied magnetic field disappears inside the su- 
perconductor. In our dual superconductor, the penetration length A will 
measure the distance within which the electric field and the magnetic current 
vanish outside the flux-tube which develops, for example, between a quark 
and an antiquark. The correlation length is related to the distance within 
which the scalar field acquires its vacuum value S = v. It is also a measure 
of the energy difference, per unit volume, of the normal and superconducting 
phase, usually referred to as the bag constant: 



B = ^ m » v2 = ^ = We (3 - 13) 



In type I superconductors (pure metals except niobium) £ > A and my > run- 
in type II superconductors (alloys and niobium) A > £ and m# > my. In 
Sect. 13.41 we shall see that the dual superconductors which model the con- 
finement of color charge have m# ^ my. They are close to the boundary 
which separates type I and type II superconductors. The London limit (Sect. 
13. 6j) . in which it is assumed that b —>■ oo so that m# ^> my, is an extreme 
example of a type II superconductor. In type II superconductors, the only 
stable vortex lines are those with minimum flux. In type I superconductors, 
vortices attract each other whereas they repel each other in type II super- 
conductors. Useful reviews of these properties can be found in Chap. 21. 6 
(volume 2) of Weinberg's "Quantum Theory of Fields" [22] and in Chap. 4. 3 
of Vilenkin and Shellard's " Cosmic Strings and Other Topological Defects" 



3.2 The Landau-Ginzburg action in terms of 
euclidean fields 

Let us write: 

B»=(x,B) (3.14) 



and let us express the antisymmetric source term in terms of the two 
vectors E st and H st as in (|2.90|) . The action (|3.1|) can then be broken down 
to the form: 
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I ^ dt B - V X + H st y - l - (-V x B + E st y + X - (d t il> + ig X ip) W - i gX >J>*) 



(3.15) 



Since no time derivative acts on the field Xi it ac ts as the constraint j- = 0, 
namely: 

V • (-d t B - Vx + + y (W - Tp*d t ip) + g 2 Xipip* = (3.16) 

The Eg. (J3.4j) which relates the source terms to the electric charge density 
and current reads: 

V-E st = p - d t E st + V x H st = j (3.17) 

3.3 The flux tube joining two equal and op- 
posite electric charges 

Consider a system composed of two static equal and opposite electric charges 
±e placed on the z-axis at equal distances from the origin and separated by 
a distance R. The charge density is then: 

p(f^ = eS (r- R^-e6 (r- R^ R\ = \ ^ = U 0, | J 

(3.18) 

The electric current is then j^ 1 = <5 M °p with j = 0. As shown in Sect 12 .111 the 
string terms satisfy the equations: 

V-E st = p Vx£ s Vfl H st (r) = (3.19) 

Note the condition V x E st ^ 0. If this condition is not satisfied, the electric 
density decouples from the system, as can be seen on the expression ()3.22|) 
of the energy. String solutions are designed to avoid this. 
The string term E st has the form 1)2.96)1 : 

E st (f) = e f'dZSff-Z) (3.20) 
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where the integral follows a path L (the string), which stems from the point 
Ri and terminates at the point R 2 . Following the steps described in Sect. 

we can easily check that the form (J3.20j) satisfies the equation V • E s t — p 
with p given by (|3.18j) . 

When the fields are time-independent, the energy density is equal to 
minus the action density given by ()3. 15j) . The energy of the system is thus: 



(Fr 



1 \ 2 1 
"2 N + 2 



-V x B + E. 



A (y<i[> - igB^ (W>* + igBip*) + h (^* - v 2 ) 2 



st 



(3.21) 



The constraint ()3.1(ij) is satisfied with x = 0. The energy becomes the fol- 
lowing sum of positive terms: 

_ 1 



d r 



-V x B + E 



st 



+1 (yip - igB^j (y^ + igBip*) + h {^* - v 2 ) 2 



(3.22) 



This expression can also be derived from the classical energy (|3.168|) by 
making the energy stationary with respect to the conjugate momenta H and 
P, as becomes time-independent fields. 



3.3.1 The Ball-Caticha expression of the string term 



The string term (|3.2(J|) does not depend on the fields B,ip and if)*. A useful 
trick, introduced by Ball and Caticha [53| . and used in all subsequent work, 
consists in expressing the string term E st in terms of the electric field E° 
and the dual vector potential B°, which are produced by the electric charges 
when they are embedded in the normal vacuum where ■0 = 0. It is a way to 
express the longitudinal and transverse parts of the string term E st in terms 
of the known fields E° and B°. 

The electric field E° produced by the electric charges embedded in the 
normal vacuum is the well known Coulomb field: 



1 



Ri 



Ro 



e 

Air 



Ri 
Ri 



-R2 
Ri 



3.23) 
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This electric field E° can, however, also be expressed in terms of the dual 
potential B° and the string term E st , using the expression (|2.91|) : 



E° = -V x B° + E st 



(3.24) 



The idea is to use this equation in order to express the string term E st in 
terms of E° and B°. 

We can calculate the dual vector potential B° as in Sect. 12.71 Since no 
magnetic current j mag occurs in the normal vacuum, the field B° is given by 
(jHH) : 

V x (v x B°) = V x E st (3.25) 



(3.26) 



The string term E st is given by the line integral ()3.20|) so that: 



VxfVx B") = rV, •: I ciZ S { r — Z 



This is an equation for the transverse part B^ which is the only part we 
need. We have V x (v x B° J = — V 2 £?° and we can use ()2.69j) to write the 
equation above in the form: 



V?fl° = -— V r x [ dZV 2 ' 

47T J L 



f — Z 



(3.27) 



so that we can take: 



B° (f) = — V r x / dZ 
Air 



f — Z 



(3.28) 



From ()3.24j) we see that we can write the string term in the form E st = 
£° + Vx B°. We substitute this expression into the energy ()3.22|) . with the 
result: 

E p (V, i)\ = J d 3 r l - (- V x B + E° + V x B° 



(3.29) 
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Since E° is a gradient, the mixed term E°- x ^— B + vanishes. The 

field E° contributes a simple Coulomb term to the energy: 



d 3 r- El 



AnR 



+ {terms independent of R) 



(3.30) 



In the following, we neglect the (albeit infinite) self-energy terms which are 
independent of R. The energy can thus be written in the form: 



AtiR 



+ / d 3 r 



- (- V xB + VxB° U 



where B° is given by f)3.28j) . 



(3.31) 



3.3.2 Deformations of the string and charge quantiza- 
tion 

Consider the effect of deforming the string, that is, the path L which defines 
the vector B° (r) in the expression ()3.28j) . Let us deform a segment of the 
path, situated between two points A and B on the path. The corresponding 
change of B° is: 



B° (f) B° (r) + — V r x 

4:71 



1 



dZ 

< f-Z 



(3.32) 



where the contour C follows the initial path from A to B and then continues 
back from B to A along the modified path, as illustrated on Fig |2.3l The 
expression (J3.32)) is the same as the expression (J2.76|) and we can therefore 
repeat the argument given in Sect 12.91 to show that the deformation of the 
string L adds a gradient to B°: 



B° (r) -> B° (f) + —Vtt 
An 



(3.33) 



where Q is the solid angle subtended by the surface S, bounded by the contour 
C, viewed from the point f. The energy (|3.31j) is changed to: 



AttR 



+ / d 3 r 



Air 



-V x VQ 
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+ 



- (Vip - igB^j (Vip* + igBtp*} + -b (V>V>* - v 2 ) 2 



(3.34) 



We have purposely not set to zero the term V x |^Vf2j because, as shown in 

Sect 021 the solid angle Q (f) is a discontinuous function of r. It undergoes 
a sudden change of 4n as r crosses the surface S bordered by the path C 

(the shaded area in Fig. The Eq. (gSD shows that V x (vfl) is 

non-vanishing on the surface S 1 . We can, however, compensate for the extra 

term x ^Vfij in the energy (|3.34|) by performing the following singular 

gauge transformation: 



B 



(3.35) 



which reduces the energy (I3.34j) to its original form ([3.31)1 . The energy (|3.31j) 
becomes thus independent of the shape of the string. The gauge transfor- 
mation ()3.35|) is not well defined because the transformed field e~ t9e ^ip is 
a discontinuous function of f thereby making the gradient Ve" ,9e ^t/) ill de- 
fined. However, we can impose the condition: 



eg 



2rm 



(3.36) 



which makes the field e~ l9e i^ifj a continuous and differentiate function of f. 
We recover the Dirac quantization condition ()2.8()jl . Thus, deformations of 
the string can be compensated by singular gauge transformations. 

3.3.3 The relation between the Dirac string 
and the flux tube in the unitary gauge 

It is convenient to express the energy (|3.31|) in terms of the polar represen- 
tation (|3.7J) of the complex scalar field: 



E p [B,cp,S 



AttR 



+ / <fV 



- (-V xB + VxB '" 



9 2 S 2 



(B-vv) 2 + \(ys)\\b{s 2 -v 2 ) 2 



(3.37) 
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As shown in Sect 13.3.21 a modification of the shape of the string, which 
defines B Q in Eq. (J3.28|) . adds the gradient ^:Vfi to B . The corresponding 
modification of the energy can be compensated by the gauge transformation 

B^B + —(vd] tp^tp + ige— (3.38) 
47r V J An 

Because the energy ()3.37j) is invariant under the gauge transformation: 
B^B-Vp cp^cp-fi (3.39) 

we can choose the gauge (3 = <p>, which is usually referred to as the unitary 
gauge. In this gauge, the field (p vanishes and the energy ()3.37|) is equal to: 



£ p (B,S)=-^ + / cfr 



- (-V xB + Vx B ^ 2 



(3.40) 



The energy (|3.40|) . expressed in the unitary gauge, is not independent of 
the shape of the string which defines the field B°, nor should it be, because 
modifications of the shape of the string are compensated by modifications 
of the phase (p. When, for example, flux tubes joining electric charges are 
calculated by minimizing the energy (|3.40j) expressed in the unitary gauge, 
the flux tubes follow and develop around the Dirac strings. For example, in 
Sect 13. 3. 51 we shall see that the string term represents the longitudinal part 
of the electric field. In the unitary gauge, the shapes of the Dirac strings can 
be chosen so as to minimize the energy. They can subsequently be deformed, 
by re-introducing the field ip, as in the expression ()3.37|) for example. 



3.3.4 The flux tube calculated in the unitary gauge 

When the system consists of two static equal and opposite charges, the charge 
density is given by ()3.18j) and the Dirac string is a straight line joining the two 
charges. The field B° (r) is given by the expression (j3.28|) with dZ = e^dz, 
where e*( 2 ) is a unit vector parallel to the z-axis. The explicit expression can 
be written in cylindrical coordinates (App lA.fT2|) : 

B° (r) = x e (z) / ' dz' (3.41) 

4vr J_r / 9 , / /\2 

~ VP 2 + {z-z') 
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Use (lA.K)fij) to get: 

B° (r) 



Air 



dz'- 



2N§ 



(f + (z - z') 2 ) 



( 



He) 



e 1 

4-7T p 



~ R 

Z 2 



Z+f 



\ 



v/^ + ^-f) 2 ) y / (p 2 + (- + f) 2 )J 



e {6) B° (p, z) 



(3.42) 

This is an analytic expression for the field B° in cylindrical coordinates. From 
(jA.105|) we can check directly that V.-B = so that B° is transverse. Near 
the z-axis, where p is small, the field B° (p, z) becomes singular: 



B° (p, z) 



e 1 / z- \R z+\R 



An p \ \z- \R\ |z+ \R\ 



when -\ R<Z <% 
when z < — ^ and z > |- 



(3.43) 



The fields which make the energy ()3.4()j) stationary satisfy the equations: 
V x ( V x b) - V x (v x B°) + g 2 S 2 B = 



-V 2 + 26 (S 2 -v 2 )+g 2 B 7 



S = 



(3.44) 



When B° has the form (|3.42|) . a solution exists, in cylindrical coordinates, in 
which the fields S and B have the form: 



S(r) = S (p, z) 



£ if) = Ho)B (p, z) 



(3.45) 



The field equations (|3.44j) reduce to the following set of coupled equations 
for the functions B (p, z) and S (p, z): 



d 2 (B - B 

dz 2 

1 d 



g 2 S 2 B = 



p—S 
pdp \ dp 



2b (S 
45 



2 v 2 )S + g 2 C 2 B = 



(3.46) 



Using (|A.106|) , we can express the energy (|3.4U|) in terms of the functions 

S (p, z) , B (p, z) and B° (p, z) : 



£ P (B,S 



+2ir 



pdp 



dz 



1 (d{B-B QS 
2 



dz 



d_ 

2p 2 \dp 



p(B 



g 2 s* 2 ifdsy ifdsy i /a2 



(3.47) 



±oo), the fields 



We require that, far from the sources (p — > oo, z 
should recover their ground state values S = v and B = 0. Close to the 
z-axis, the field B (p, z) f=a - thereby making the electric field finite in this 
region. However, such a beriavior would make the contribution of the term 
2-^-B 2 diverge, unless S — > in the vicinity of the string. This is reason why 
the energy (as well as the string tension), calculated in the London limit, has 
an ultraviolet divergence (see Sect J3.3.5j) . 

The role of the model parameters is made more explicit, if we work with 
the following dimensionless fields and distances: 

B(p,z)-B° (p,z) =vl(x,y) S(p,z) =vs(x,y) B° (p, z) = vb (x,y) 

x = gvp = m v p y = gvz = m v z (3.48) 

where m# and my are the Higgs and vector masses ({SUP} and l3.1ljl . respec- 
tively equal to the inverse penetration and correlation lengths (j3.12j) . The 
energy acquires the form: 



t-^2. 



AnR g 2 



+- 



xdx 



dy 



1 (dl_ 
2 



\^(^xl" 
2x 2 V dx 



ds 
dx 



2\dy) + 8m 2 v [ } 



(3.49) 



where = 

9 

The flux tube obtained by minimizing the energy (|3.49|) is calculated and 

The flux 



displayed in the 1990 paper of Maedan, Matsubara and Suzuki 
tube is very similar to the one obtained in lattice calculations, illustrated in 
Fig l3.ll A more detailed comparison to lattice data is made in Sect J3. 41 More 
recently, progress towards analytic forms for the solutions has been reported 
in the 1998 paper of Baker, Brambilla, Dosch and Vairo jHE] ■ 
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Figure 3.1: The lines of force of the electric field between two static color- 
electric SU (2) sources obtained from a lattice calculation in the maximal 
abelian projection [3*T] . 



3.3.5 The electric field and the magnetic current 

The electric and magnetic fields are given by ([2.91)1 : 

E = — V x B + E st H = (3.50) 

The longitudinal part of the electric field is the string term E st and it is given 
by: 

V • E = V • E st = p (3.51) 

A magnetic current (|2.94jl is produced by the transverse part of the electric 
field: 

Jmag = ~V X E (3.52) 

The relation j mag = —V x E is sometimes referred to as the "Ampere law". 

When the fields have the cylindrical symmetry ()3.45j) , the electric field E 
and the magnetic current j mag , given by ()3.5()j) and ()3.52j) . are: 

dB 1/9 \ 
E(p,z) = -VxB + WxB + E° = e (p) — - e {z) - 1^-pB) + E° 

Jmag (p, z) = V X X Bj - V X X S ) 
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_ fd 2 (B-B°) did, nx \ , 

= -^{- A 8^ + a- p - p ir/ {B - B) ) (3 ' 53) 

We see that the magnetic current circulates around the z-axis. This is why 
the flux tube is often called a vortex. 



3.3.6 The Abrikosov-Nielsen-Olesen vortex 

When myR 3> 1 and m^R ^> 1, that is, when the distance which separates 
the electric charges is much larger than the width of the flux tube, in regions 
of space where both p and z are much smaller than R, the source term (|3.42j) 
reduces to B° = — e^)^-, and this in turn implies that V x B° = 0. In this 
region of space, the fields ()H.45j) . which are solutions of the equations ()3.44|) . 
become independent of z and they acquire the simpler form: 

S(r) = S(p) B (f) = e w B (p) (3.54) 

The electric field (I3.53j) points in the z-direction: 

1 c) 

E (f) = -V x B = -e (z) -— ( P B) (3.55) 

pop 

and a magnetic current j mag circulates around the z-axis: 
-> -> -> did 

Jmag = "V x E = -e (e) — -—pB (3.56) 

A flux tube is formed which is referred to as the Abrikosov-Nielsen-Olesen 
vortex, studied in the 1973 paper of Nielsen and Olesen pQ. It is the ana- 
logue of the vortex lines, which were predicted to occur in superconductors 
by Abrikosov [HZj- Being particularly simple, the Abrikosov-Nielsen-Olesen 
vortex has been extensively studied. In Sect. 13.41 we shall see that lattice 
simulations confirm the formation of an Abrikosov-Nielsen-Olesen vortex be- 
tween equal and opposite charges. 

The energy density becomes independent of both z and 9, so that, in 
terms of the dimensionless fields and distances ()3.48|) . the energy per unit 
length along the z-axis reduces to: 

d£ R (b,s) 2 f°° , / 1/3 \ 2 1 2l2 1 fds\ 2 mjj , 2 , 

(3.57) 
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where b(x) = vB (z) and where v 2 = ^f. The fields s (x) and b (x) which 
make the energy stationary are the solutions of the equations: 



1 d d o, m 



-— x— s + s 2 b+ — %r is 2 - 1) s = (3.58) 
x dx dx Zmt ' 



The boundary conditions are: 



b -> - s -> 6^0 s -»• 1 (3.59) 

The constant a can be determined from the flux of the electric field crossing 
a surface normal to the vortex. 

The expression (|3.57|) is interpreted as the string tension, that is, the 
coefficient a of the asymptotically linear potential aR which develops between 
static electric charges embedded in the dual superconductor. The string 
tension depends on the two parameters v and the ratio of the Higgs 
and vector masses. When m# = my, that is, when the system is on the 
borderline between a type I and II superconductor, analytic solutions of the 
equations of motion have been found by Bogomolnyi [SB] • In this Bogomolnyi 
limit, the string tension is given by the expression: 



\fo=.\—m 2 m = ran = ray (3.60) 
V 9 

The stability and extensions to supersymmetry have also been investigated. 
For a review, see the Sect. 4.1 of the useful book by Vilenkin and Shellard 
jSH]. The interaction between vortices is discussed in Sect. 4. 3 of that book. 
When m H ^> m v (type II superconductors), vortices repel each other. When 
m H m v (type I superconductors) an attraction between vortices occurs. 

• Exercise: Consider two equal and opposite electric charges ±e em- 
bedded in a dual superconductor. Assume that they are separated by a 
large distance L, such that a flux tube is formed, the energy of which is 
proportional to its length L. We can then write the energy of the flux 
tube in the form S = ae 2 L. Study how the energy of the system varies 
as a function of the electric charge e, for a fixed value of L. Show that, 
on the average, the energy increases linearly (and not quadratically) 
with e, because several flux tubes can form. 
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3.3.7 Divergencies of the London limit 



The London limit is the extreme case m# 3> my of a type II superconductor. 
In this limit, the energy ()3.57|) is be minimized when the field s (x) maintains 
its ground state value s — 1 (S — v) for all values of x. The equation for 
b (x) reduces then to: 

did.,. 
-- — —(xb)+b = 

ax x ax 

which is equivalent to the following equation for the field C (z) : 

-^---^-(pB)+m 2 v B = (3.61) 
dp pap 

The solution, which vanishes far from the vortex is: 

b (x) = aK x (ar) or B (z) = °^K X (m v p) (3.62) 
However, for small values of x, we have: 

x(Ki (x)f = x- 1 + (-In- -ln2 + 7- - ] x + O (x 3 ) (3.63) 
\ x 2 J 

so that the integral of ^s 2 b 2 , in the string tension ()3.57|) . produces a logarith- 
mic divergence at small x. This is the origin of the divergence obtained in 
the analytic expression (j3.106|) of the string tension in the London limit. The 
electric field, far from the sources, is E = — V x B and it has the singular 
behavior E (f) ~ —e( z )^ In {rayp) when p — > 0. This singular behavior does 
not occur in the Landau-Ginzburg model and Fig. 13.21 shows that it is also 
not observed in lattice simulations. 



3.4 Comparison of the Landau-Ginzburg model 
with lattice data 

In 1998, Bali, Schlichter and Schilling compared the flux tube formed 
by two static sources in a SU (2) lattice calculation in the maximal abelian 
gauge (see Chapt0J) with the Abrikosov-Nielsen-Olesen vortex formed in the 
Landau-Ginzburg model of a dual superconductor. They measured both the 
electric field and magnetic current. The flux tube is depicted in Fig l3.ll 
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Figure 3.2: The profiles of the electric field E and of the magnetic current, 
marked kg on the figure, are plotted as a function of the the distance x from 
the center of the flux tube. The results are obtained from a pure gauge 
SU (2) lattice calculation in the maximal abelian gauge |33j . The curves are 
obtained from the Landau- Ginzburg model of a dual superconductor. 
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Figure 1331 shows that the relation j mag = -Vx£ between the electric field 
(|3.55|) and the magnetic current (|3.56|) is satisfied, and they checked that a 
magnetic current circulates around the flux tube. They also measured the 
profiles of the electric field and of the magnetic current and compared it to 
the corresponding expressions (j3.55J) and (j3.56J) obtained for the Abrikosov- 
Nielsen-Olesen vortex. Figure shows the fit obtained with the parameters: 

m v = gv = 1.23 GeV, m H = 2vVb = 1.04 GeV (3.64) 

The values of my and uih show that system is a type I superconductor, 
but close to the border separating type I and type II superconductors. A 
similar conclusion was reached by Matsubara, Ejiri and Suzuki in an earlier 
1994 paper [SH] for both color SU (2) and SU (3). The fit appears to be 
good enough to be significant. Note that the electric field behaves quite 
regularly when x — > 0, that is, close to the z-axis. This is in contradiction 
with the electric field calculated in the London limit. In a 1999 paper j33], 
Gubarev, Ilgenfritz, Polikarpov and Suzuki fitted the same lattice data with 
the parameters: 

g = 5. 827 ±0.004 m v — 1.31 ± 0.07 GeV m H — 1-36 ± 0.01 GeV 

(3.65) 

The data are taken with a lattice spacing fitted to the observed string tension 
y /cr su(2) = 440 MeV. The Landau-Ginzburg dual superconductor gives a 
string tension equal to: 

\fo = 400.1 ± 53.0 MeV ~ O.Ol^/asx^iy (3.66) 

A more recent fit reported in the 2003 paper by Koma, Ilgenfritz and Suzuki 
jnOl yields the following parameters: 

m v = 953(20) MeV m# = 1091(7) MeV (3.67) 

In this paper, the authors find a small dependence of the fitted coupling con- 
stant g on the distance separating the quark and antiquark, which is com- 
patible with antiscreening of the effective QCD coupling constant, derived 
from the Dirac quantization condition e e ff = 

Suganuma and Toki [UT] argue in favor of a different set of model param- 
eters, namely: 

m v = 0.5 GeV m H = 1.26 MeV (3.68) 
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Figure 3.3: The figure is a lattice confirmation of the "Ampere law" j mag = 
— V x E. The profiles of the magnetic current, marked k on the figure and of 
V x E, marked curlE on the figure, are plotted as a function of the distance 
x from the center of the flux tube. The results are obtained from a pure 
gauge SU (2) lattice calculation in the maximal abelian gauge 



53 



which would make the QCD ground state is a type II superconductor. They 
do not, however, fit the profiles of the Abrikosov-Nielsen-Olesen vortex. 

The scalar field S = which acts as an order parameter in the Landau- 
Ginzburg model, does not, as such, specify the nature of the monopole 
condensation, assumed to occur in the QCD ground state. Similarly, the 
Landau-Ginzburg model of usual superconductors was not direct evidence 
of electron-pair condensation, which was postulated and discovered six years 
later |H21 ESI 

3.5 The dielectric function of the color-dielectric 
model 

The Abrikosov-Nielsen-Olesen vortex may also be described in terms of the 
color dielectric model of T.D.Lee [HI]. The model is described by the action: 

h (A a) = J d'x (a) (d ^Af+ l - (da) 2 - \ {a 2 - v 2 f - j ■ A^j 

(3.69) 

where k (a) is a dielectric constant chosen such that is decreases regularly 
from 1 to as o varies from zero to v and where is a source for electric 
charges and currents. The equations of motion are: 

d ■ k (d A A) = j 

— d 2 a — - (<9 A A) 2 k (a) — 2b (a 2 — f 2 ) cr = ( k' (x) = (3-70) 
2 \ oa (x) J 

If we define the field strength tensor to be = k (a) (d A A) then 
the equation of motion for N 1 becomes equivalent to the Maxwell equation 
d ■ F = j. The system also develops a magnetic current 1 , given by: 

3 ma9 = d ■ F = (Ok) -WAA (3.71) 

Let us write: 

A» =(<!>,£) f=( P j) (3.72) 

1 I am indebted to Gunnar Martens for showing this to me prior to the publication of 
his calculations. 
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The action can be broken down to: 

d 4 x | -k (-d t A - V0) -- K (vx A) 

+\ id t af - \ (V*) 2 -\{a 2 - v 2 f - p0 + j ■ l} (3.73) 
The field imposes the constraint: 

V-k (-dtA-V^ =p (3.74) 
In the presence of static charges: 

p(x)=p(r) j( x )=0 

the source j M is time independent and the fields A^ and o may also be 
assumed to be time independent. The energy in the presence of static sources 
is: 

£ p (A,<f>,a) = J d 3 r (-± K (V0) 2 + \k (v x a)* + \ (w)' + \ (o 2 - v 2 \ 

(3.75) 

and the constraint is: 

- k (Vaj ■ (V^) - kV 2 = p (3.76) 
The energy is minimized when: 

VxA = (3.77) 
in which case the energy reduces to: 

£p(M = /^(-^(v0) 2 + ^(Va) 2 + ^(a 2 -t; 2 ) 2 -p0^) (3.78) 

The field o which minimizes the energy satisfies the equation: 

- VV + 26 (a 2 -v 2 )a-^K (v^) 2 = (3.79) 

Consider two static charges on the z-axis, placed symmetrically with re- 
spect to the origin: 

p(f) = - e s(f-R^ +e s(f-R 2 ^j (3.80) 
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If the charges are far apart, then, close to the x — y plane, we expect the 
electric field to be constant and directed along e^ z y A solution of equations 
(I3.76J1 and ()3.79|) exists in cylindrical coordinates, of the form: 

a (f) = a (p) <p (f) = az (3.81) 

and the function a (p) is the solution of the equation: 

- --r ( p-r) + 2b (° 2 -v 2 )v~ ^a 2 K = (3.82) 
pap \ dp J 2 

In this region, the electric and magnetic fields are: 

E{p) = -kV^= -e w o«(p) H = (3.83) 

In this sense, lattice calculations, as well as the Landau- Ginzburg model 
which agrees with them, determine the function k [a (/?)] of the dielectric 

model. The magnetic current j!^ ag = \p m ag-,jmagj is suc h that p mag = and: 

Jmag = - V x E = (v^j x = k' (v ' aj x (v^) = -e>)K^ (3.84) 

so that a magnetic current flows around the vortex which may be viewed as 
an Abrikosov-Nielsen-Olesen vortex. 

There is however an important difference between the color-dielectric 
model described by the action ()3.69|) and the Landau- Ginzburg model. If 
the system, described by the color-dielectric model, is perturbed or raised to 
a finite temperature, the expectation value of the field a may be shifted to, 
typically, a lower value and the color dielectric function no longer vanishes, in 
spite of having been fine-tuned to do so in the vacuum at zero temperature. 
The model therefore describes the physical vacuum in terms of an unstable 
phase, which is not the case of the Landau- Ginzburg model which describes 
the superconducting phase within a finite range of temperatures. 

3.6 The London limit of the Landau- Ginzburg 
model 

The London limit of the Landau Ginzburg model is obtained by letting b — ► 
oo, which means that the Higgs mass m# is much larger than the vector 
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meson mass my. In that limit, the scalar field maintains its vacuum value 
S = v and the model action (|3.8|) reduces to: 



= J d 4 x^(dAB + G) 2 + ^(B + dip) 2 ^j (3.85) 

where my, the mass of the vector boson, is equal to the inverse penetration 
depth (|3.12|) . and where the source term G is related to the electric current 
j by the equation: 

d-G = j (3.86) 

Because of its simplicity (the action is a quadratic form of the single field 
5 M ), the London limit has been extensively studied. It yields simple analytic 
forms for the confining potential, field strength correlators, and more. How- 
ever, as explained in Sect. 13.3.41 the linear confining potential develops an 
ultraviolet divergence because, in the London limit, the scalar field S is not 
allowed to vanish in the center of the flux tube. For the same reason, the 
electric field develops a singularity in the center of the flux tube and this is 
not observed in lattice simulations (see Sect. I3.4J) . 

We shall study the London limit for a system consisting of two static 
electric charges. The electric current is then: 

f = (p(f), 0,0,0) #f = -StfPp 

p (f) = eS (r — r*i) — e5 (r — r 2 ) (3.87) 

We work in the unitary gauge and we assume a straight line Dirac string, 
running from one charge to the other. In this case, we can use the form 

flnnnp : 

G = ^—nAj (3.88) 
n ■ a 

with a space-like vector = (0, ft) and with n parallel to the line joining 
the charges: 

n = f\ 2 — f\ 

n M = (0,n) n 2 = -n 2 n-j = (n-d)j = n-Vp (3.89) 
The equation of motion for the field is: 

d ■ (d A B + G) + m 2 v B = (3.90) 
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If we decompose B^ into longitudinal and transverse parts using the projec- 
tors (jA.lOj) . we find that B^ is transverse and equal to: 

B --wh? v d - v (3 - 91) 

We can then eliminate the field B^ from the action, which reduces to: 

I i =Jd' X (- 1 -( 9 .G) w l^-( 9 .G)- 1 -G^ (3.92) 

3.6.1 The gluon propagator 

When G has the form (J3.88)) . we have: 

(d ■ Gf = -L (d ■ ^XjY = 1 £ ^d a n ph G 2 = - (\ (n A j 

v ' n ■ o v ' n ■ o \n ■ o 

(3.93) 

Substituting, the action ()3.92|) becomes: 



L = I d 4 x ( ~j 3 (e^dmi) — 1 „ - 1 9 (e„ 456 5 4 n 5 j 6 ) — - (n A j) — * (n A j 

We can use (jA.8|) to evaluate e^ 123 e^ 5e . Remembering that = 0, a 

straightforward, albeit risky calculation yields the action: 

h = I ^ (ij , ( d2n2 ~ ( n ■ 9 f) o 2 ^ 2 -? - \ ( n ■ j) , 9 ~2a2 rZT2 ( n • j) 



2 J ( n -df K v J/ d 2 + m 2 v J 2 V JJ ( n • df d 2 + 

(3.94) 

- j d'x (-\ (n ■ j) —L, (n • j)) (3-95) 

We group together the terms which depend on n - d and those which do not, 
to get: 

f 4 . / 1 <T 1 n 2 ( m 2 v W n"n» 

(3.96) 

Since j M is a source term for the gluon field A^, the London limit of the gluon 
propagator in the dual superconductor can be read off the expression (|H.9(jj) : 

1 f 1 n 2 ( m 2 \ ( „„ n*V\ 
^ = _ o aTTTT " o—^2 ( *rr=T " " -2- (3-97) 



2<9 2 + m 2 / 2( n -9) 2 V«9 2 + 



my / \ nr 
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3.6.2 The energy in the presence of static electric charges 
in the London limit 

When the system consists of static electric charges, the fields are time- 
independent and the energy density is equal to minus the lagrangian. The 
energy obtained from the action l3~§o1 is thus: 



d r U o 



1 9 



1 n 2 



2d 2 + m 2 v 2 ( n ■ d) 2 \d 2 + my 
( 



d 3 r-p 



1 



n 



m 



v 



\ 



-V 2 + m 2 y 
( 



n ■ Vp 



(2tt 



If we substitute the form 



1 



-V 2 + m^ 



n 2 
\ 

I 



J, 



P 



n 



fi- 



le 2 + m 2 v (k 2 + m 2 v ) ^ _ £ 

of p into the energy (|3.98jl . we obtain: 



(3.98) 



/ 



£ P = - I d nd r 2 p (fi) v (fi - f 2 ) p (f 2 ) 



(3.99) 



with: 



v r 



(27TY 



d 3 ke- ik ^ 



\ 



m\ 



fc 2 + < (k 2 + m 2 v ) ( f ^ 2 



(3.100) 



where we set n = r in accordance with 



3.6.3 The confining potential in the London limit 

The first term of the potential (|3.100|) is a short ranged Yukawa potential: 



-p 2 r i -p 2 



vsr [r 



d ke 



-myr 



(27r) 3 J k 2 + my Arrr 



Consider the second (long-range) term: 

-e 2 

vlr (r) = - — — / d 3 ke~ 



(27T) 



ik-r 



m 2 v 



(k 2 + m 2 v ) / f .- 2 



(3.101) 



(3.102) 
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We write k ■ f = kr cos 9 and the long range potential becomes: 

Vrn (jf\ - ^ [ J3L -tJfcrcosfl 

LR [ ] ~ (2vr) 3 J k 2 (k 2 + m 2 v ) cos 2 9 

4 " r dk-^L_ C dx C0S{1 T x) (3.103) 



x 2 



(2vr) 3 7 (P + m^jQ 

The integral diverges both at small x and at large k. The divergence at small 
x contributes an infinite term which is independent of r. This can be seen 
by making a subtraction, that is, by evaluating v^r (r) — vlr (r ). Indeed, 
we have: 

f 1 cos (krx) — cos (kr x) , , , s i r rwr \ 

/ z ax = — cos fcr — krbi (kr) + cos kr + /cr M f Atq ) 

Jo x 

(3.104) 

so that: 

4vre 2 2 /- Ar , 1 



wl/? (r) = (f ) - ——^m v r / */^- 2^ (- cosy - ySi (y)) 

(27r) Jo (2/ + < r ) 

(3.105) 

Except for regions where myr < 1, we can approximate the function cosy + 
ySi (y) by its asymptotic value |y and, adding the short range part (|3.101|) . 
the potential becomes: 

vlr (rO = v LR (r ) - ^e~^ + ^£ r In (3. 106) 

47rr 87r m^ 

where a sharp cut-off A was introduced to make the /c-integral converge at 
large k. As discussed in Sects J3.3Tl and l3~ol the ultraviolet divergence is an 
artifact of the London limit. 

Equal and opposite electric charges are thus confined by a linearly rising 
potential. Let us define the string tension a by writing the potential in the 
form V (r) = — r v — h ar. The London limit of the Landau- Ginzburg 
model produces a string tension equal to: 

2™ 2 A2 1 ™2 ^2 



em 



v A 2 + m v _ n 2 7i 2 A 2 + 



a = — J! In ^ = — t, 2 In ^ (3.107) 

8n m v 2 m^ 

where we used the relation eg = 2mr between the magnetic and electric 
charges. The string tension may be compared to the prediction (|3.57J) ob- 
tained in the Landau- Ginzburg model, without taking the London limit, 
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namely: 



(3.108) 

The latter depends on two parameters, namely v and the ratio of the 
vector and Higgs masses. The cut-off in the expression obtained in the Lon- 
don limit mimics the missing vanishing of the scalar field in the vicinity of 
the vortex. 



3.6.4 Chiral symmetry breaking 

The gluon propagator (|H.97jl has been used by Suganuma and Toki as an 
input for a Schwinger- Dyson calculation, in which the quark propagator is 
dressed by a sum of rainbow diagrams in the Landau gauge [S3]. The quark 
gluon propagator was found to be strong enough to produce spontaneous 
chiral symmetry breaking. The dependence of the gluon propagator (|H.97j) 
on the direction of the vector n fl was averaged out. It was further modified at 
large euclidean momenta so as to make it merge with the perturbative QCD 
value, and at small momenta so as eliminate its divergence as k — > 0. As a 
result, the statement that the Landau-Ginzburg model explains the observed 
chiral symmetry breaking is only qualitative. However, they did observe that, 
when the vector mass my = gv was small enough, both confinement and 
chiral symmetry vanished. More recently, the relation between monopoles 
and instantons has been studied in Abelian projected QCD, both analytically 
and by lattice simulations [HHl EH EH EH| • The topic is reviewed by Toki and 
Suganuma [6T] . 



3.7 The field-strength correlator 

In these lectures, we have expressed 4-vectors and tensors in a Minkowski 
space with the metric g^ v = diag (1, — 1, — 1, — 1), such that detg = — 1. 
Such "Minkowski actions" are expressed in terms of vectors and tensors in 
Minkowski space, and they are suitable for canonical quantization and for var- 
ious representations of the evolution operator e lHt , by means of path integrals 
for example. However, present day lattice calculations are limited to evalua- 
tions of the partition function Tre~^ H , which is expressed in terms of a func- 
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tional integral of a Euclidean action. The latter is written in terms of vectors 
and tensors in a Euclidean space with a metric g^ u = diag (1, 1, 1, 1) = 5^ u , 
such that detg = +1. A Euclidean action is obtained when the functional 
integral for the partition function is derived from the Hamiltonian of the 
system. Crudely speaking, the Euclidean action which is thus obtained, is 
related to the Minkowski action by making the fi = components of vectors 
and tensors imaginary. This is summarized in the table (jB.lJ) of App[Bl The 
reader who is not familiar with the use of Euclidean actions, is referred to 
standard textbooks jTHl I7T| 172] . 

The Euclidean form of the Landau- Ginzburg action ()3.8|) in the unitary 
gauge, as given in App[BJ has the form: 

h (B, S) = Jd^x Q (dAB + G) 2 + 9 ^B> + \ (dSf + \b (S 2 - „ 2 ) 2 ) 

(3.109) 

In a 1998 paper [56], Baker, Brambilla, Dosch and Vairo suggested to use 
this Landau- Ginzburg action in order to model the field strength correlator 
which is observed in lattice calculations, in the absence of quarks. Formally, 
the source term G appearing in the action, can be used as a source term for 
the dual field tensor F = d A B. 

Let us illustrate the method by restricting ourselves to the London limit 
(Sect. 13. 6|) in which the field S maintains its ground state value S = v. The 
euclidean action reduces then to: 

Ij (B) = J d A x Q (d A B + G) 2 + ^B^j (3.110) 

The equation of motion of the field B^ is: 

-9- (dAB + G) +m 2 v B = (3.111) 

By considering the longitudinal and transverse parts, defined in App. IA.3.1|) . 
we can solve this equation in the form: 

B= * 2 d-G (3.112) 
—a 2 + my 

Substituting back into the action, we obtain: 

k = Jd'x (-1 (B-B) (a ■ a) + \g>) (3.113) 
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We may drop the second term which is simply due to the fact that the source 
term appears quadratically in the action. Use (|A.48|) to write the action in 
the form: 

h = f d^-G (-^-t) G=*[d*kG (*) G (-*) 



2 \-d 2 + m 2 v J 2(27r) 4 i w \k 2 + 



m 



v 



(3.114) 

where K^^p is the differential operator defined in (|A.42|) . In the Euclidean 
metric, it has the form: 

K^ a p = {Snadudp - Svad^dp + 5 u pd^d a - 5^d u d a ) (3.115) 

and its properties are listed in (jB.8|) . From the form (|3.11U|) of the action, 
we see that G may be considered as a source term for the dual field strength 
F = d A B . The Fourier transform of the dual field strength propagator can 
be read off (jSZH: 

/ d'xe-^ (T [F (x) F (0)] > = j^L. (3.116) 

Let us define the function: 

L (x) = -^— A [ d 4 k e ik - x - 1 2 = -^— A — i°° , 9 ^ 2 Ji (kx) dk 
V ' (2tt) 4 J k 2 + m 2 v (2vr) 4 x J k 2 + m 2 v V ; 

2 (^=Fl) (3-117) 



(27r) wyx 

The dual field strength propagator can be written as: 

(T[F(x)F(0)}) = ^jd 4 k w ^- = K^L(x) (3.118) 

Since F (x) = eF (x) and since, in the euclidean metric, we have eKe = E = 
G — K, the field strength propagator is: 

(T [F (x) F (0)]) — (G- K) d 2 L (x) (3.119) 

where G^ a/ 3 = 5^ a 5 u/ 3 - 5^5^. 

In order to conform to the conventional notation found in the literature, 
we note that L is a function of x 2 = x^x^ so that d^L = 2x^jj^. Substituting 
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the form (|3.115|) into the expression (|3.119|) . we find that the field strength 
propagator can be written in the form: 

(T [F^ (x) F aP (0)]) = (5^5„(3 - 5tf5 va ) d 2 L 

- 2 (5^ a d u Xf3 - Svad^Xp + Svpd^Xa - StfdvXa) (3.120) 

The contention of Baker, Brambilla, Dosch and Vairo 56 3 is that this field 
strength propagator can be compared to the gauge invariant field strength 
correlator (g 2 F tlu (x) U (x, 0) F a p (0) U (0, x)) which is measured on the lattice 
[73*1 [Jj] and which is usually parametrized in terms of two functions D (x 2 ) 
and Di (x 2 ): 

(g 2 F^ v (x) U {x, 0) F a(3 (0) U (0, x)) = (S^ a 5 u(3 - S^5 ua ) g 2 D (x 2 ) 

+ - {SftadyXp - S ua d^Xf3 + SvpdpXa - 5^d u x a ) g 2 D 1 (x 2 ) (3.121) 

By comparing the expressions (|3.12(J|) and (|3.121|) . the Landau-Ginzburg 
model makes the following predictions, in the London limit, for the func- 
tions D (x 2 ) and Di (x 2 ): 

g 2 D (x 2 ) =d 2 S = 5 (x) - m 2 v L (x) =*(*)- 
m v f 2Ki (m v x) 



2tt 2 x 2 V myx 



+ K (m v x)\ (3.122) 



g Di (x 2 ) = -4— = —f^ U V ' + K {m v x) 3.123 
ax z 2ir z x z \ myx J 

The cooled lattice data can be fitted, in the range 0.1 fm < x < 1 fm 
with the parametrization: 

D [x J = Ae T o H — t p L>i (x J = -oe T f H — ' .r = | »' 



A = 128 GeF 4 5 = 27 Gel/ 4 a = 0.69 b = 0.46 

T g = 0.22 /m T p = 0.42 fm (3.124) 

, 

The t p terms are negligible when x > 0.2 fm and, since the London limit 
is unreliable at small x, we neglect these terms for the comparison. Even so, 
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the fit to the shape is qualitative at best. The exponential decrease of the 
correlator allows us to make the identification my ~ ^- = 0.9 GeV, which is 
of the same order of magnitude as the values quoted in Sect J3.4I and obtained 
by fitting the profiles of the electric field and magnetic currents to lattice 
data. The reader is referred to the 1998 paper of Baker, Brambilla, Dosch 
and Vairo jSHI for the improvements obtained beyond the London limit. 

3.8 The London limit expressed in terms of 
a Kalb-Ramond field 

The confinement of static electric charges was modeled in Sect 13.61 in the 
London limit of the Landau- Ginzburg model. The same confining force is 
obtained in terms of the following model action: 

Ijj (A, $) = J d A x (~ (d ■ $) 2 - ^ (d A A - m$) 2 - j ■ A — mG ■ <s\ 

(3.125) 

which is expressed in terms of an antisymmetric tensor field $^ = — <& uti and 
its dual <£>, often referred to as a Kalb-Ramond field. In this action, m is a 
mass parameter, A^ is the gauge field which couples to the electric current 
j M and J^ v = —J UfJi is an antisymmetric source term for the field $. The 
latter is introduced so as to maintain the gauge invariance discussed below. 
The action (J3.125|) was studied in 1974 by Kalb and Ramond in the context 
of interactions between strings [73]. The duality transformation which leads 
to the use of Kalb-Ramond fields is proposed in the 1984 and 1994 papers of 
Orland [TH] [77] . Confining lagrangians are also expressed in terms of Kalb- 
Ramond fields in the 1996 paper of Hosek J7H]. In 2001, Ellwanger and 
Wschebor proposed a similar action to model low energy QCD and they 
showed that it confines [75] [HD]. A linear confining potential was also derived 
from it in the 2002 paper of Deguchi and Kokubo [81J. 

Let us show that this action is equivalent to the London limit of the 
Landau- Ginzburg action. The kinetic term of the action ()3.134j) is often 
written in the literature in terms of the antisymmetric tensor: 

Fafr = d a e^ + dpQ ia + 9 7 6 Q/3 F 2 = \F aM F a ^ (3.126) 
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It is simple to check that: 

" \ (d ■ e) 2 = i W"* 7 = \F 2 (3.127) 
so that the action ()3.134j) is often written in the form [7H] : 

^ (6) = J <fx Qf 2 - ^0 2 - mGej (3.128) 
3.8.1 The double gauge invar iance 

The action ()3.125j) is invariant under the usual abelian gauge transformation: 

A^A + (da) (3.129) 
It is however, also invariant under the joint gauge transformation: 

A^A + L $^$ + -(9AL) (3.130) 

m 



because d ■ d A L = 0. In the presence of the sources and J 1 *", the double 
gauge invariance is maintained, provided that the sources satisfy the (com- 
patible) equations: 

3-j = d-G = j (3.131) 

The second gauge invariance relates the source term J to the electric current 
j. We can choose L = —A so as to write the action in a form in which the 
double gauge invariance is explicit: 

Ij (A, $) = J d 4 x (d- $) 2 - ^ (d A A — m$) 2 + G{d A A — m$)^ 

(3.132) 

We can define an antisymmetric field Q^ u : 

dAA-m$ = -mQ 6 = $ — —d A A (3.133) 

m 

The field Q^ v has the property of being invariant under the gauge transfor- 
mation (j3.13(Jj) . Noting that d ■ $ = d ■ 6, the action (J3.132|) can be written 
functional of G alone, namely: 

i 5 (e) = y d 4 x (~ (o ■ e) 2 - ^e 2 - mcej (3.134) 
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This is equivalent to a choice of gauge in which the gauge field is absorbed 
by the field 9^. 

We can use (jA.49|) to write the action in the form: 

I. (6) = J d 4 x Qe (Ed 2 - m 2 ) 9 - mGej = J d A x Qe (E (d 2 + m 2 ) - m 2 K) 9 - mGQ^j 

(3.135) 

where K and E are the longitudinal and transverse projectors (|A.42|) . The 
action is stationary with respect to variations of 9 if: 

(E(d 2 + m 2 ) -m 2 K)e = mG 9 = f — - - - ^ ] mG (3.136) 
v \ / / \d 2 + m 2 m 2 J 

We can eliminate 9 from the action, to get: 

I 3 = - fd*x l - G ( f E -k) G 
3 J 2 \d 2 + m 2 J 

=/ A H (8 - S) (^) (9 ' 5) -H (3i37) 

This form is identical to the expression (|3.92|) obtained in the London 
limit of the Landau- Ginzburg model, with m = my. The source term G 
satisfies the equation d • G = j in both cases, and a straight line string 
G = j^^Qj (n A j) leads to the same confining potential ()3.106|) . This is not a 
coincidence, because we shall next display a so-called duality transformation, 
by means of which the Landau- Ginzburg model can be expressed in terms of 
a Kalb-Ramond field (not only in the London limit). 

3.8.2 The duality transformation 

Let us show that the Landau- Ginzburg model can be expressed in terms of 
an antisymmetric tensor field 9. For this purpose, we add to the Landau- 
Ginzburg action (|3.8|) the term: 

where m is a constant mass. Adding such a term is permissible, because the 
equation of motion of the field 9 simply makes the added term vanish. The 
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resulting action is: 

,2\2 



J, (0, B,S,<p) = j d 4 x{~{dAB + G) 2 + ^(B + d V f + ~ («9S) 2 - ift (S 2 - ^ 2 ) : 

+^ (® + ^( 9AS + (5 )) 2 } ( 3 - 139 ) 

The added term is chosen such that the term — \ [d A B + G) 2 cancels. We 
are left with the action: 

Ij (6, B, S,<p) = J d 4 x 

I !^!e 2 + mO {d A B) + mG0 + 9 — (B + dip) 2 + - (OS) 2 - -b (S 2 - v 2 ) 2 ) 
I 2 2 2 2 1 

(3.140) 

The identities (|A.28j) and (|A.46|) allow us to write: 

-m J d 4 xC (d ■ 6) = m J d 4 x (B + dip) (d ■ 6) (3.141) 

so that the action can be cast into the form: 

I j (e,B,S,<p) = jd'xi^f- ^B + dip-m^^- 

~Wsi ( d ' ^ + + m0 ^ + \ (ds)2 ~ \ b ^ - (3 ' 142) 

The first term of the action can be dropped because it vanishes when the 
equation of motion for the field B is satisfied. The remaining action is: 

J, (6, S) = J d*x (-J^ (d-Q) 2 + ^Q + mGQ 2 + \ (dS) 2 - \b (S 2 - v 2 ) 2 

(3.143) 

This action, expressed in terms of the Kalb-Ramond field G, is identical to 
the action ()3.8j) of the Landau-Ginzburg model. In the London limit, where 
S — v, we can choose m = v and the action becomes identical to the Kalb- 
Ramond action (|3.134j) . 

The transformation of the action Ij (0, S), expressed in terms of a Kalb- 
Ramond field Q^, into an action Ij (6, B, S, ip) which depends on the "dual" 
vector field is called a duality transformation. It is widely used in the 
literature [HHH3 , [HI EH • |Z3ISQ|- 
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3.8.3 The quantification of the massive Kalb-Ramond 
field 

To visualize the physical content of the action (J3.134j) . let us express the Kalb 
Ramond field in terms of two cartesian 3- dimensional fields e and h: 

e i = _ Q oi = ~ £ W k e jk = -9 0i = -^e 0ilk Q jk (3.144) 



It is an easy exercise to check that: 

- I (d ■ e) 2 = -\ (V ■ hf + \ (d t h + V x ef (3.145) 



The source term G^ v can, in turn, be written in terms of the two cartesian 
3-dimensional vectors E st and H st defined in (J2.9U|) : 

E\ t = -G 0i = \e 0i3k G jk Hi t = -\^ jk Gjk = ~G 0i (3.146) 

G- $ = -e-E st + h-H st V-E st = p - d t E st - V x H st = j (3.147) 

where j M = (^p, j\ is the electric current. 

The action ()3.134|1 can thus be broken down to: 



m _n i / -> -a L m 



2 



I p3{^ h ) =J d'x^-^h + Vxe) +—?--[V-h) -—h> 

+me ■ E st - mh ■ H st j (3. 148) 
No time derivative acts on the field e so that it acts as the constraint: 

V x (d t h + V x ej + m 2 e = -mE st (3.149) 

The conjugate of the field h is: 

jf = d t h + Vxe (3.150) 
and the conjugate of the field 5 is: 

P=(d t S) (3.151) 



69 



Taking the constraint (|3.149j) into account, the hamiltonian, or classical en- 
ergy, is: 

n(n,h\= f d 3 r (n ■ (d t h) + P (d t S) - L 



d 3 r 



r vf 2 



— + ^ (V x 7? + mE st y + 1 (v • Ky + ^m 2 h 2 + mh ■ H st 



(3.152) 

In the absence of sources, the classical energy is a sum of positive terms. 



3.8.4 The elementary excitations 

Let us expand the hamiltonian (|3.152|) to second order in the fields around 
their vacuum values S = v , tt = h = 0, in the absence of sources. The 
Hamiltonian reduces to: 



H[T?,h 



d 3 r 
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2 1 



, .Vxjfl +-[V-h) + -m 2 h 2 
2 2m 2 V / 2 V / 2 



2 1 



d r 



^tt t (-V 2 + m 2 ) 7? T + X -m 2 h 2 T + y + ^l (-V 2 + m 2 ) h L 



(3.153) 

where we used flA.84|) . 

As in Sect J3 .1(11 we expand the field h and its conjugate tt on the plane 
wave basis (|3.169j) : 



(ri 



ka 

;r ( /■ ) -= — = — (In 



ka 



k(1 ) + 4a l rl /) = ^ W + ^ ^ 



fca ^ a %a - al a (ka \ri ^ = frf (f) + fff (f) 



(3.154) 

where the transverse and longitudinal parts correspond respectively to the 
a = 1, 2 and a = 3 contributions. If we assume boson commutation rules 



kai k'b 



SkpSab, the fields h and n become quantized: 
[ir* (r) , h j (r*)] = 



(3.155) 
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Upon substitution of the expansion ([3.154]) into the hamiltonian ([3.153]) . we 
find that the following coefficients a^ a cancel the [c^^a ^ + at a} ^ j terms 
in the hamiltonian: 

1 

«fc,a=l,2 — T 

(k 2 + m 2 ^ 4 

With these coefficients, the hamiltonian ([3.153]) acquires the diagonal form: 
H = \l& + m 2 [a\ lbHa + (3.157) 

ka 

so that the elementary excitations of the massive Kalb-Ramond hamiltonian 
consists of three vector particles of mass m, which is exactly the same as 
the elementary excitations of the Landau- Ginzburg model in the London 
limit. However, in the limit m — ■> of vanishing mass, the transverse field 
disappears and only one particle of mass m subsists. 

• Exercise: Let be the vector which is the dual form of the tensor 

F" = ^e^F aM F a ^ = e a ^F^ F 2 = —F 2 (3.158) 

Note that the duality transformation F — > F of tensors with an odd 
number of indices is reversible without a change in sign (App. IA.5j) . 
Show that the components of the vector F M are: 

p° = -V • h F i — d t ti + (y x eV (3.159) 

Show that, when m — 0, the lagrangian ()3.148|) can be expressed in 
terms of the four fields F^. The quantification of the fields in this limit 
is discussed in the 1974 paper of Kalb and Ramond [75] . 

3.8.5 The Nambu hierarchy of gauge potentials 

The kinetic term of the Kalb-Ramond action ()3.128j) is written in terms of 
an antisymmetric field tensor F a ^ and its associated gauge potential O^. 



a k ,a=3 = — I k +m 
m 



i 

2 \ 4 



(3.156) 
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Nambu has displayed a hierarchy of field tensors F^ u> F a/37 , ... which can be 
expressed in terms of gauge potentials A^, A^ u , ... [54] : 

Fp,v = d^Ay - d u A^ F^x = d^A uX F^ X p = d^A uXp (3.160) 

cycl cycl 

The respective lagrangians: 



L = --F,^ + gA,f 



12m 2 
1 



4877V 



F^xF^+gA^f 
F,„x P F^ Xp + gA, uX j^ x 



(3.161) 



are invariant under the gauge transformations: 



An -> A^+id^A) A„ v -> A^+{d^K - 9„A M ) -> A^a+J^ (^A^ a ) 

(3.162) 

provided that the currents are conserved: 

8^ = 8^ = 8^ = (3.163) 



3.9 The hamiltonian of the Landau-Ginzburg 
model 

Consider the action ()3.15|) . In the unitary gauge, ip = S is a real field and 
the action reduces to: 



d x 



I ^ dt B - \7 X + 2 - I (-V x B + E st N 2 



9 2 S 2 



5^)' -5* 



The constraint imposed by \ is: 

g 2 S 2 X + V ■ l-d t B-V X + H st ) 



(3.164) 



(3.165) 
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The conjugate momentum to the field B: 

51 



SldtB 



d t B + V X -H 



St 



-H 



(3.166) 



The minus sign on the right hand side is inserted in order to conform with 
the definition (gHH}) of the magnetic field. The momentum conjugate of the 
field S is: 

P=jjtf = m (3.167) 



Taking the constraint (j3.165|) into account, the hamiltonian, or classical en- 
ergy, is: 

H (H, B, P,S)= [ d 3 r (-H (d t B) + P (d t S)) - I 



1 1 / \ 2 n 2 Q 2 1 



2 2 V J 2 
In the absence of sources, the energy is a sum of positive terms 



V • // ) — H ■ H st 

(3.168) 



3.10 The elementary excitations of the Landau- 
Ginzburg model 

Vector fields can be expanded in the plane wave basis: 



ri 



V 



e ~ e \ a ) [ k ) (3-169) 



where the vectors et a =i 2 3) are orthogonal unit vectors with etg\ parallel to k: 



6(3) = ^ e( a ) x e( 6 ) = e a b c e( c ) 



(^123 — —^213 — 1) 



(3.170) 
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and where V is the volume in which the plane waves are normalized. The 
basis (|3.169j) is complete and orthogonal: 



ka i 



The fields B and H can be expanded in this basis: 



k'b ) = S^,S ab 
(3.171) 



B l (r) 



jfE^( 

ka 



ri 



ka)Ha + 4a l f V ) = ^ W + ^ ^ 



A;a ^ ag a — at ^/ca \ fi 



HI (f) + fl* (r) 



(3.172) 

The longitudinal and transverse parts Bl and i?r of B correspond respec- 
tively to the contributions of a = 3 and a = 1, 2. We have V x = and 
V • B T = 0. Similarly for the longitudinal and transverse parts of H. 
If we impose boson commutation rules on the a^ a and a\ a coefficients: 



— * —» 

the fields B and H become quantized: 

[W (f) , B j {r)] = i (ri \r j ) = i5 ij 5 (f - r) 
The scalar field may be similarly quantized: 

s (*o = -4 E«* (( r 1*) 6 s + 4 p)) p (*o = t4 E 



(3.173) 



(3.174) 



6t /fc 



[P(f),5(^)] = (r|f) =<J(f-0 (3.175) 

The elementary excitations of the vacuum, in the Landau- Ginzburg model, 
may be obtained by expanding the fields to second order in the vicinity of 
their vacuum values S = v,B = H = and in setting the sources E st and 
H st to zero. To second order in the fluctuating parts, the hamiltonian can 
be written in the form: 

H (il, C, P, s) = J d 3 r +^H* + ^B T (-V 2 + m 2 v ^j B T + (-V 2 + m\) H 
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(3.176) 



where we used (IA.84J) and where m# and my are the Higgs and vector 
masses ()3.10|) and (j3.11|) . If we substitute the expansions ()3.172|) into the 
hamiltonian ()3.176|) . we can choose the coefficients a^ a such that the terms 

[ar a t + at a - ) vanish, and the hamiltonian reduces to: 

V ka ~ ka ka -ka) ' 

ka k 

The elementary excitations of the Landau- Ginzburg model consists of three 
vector particles with mass my and a scalar particle with mass m#. 

The appearance of massive vector particles with masses barely higher 
than 1 GeV (see the values of the estimated vector and Higgs masses in Sect. 
13.4)1 has been invoked to criticize the dual superconductor model. Indeed, 
as we shall see in Sect. 15.11 the Higgs field is not a color singlet so that 
the model predicts the existence of freely propagating particles with non- 
vanishing color. 

• Exercise: show that the coefficients a ka are: 

atk,a=i,2 = (k 2 + m v ) 4 a kya=3 = ^- (k 2 + m v ^ 4 (3.178) 

Show that, in the limit my 0, the field Hl vanishes and that V--B = 
so that only two massless vector particles remain. 



3.11 The two-potential Zwanziger formalism 

In sections 12.21 and 12.111 we saw that, in the presence of both electric and 
magnetic currents, electrodynamics could be expressed either in terms of a 
gauge potential N 1 associated to the field strength F = d A A — G or a 
gauge potential associated to the dual field strength F = d A B + G. 
The corresponding actions ()2.13|) and (J2.88J1 involve non-local string terms. 
In 1971, Zwanziger wrote a beautiful paper in which he proposed a local 
lagrangian which uses two potentials A^ 1 and B^ jlH]- We shall apply his 
formalism to the Landau-Ginzburg action of a dual superconductor. The 



75 



theory was further developed, in particular its Lorentz invariance, by several 
authors. For a review see the 1979 paper of Brandt, Neri and Zwanziger [Ho] . 
The Lorentz invariance is also discussed in the more recent 1998 paper of 
Gubarev, Polikarpov and Zakharov Other aspects of the formalism are 
discussed in the 1976 and 1979 papers of Blagojevic and Senjanovic [8*71 IH""| . 
We shall show that, when applied to a dual superconductor, the Zwanziger 
action leads to the form (|3.8|) of the Landau- Ginzburg action. 

3.11.1 The field tensor F lw expressed in terms of two 
potentials and B lL 

Zwanziger writes the components n ■ F and n ■ F of the field strength tensor, 
along a given fixed 4- vector n M , in terms of two potentials, namely: 

n-F = n-(dAA) n ■ F = n ■ (d A B) (3.179) 

He then expresses the field strength tensor F in terms of these components, 
using the identity (|A.38|) : 

F = (n A (n ■ F) - n A (n - 

= -J (n A (n • (d A A)) - n A (n • (d A £))) (3.180) 
and the dual tensor is: 

F = ^ (n A (n ■ (d A A)) + n A (n • (<9 A £))) (3.181) 

The Maxwell equations d ■ F = j and d ■ F = j mag can be satisfied with 
potentials A^ and which satisfy the equations: 

^ ((n ■ df A - (n ■ d) d {n ■ A) - n (n ■ d) (d ■ A) + nd 2 (n ■ A) + (n ■ d) (n ■ dT^B)) 

\ ((n ■ df B - (n • d) d (n ■ B) - n (n ■ d) (d ■ B) + nd 2 (n ■ B) - (n ■ d) (n ■ JXA)) 

(3.182) 

where we used successively (IA.16J) . ()A.15|) and ()A.19|) . 
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• Exercise: Choose to be space-like and define: 

A>* B tx =(x,B^j n M = (0,n) n-n = l (3.183) 

Show that (I3.18U)) and (13.181)) break down to: 

- ( n . F)° = n ■ E = - [n ■ (d A A)] = n ■ (- V0 - d t Xj 

(n • F) 1 = - (n x My = [n ■ (d A A)f = - (ft x Jvxl)) 

(n • F)° = -n ■ H = [n ■ (d A B)}° = n ■ - d t B^j 

{n ■ Tf = (n x jf) = [n ■ (d A B)f -|nx|vxfl)) (3.184) 

so that the potential A describes the longitudinal part of the electric 
field and the transverse part of the magnetic field, whereas the po- 
tential B describes the longitudinal part of the magnetic field and the 
transverse part of the electric field. In this instance, the longitudinal 
and transverse parts of the fields are defined relative to the vector n 
and not relative to V as in (|A.84|) . 



3.11.2 The Zwanziger action applied to a dual super- 
conductor 

The Zwanziger action applied to a dual superconductor is: 

1 



I J (A,B,S, V ) = J tf 



x 



2n 2 



(n • (9 A A)) - j ■ A 



2n 2 



(n ■ d A A) ■ {n ■ d A B) + — (n ■ d A B) ■ {n ■ d A A) 



' (n ■ (8 A B)) 2 + ^ (B + d<p) 2 



2n 2 



+ 1 -(dS) 2 -h(S 2 -v 2 ) 2 



(3.185) 



The first line describes the gauge field A^ and its coupling to the electric 
current j M . The second line is an interaction between the two gauge potentials 
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and B^. The third line describes the interaction of the gauge field B^ 
with the complex scalar field Se igip of the Landau- Ginzburg model. The last 
term describes the dynamics of the order parameter S. 

The action (j3.185|) is invariant under the following gauge transformation 
of the field A^\ 

A^A + (da) (3.186) 
It is also invariant under the joint gauge transformation: 

B^B + (d(3) (p^np-p (3.187) 

The electric current j M might be caused by quarks, in which case the term 
—j ■ A would be replaced by the term: 

q [i (<9 M + ieA) 7m - m) q (3.188) 

When this is done, the action becomes a local action for a system of confined 
Dirac particles with electric charges. 
The following are useful identities: 

J d 4 x (n-dAA)-(n- d~AB) = - J d 4 x (n ■ d A B) ■ (n • d~AA) 

= / d 4 x e^A ^n v (n ■ d) d a B/3 

= - J d 4 xA(n-d)[n- dT^B] = + J d A x B (n ■ d) [n ■ d7^A\ (3.189) 
and: 

r d 4 x [n ■ {d A A)] [n ■ {d A A)} 



d A x [A (n ■ df A - A(n ■ d) d (n ■ A) - An (n ■ d) (d ■ A) + (A ■ n) d 2 (n ■ A)] 

(3.190) 



3.11.3 Elimination of the gauge potential 

Let us show how to eliminate the field A^ in order to reduce the action 
(|3.185J) to the form ()3.8|) of the Landau- Ginzburg action. Because the theory 
is invariant under the gauge transformation A — > A+ (da), we can, following 
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Zwanziger, add a gauge fixing term ^2 [d (n ■ A)] 2 to the action. After doing 
this, and in view of the identities ()3.189j) and the action ()3.185j) can 

be written in the form: 



Ij (A, B, S» = J d 4 x 



-ApM^K + A^ 



1 



n- 



(n-d) (n-dABf-^ 



9 2 S 2 



{B + d V f + l -{dS?- l -b{S 2 -v 2 ) 2 



where M^ v is the matrix: 



^ (n ■ (d A B)) 2 



(3.191) 



AP" = —( n -d) [{n ■ d) - d"n u - n^d" 



(3.192) 



The inverse matrix is: 



n 



(n ■ d) 

The field A 11 satisfies the equation of motion: 



(9» v ~ ~2 n ^ ~ Jp^^ (3.193) 



M"M„ = - 



— (n-d) (n-dABY -f 



(3.194) 



We can use the inverse matrix ()3.193|) to eliminate the field A^ from the 
action, which becomes: 



Ij (B,S,<p) = J (fx 



\(dAB) 2 



1 l n 2 

(d A B) aN n A j - -f 



(n-d) J 2 J (n-dY 



n 



9 2 S\ n , , 2 , l, Qm2 1 



} {B + d V y + -{dsy--b(s 2 -v 2 y 



where we used the property: 



d v (n ■ d A B) =0 n u (n ■ d A B) =0 



(3.195) 



(3.196) 



as well as the second Zwanziger identity (jA.39|) for the d A B terms. The 
action can finally be reduced to the form: 



0) 



n A j 
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+ ^ {B + d V f + \{dSf-\b{S 2 -v 2 ) 2 



(3.197) 



This is precisely the form f!3.8|) of the Landau- Ginzburg action, with a 
straight-line string term G = j^j^n A j. 

There is however an apparent difference. The Landau- Ginzburg action 
(J3.8j) contains a non-local string term G which we may choose to be a string 
which stems from a positive charge and terminates on a negative charge. In 
the local Zwanziger action ()3.185|) . the vector n M is fixed and independent of 
the position of the charges. This is a source of difficulties, because, when the 
Zwanziger action is used with classical fields, it breaks Lorentz invariance, as 
discussed at the beginning of Sect 13. Ill 



Exercise: Consider projectors K^^p and E^^p defined, not in terms 
of as in (|A.42|) . but in terms of the given vector n M : 

K^ap = {9^n u np - g ya n^np + g u pn^n a - g^n y n a ) 



rr 

(3.198) 

Check that the projectors K and E satisfy the relations: 

K 2 = K E 2 = —E KE = K—E — G e 2 = -G (3.199) 
Check that: 

1 1 ; — =r 

KF = —n A (n ■ F) EF = —n A (n-F) (3.200) 

n z n 

Show that the Zwanziger identity (|A.38J) can be expressed in the form 
F = (K — E)F. Show that Zwanziger expresses the "longitudinal" 
part KF of the field strength F^ v in terms of a vector potential 
and the "transverse" part EF in terms of the potential B^: 



KF = K (d A A) EF = —Ed A B (3.201) 

Exercise: Show that, if S^ v and T^ v are antisymmetric tensors , we 
have: 

SKT = — (n-S)(n-T) SET = SKT = — (n-S)(n-f) 
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SKT = —in- S)(n-T) = -SET SET = -SKT = — - (n ■ S) (n 

(3.202) 

Use d ■ d A A = to check that: 

d ■ (KJXA) = 8 ■ (eFXA) (3.203) 

Show that the Zwanziger action 1)3.185)1 can be written in the form: 

I jijmag (A,B) = jd'x j-i(«9 AA)K (d AA) -\{d AB) K (d A B) 



-- (d A A) Kd A B + - (d A B) Kd A A - j • A - j mag ■ B \ (3.204) 
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Chapter 4 

Abelian gauge fixing 



The formation of monopoles and their condensation in the QCD ground state 
is a feature which is related to abelian gauge fixing, discussed in this chapter. 
The gluon field acquires a singularity in the vicinity of points in space where 
abelian gauge fixing fails and magnetic monopoles are formed there. The 
ideas discussed in this chapter can be found in the 1974 and 1981 seminal 
papers of Polyakov jHHj and 't Hooft [OHl EE]- It is also very instructive 
to read the Sect. 23. 3 (vol.2) of Weinberg's Quantum Theory of Fields |52~] . 
The formation of monopoles in QCD is still a subject of occasional debate 
The choice of the abelian gauge is, of course, not unique, and nor 
is the corresponding definition of the monopoles. The recent gauge invariant 
definition of monopoles, proposed by Gubarev and Zakharov jOUESI, i s n °t 
discussed in these lectures. 

The dynamical formation of monopoles in the QCD ground state is not 
explained by the 't Hooft construction. A remarkable feature has however 
been confirmed by lattice calculations P3 1201 1211 123 , [21 123 121 , namel y the 
previously surmised condensation of monopoles in the QCD ground state. 
The lattice calculations evaluate the vacuum expectation value of an opera- 
tor which creates a magnetic monopole in an Abelian gauge. It is found that 
the vacuum expectation value of this operator is non-zero in the confining 
phase and zero in the deconfined phase, thereby signaling the condensation 
of magnetic monopoles in the confining phase. This is one of the main moti- 
vations for a phenomenological description of the QCD ground state in terms 
of dual superconductors. The lattice calculations suggest that the conden- 
sation of monopoles in the QCD ground state is remarkably independent of 
the gauge fixing condition. 
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The gluon field is a vector in color space. Of course, it would be nice 
if we could make a gauge transformation at every space-time point x which 
would rotate the gluon field so that only its diagonal components A% and A$ 
would remain. This would reduce QCD to an abelian theory. The trouble is, 
of course, that the gluon field has four components A^ and that it is only 
possible to align one component at a time. This is why a scalar field is most 
often used to fix a gauge. 

Let $ (x) be a scalar field in the adjoint representation of SU (N c ), which 
means that the field is a vector in color space with N% — 1 components <3> a (x). 
The field can be written in the form: 

<S>(x) = <5> a (x)T a (4.1) 

where T a are the N% — 1 generators of the SU (N c ) group. They are equal 
to one half of the Pauli matrices in the case of SU (2) and to one half of the 
Gell-Mann matrices in the case of SU (3). The field $ (x) does not have to 
be one of the fields appearing in the model lagrangian. The choice of <3> is 
not innocent and will be discussed below. The orientation of the vector $ (x) 
in color space, at the space-time point x, defines a gauge. Different choices 
of $ lead to different choices of the gauge. 

Consider how this is done in practice. Local rotations in color space are 
generated by operators of the form: 

n (x) = e iXa{x)Ta (4.2) 

The operators Q (x) are elements of the color SU (N c ) group. A local rota- 
tion in color space is called a gauge transformation. The generators T a are 
traceless hermitian N c x N c matrices, so that the field <3> (x) = $ a (x) T a may 
be viewed as a traceless matrix in color space. We can always perform a 
gauge transformation (a rotation of the vector <3>) so as to diagonalize the 
matrix <&(#). This means, that there always exists a rotation Q(x), such 
that: 

Q, (x) $ (x) fi f (x) = diag (Ai (x) , A 2 (x) , \ Nc (x)) (4.3) 

The gauge in which $ (x) is diagonal is called an abelian gauge. The abelian 
gauge depends, of course, on the choice of the scalar field $. 

When N c = 2, the abelian gauge is obtained by aligning the vector $ a (x) 
along the color 3-axis: 

$ (x) - n (x) $ (*) ^ ( X ) = $' 3 ( X ) T 3 = U P °_ $/ \ (4.4) 
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When N c = 3, two generators of the SU (3) group are diagonal, namely: 





(4.5) 



and the abelian gauge is obtained by aligning the vector <3> a (x) along T 3 and 

T 8 axes: 

$ (x) -> n (x) $ (x) ft f (x) = $3 (x) T 3 + $ 8 (x) T 8 



( $ 3 + 75*8 



V 
















2_d)' 

8 



\ 
/ 



(4.6) 



Further gauge transformations, generated by two diagonal generators T 3 and 
T 8 , leave the diagonal form of <3> (x) invariant. Such rotations have the form 
e ix3T 3 +ix 8 Ts = e ix3T 3e i X 8T 8 and they belong to the residual U (1) x U (1) sub- 
group of SU (3), called the maximal torus subgroup of SU (3). 



4.1 The occurrence of monopoles in an abelian 
gauge 

There are points in space where the abelian gauge fixing becomes ill defined. 
We shall see that such points, which are sometimes referred to as topological 
defects, are sources of magnetic monopoles. 

4.1.1 The magnetic charge of a SU (2) monopole 

Consider the SU (2) case first. Let Q (x) be the gauge transformation which 
brings the field $ (x) = $ a (x) T a into diagonal form: 

$ = $ a T a - nm^ = AT 3 = ^ q - A ) A = v $ i + $ 2 + $ i ( 4 -7) 

The eigenvalues A (x) of the matrix $ (x) are, of course, gauge indepen- 
dent 1 . A degeneracy of the eigenvalues of $ (x) occurs when A = 0. At any 

1 They do, however, depend on the choice of the field <I> (x). 
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one time, this implies that all three components $ a =i,2,3 {r) should vanish, 
and this can only occur at specific points f = r in space such that: 

$i(r ) = $ 2 (r ) = $ 3 (r ) = (4.8) 

The three equations determine the three components (x ,y ,z ) of the 
vector r . At the point r , defined by the equations ([4.8)1 . it is not possible 
to define the gauge and we shall that the gluon field develops a singularity 
at that point. 

In the vicinity of the point r , we can express $ (f) in terms of a Taylor 
expansion: 



<9$ 

$(f) = $ a (f)T a = T a C , ab (^-x 06 ) C a 6 



(4.9) 



r=r 



The matrix C a b defines a coordinate system in which the field $ (f ) has the 
form: 

$(?)=a? a T a x' a = C ab (x b -x ob ) (4.10) 

In this coordinate system, the solution r of equation (|4.8|) is placed at the 
origin, and the field $ (r) has the hedgehog shape displayed in Eq. (|4.10|) . In 
the following, we work in this coordinate frame and drop the primes on x'. 

Let (r, 9, <p) be the spherical coordinates of the vector r (see App lA.6."3j) . 
In spherical coordinates, the hedgehog field $ (r) = x a T a is represented by 
the matrix: 

$ if) = x a T a = T\T sin 9 cos (p + T 2 r sin 6 1 sin (p + T 3 r cos 9 

r ( cos 9 e~ iip sin 9 \ , , 

~2Ve^sin^ -cos# J ^ J 

The matrix which diagonalizes $ is: 

n {o, v) = ( e " cos l e) (fl, *0 = ( e ~" T 1 " sin % 

y \ -sin I e-^cosf / v ^ ; \ smf e^cosf 



(4.12) 



Indeed, we can check directly that: 



n<&nt = ^( I \ ) =tt 3 (4.i3) 
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Now consider how the gluon field transforms under the same gauge trans- 
formation 2 : 



a„ = -> n (a^ + A = A a T a n (l+ ^vj fi f 

(4.14 ) 

The expression of the gradient V in spherical coordinates is given by (|A.112j) . 
The vector fiVfi^ is: 

■* = s - ("I"*) + (">') + ^ ("I"') (415) 

where the unit vectors e are defined in (|A.110|) . From the explicit expression 
(11121 of tt, we find: 

ft-^fit = o 

j9_ Q t _ £ / -cos0 - 1 e^sinfl 
<V ~ 2 V e-^sin^ cos#+ 1 

= — % (1 + cos #) T 3 + z sin 6* cos y?Ti — i sin 6 1 sin <^T 2 (4-16) 

so that: 

— VtWVt ] = - \-e e T 2 e %ip - e j~ + C °^ T 3 + e v - (cos cpTx - sin^T 2 ) ] 
ie e \ rsmu r J 

(4.17) 

The terms are all regular except for the term: 



I(nvnt) =-V±i^r s (4.i8) 

which becomes singular when 9 — > 0, that is, on the positive z-axis 



2 The gauge transformation 13.9|l . defined in section mi corresponds to a rotation 
f2 = e* 9/3 , in which the angle (3 is multiplied by the magnetic charge g. In the gauge 
transformation i|4.2[l . the angles Xa are not multiplied by the coupling constant e. This is 
why a factor - appears in the gauge transformation (|4.14() whereas no such factor appears 
in EHJ). 
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Thus, in the abelian gauge, obtained by diagonalizing the field $ (x), the 
gluon field can be separated into a regular part A R and the singular part 
flUED : 

A = A a T a = A*T a - UJ-^-T, (4.19) 

e r sin d 

Note that only the diagonal (abelian) 'part of the gluon field acquires a sin- 
gular form. The singular part (j4.18|) has exactly the form ()2.67j) which a 
gauge field acquires in the vicinity of a magnetic monopole situated at the 
origin, with a Dirac string running along the positive z-axis. By comparing 
the expressions (j4.19j) and (j2.fi 7)) . we see that the magnetic charge of the 
monopole is equal to: 

47T 

9 = T 3 (4.20) 

e 

Here e is the color electric charge, that is, the QCD coupling constant. This 
is another instance of the Dirac quantization condition ()2.80j) . In a way, the 
result (|4.2U|) is promising for low energy phenomena because perturbation 
theory points to a divergence of the QCD coupling constant e at low energy 
and we may expect g to be better behaved. The way the running coupling 
constant e and g manage to maintain a constant product eg is discussed in 
2001 papers of the Russian group [HElinij- In short, we have shown that, in 
the vicinity of points where the eigenvalues of the matrix $ (x) are degener- 
ate, that is, at points where the abelian gauge is ill defined, the abelian part 
of the gluon field behaves as if a monopole with magnetic charge g = — ^T 3 
was sitting there. 

• Exercise: Consider the gauge transformation f2' = fi {6 + 7r, <p). Show 
that it also diagonalizes <3>, such that Q'QQ t = — rT 3 and that the 
singular part of the transformed field Q! A^Vt'^ is -n y 1 ~ s -° s / ^3 which 
indicates the presence of a magnetic charge g = ^-T 3 with a string 
running along the negative 2-axis. 

4.1.2 The magnetic charges of SU (3) monopoles 

In the case of SU (3), there are two diagonal generators, namely T 3 and T 8 , 
given by ()4.5|) . The abelian gauge is the one in which the field $ = $ a T a 
acquires the diagonal form: 

$ = $ a T a -> Vim^ = aliag (A x A 2 A 3 ) A x + A 2 + A 3 = (4.21) 
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Monopoles will occur at points in space where two eigenvalues become de- 
generate. Indeed, consider the case where the first two eigenvalues are degen- 
erate: Ai = A2 = | and A 3 = —A. When the two eigenvalues Ai and A2 lie 
close to each other, the matrix $ may be considered as diagonal in all but 
the SU (2) subspace defined by the almost degenerate eigenvalues: 



1 / A + s 3 £ X - ie 2 n \ 3 
$~~ £l + ie 2 A - £ 3 ' = ^$ a T a + $ 8 T 8 (4.22) 
Z \ -2A / a=i 

$1 =e x $ 2 = e 2 $3 = £3 $8 = v^A (4.23) 

Consider the rotation (or gauge transformation) Q which brings this almost 
diagonal matrix ()4.23|) into diagonal form: 

3 / 3 \ / e 

^ g + ^$„T^(] U 8 r 8 + ^$ s r Ot = $ 8 T 8 + -e 

0=1 V o=l / \ 

(4.24) 

The rotation Q simply orients the vector ^o=i &a,T a in the T 3 direction so 
that the eigenvalue e is: 



'$2 + $ 2 + $ 2 ^ 4 25 ) 

Degeneracies of the eigenvalues will occur at points r in space where: 

$i(r ) = $ 2 (r ) = $ 3 (r ) = (4.26) 

These three equations define the three components (xo, yo, zq) of the position 
vector r*o, as in the SU (2) case. In the vicinity of the point r*o, the field $ (f) 
acquires the hedgehog shape ()4.10|) : 

$(f) = ^x / a r a + $ 8 (f)r 8 .< = ^c a 6K-x 6) c c d ® a 



.... •]. dXb r r 

(4.27) 

We work in the coordinate frame x' a and drop the primes. The degeneracy 
point is then placed at the origin of coordinates. In spherical coordinates, 
the field is: 

$ (r) = x a T a + $ 8 (r) T 8 = 7\r sin # cos<£> + T2r sin 6 sin (p+T 3 r cos^ + ^s (r) T 8 
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/ cos 9 e~ ilp sm9 
= -{ e^sinfl -cos 9 | + $ 8 (f) T 8 (4.28) 

2 \ 

The gauge transformation, which brings this matrix into diagonal form is: 

/ e^cosf sinf \ 
n (9, tp) — I - sin | e"^ cos | (4.29) 

v o i ; 

Under this gauge transformation, the gluon field becomes: 

A = A a T a -> ft + 1 ft f (4.30) 

The calculation of ftVft* proceeds exactly as in the SU (2) case and we find: 

— ftVST 1 " = - ( -e e T 2 e l,fi - e v 1 + T 3 + e„- (cosy?^ - sin</?T 2 ) J 
«e e \ rsmt/ r / 

(4-31) 

The second term becomes singular in the vicinity of the positive z-axis. The 
gauge transformed gluon field A can therefore be separated into a regular 
part A R and the singular part : 

A = A a T a = A*T a - -e,i±^T 3 (4.32) 

e rsmfc* 

Thus, in the vicinity of points where the first two eigenvalues coincide, the 

— * 

diagonal gluon A 3 feels the presence of a monopole with magnetic charge 
9 = -^T 3 . 

Consider next the case where the last two eigenvalues are degenerate: 
A2 = A3 = A and Ai = —A. When the two eigenvalues A2 and A3 lie close to 
each other, the matrix $ may be considered as diagonal in all but the SU (2) 
subspace defined by the almost degenerate eigenvalues: 

l I -2A \ 3 

$ ~2 X + £s £l ~ i£2 ) =®8t 8 + J2 Sata ( 433 ) 
\ si + ie 2 A - £ 3 / 0=1 

$1 —ei $ 2 = £2 $3 = £3 $8 = v^A (4.34) 
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where we defined: 




:/;, = - | o o 1 I . t 2 = T 7 = ^ | o o 




-i 
i 



From here we proceed as in the previous case, with the replacements T a — > t a 
for a = 1,2,3 and 8. The rotation Q, which brings the matrix 1)4.34)1 to the 
diagonal form, is: 




3 / a \ / 

J2 ®ata -► A $8*8 + ^ $ «*a ^ = $8*8 + - £ 
a=l \ a=l / V - 



It transforms the gluon field to: 

1 = A a T a -> + ^V^ fi f 

= fMfit + i ( -e e t 2 e iv - e^ 1 + CQ ^ t 3 + e„~ (cos^i - sin </?i 2 ) ] (4.37) 
e \ rsmv r J 

The second term becomes singular in the vicinity of the positive z-axis. The 
gauge transformed gluon field A can therefore be separated into a regular 
part A R and the singular part : 



?i< 1 _ l + cos6> 

A = A a T a = A*T a - -n 9 — -t 3 

e r sin 6 



IR U 1 + cosg [ 1 V3 \ 

= A a T a --n tp — - --T 3 + — T 8 4.38 

e rsmo \ 2 2 I 

where we expressed t 3 in terms of T 3 and T 8 . Again, it is only the diag- 
onal gluons which become singular in the vicinity of the positive z-axis. 
Thus, when the last two eigenvalues coincide, the diagonal gluon — |v4 3 T 3 + 
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-jpAgTg feels the presence of a monopole with magnetic charge g = — —t 3 = 

(-^3 + |T 8 ). 

In the final case where the first and third eigenvalues are degenerate, the 
gluon field transforms to: 

A = A a T a -+ A«T a + U*±?£ f \T 3 + f T.) (4.39) 

e rsintl 12 2 / 

and it is the diagonal gluon —^A 3 T 3 — ^A 8 T S which feels the presence of a 
monopole with magnetic charge g = \\T 3 + ^T 8 j. 

These results may be summarized by saying that the topological defects of 
abelian gauge fixing, in the case of SU (3), are sources of magnetic monopoles, 
with magnetic charges equal to: 

g = — K ■ T) (4.40) 

e 

where u>o=i,2,3 are the root vectors (jC.14|) of the color SU (3) group and T 
is the vector (T 3 ,T 8 ). It is this observation which suggests the form of the 
SU (3) Landau- Ginzburg (Abelian Higgs) model presented in Sect 15. II 



4.2 The maximal abelian gauge and abelian 
projection 

The choice of the field <&(x), used to fix the gauge, is far from being a 
trivial problem and the reader is referred to the 1981 paper of 't Hooft for a 
discussion of some appropriate and inappropriate choices [HI] • Many different 
choices have been used. Some are defined in terms of Polyakov loops on a 
lattice [121 12II1 1211 j some in terms of the lowest eigenvalue of the covariant 
laplacian operator 98 , some in terms of a "maximal abelian gauge" [99 . 
Most of these choices are defined on the lattice. A discussion of abelian 
gauge fixing on the lattice is beyond the scope of these lectures and we limit 
the discussion to a brief description of the maximal abelian gauge, from which 
92 % of the full string tension is obtained in the SU (2) case [TSJ. A useful 
account of evidence for the occurrence of monopoles obtained from lattice 
calculations in the maximal abelian projection, can be found in the 1997 
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Cambridge lectures of Chernodub and Polikarpov |lU(Jj . The reader is also 
referred to the Sect. 4. 10 of the extensive 2001 Physics Report of Bali |101j 
and he may find it instructive to consult the 1999 thesis of Ichie |1021 1103j 
as well as the recent 2003 paper of Chernodub [2*o] . 

Consider color SU (3). The maximal abelian gauge attempts to minimize 
the off-diagonal gluons. The gluon field can be expressed thus: 

3 

A» = A»T a = A»T 3 + A»T 8 + ^ C?E a + C»E_ a (4.41) 

a=l 

In this form, the diagonal generators T 3 and T 8 are explicit and the charged 
non-diagonal gluons C a and C* are expressed in terms of the generators E± a 
defined in (jC.8|) . Let us represent the diagonal generators T 3 and T 8 by 
the two dimensional vector H = (T3, T§). The commutator of the covariant 
derivative = + ieA%T a with the diagonal generators H i= i 2 can be 
expressed in terms of the root vectors w a , defined in (JC.14)) : 

3 

[D", Hi] =teJ2 w ai (CTE a - C»E_ a ) (4.42) 
0=1 

The commutator H] singles out the off-diagonal part of the gluon field. 
Let us calculate the trace: 

2 23 
R = trJ2 H i} [D", H i\ = 2 V ib tr (C^E a - C a ^E_ a ) {C?E h - CgE„ b ) 

(4.43) 



i=l i=l a, 6=1 



It is easy to check that: 

trE a E b = trE a E_ b = 5 ab N c (a, b) > (4.44) 
so that, using (|C.15|) . we find: 

3 

i?(a;) = 2e 2 iV c ^|^(x)| 2 (4.45) 

a=l 

We see that, in a gauge which minimizes the field R (x), the intensity of the 
charged gluons C*^ (x) and C£ (x) is minimized. 
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-Let us seek this gauge. Let Rn be the field R obtained by performing a 
gauge transformation — > QD^Q* 1 of the covariant derivative: 

2 2 
R n =trY^ [MD^, Hi] [QD^Q\ H { ] =tr^ [D m ttfHiSl] [£> M , Q^Hitt] 

i=l i=l 

(4.46) 

We want R (x) to be stationary with respect to infinitesimal gauge transfor- 
mations of the form Q ~ 1 + ix, where x = XaX a - To first order in Xi we 
have Q^HiQ = Hi — i [x, Hj\ so that the first order variation of Rq is: 

2 2 

itff = -2itr Pm> [X, P", ffi] = 2itr X £ lH, [D„ [£>", Hi]]] (4.47) 
i=i i=i 

If this is to vanish for any x — XaT a , we must have: 

[H i ,[D (t ,[D' t ,H i ]]] = (4.48) 

In the SU (2) case, the vector Hi has only one component Hi = T3. The 
maximal abelian gauge is the one which aligns the vector: 

$ = T 3 ]] (4.49) 

along the T 3 axis. 

Abelian projection in the continuum consists in making the corresponding 
gauge transformation of the gluon field and in retaining only the diagonal 
part. For example, in the expressions (j4.19|) and (|4.32|) . this means setting to 
zero the non-diagonal parts of A^T a . The monopole singular part is retained 
in this process. 

In the case of SU (3), the condition ()4.48|) reads: 

[T 3 , [£>„, [D», T 3 ]]] + [T 8 , [Dp, [D", T 8 ]]] = (4.50) 

Maximal abelian gauge fixing is more subtle in this case and the reader if 
referred to the interesting 2002 paper by Stack, Tucker and Wensley |104j . 

4.3 Abelian and center projection on the lat- 
tice. 

On the lattice the gluon field does not appear explicitly and, instead, the 
action is expressed in terms of link variables. Abelian gauge fixing and center 
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projection on the lattice is usefully reviewed in the 1997 paper of Del Debbio, 
Faber, Greensite and Olejnik [105J. In SU (2) , the maximal abelian gauge is 
the gauge which maximizes 

4 

r = E Tr fa u *,*+^ u l*+») ( 4 - 51 ) 

X [1=1 

so as to make the link variables U XtX+fl as diagonal as possible. Under a gauge 
transformation, generated by the SU (2) element Q(x) = e %Xa ^ Ta , the link 
variable transforms as: 

U x>x+Il -> fi (x) £4,*+^ (x + n) (4.52) 

and i? transforms as: 

4 

X [1=1 

(4-53) 

The maximal abelian gauge is then defined by the SU (2) element Q (x) = 
e iXa(x)T a j n -yy^i^ fa e angles a a (x) are chosen so as to maximize Rq. 

Abelian projection means the replacement of the full link variables by 
Abelian links A according to the rule 

U = oqI + ia- a -> A= = l0 (4.54) 

a/g^ + g^ V u e / 

where A stands for "Abelian" (and does not designate the gauge field). In 
these expressions, we have omitted the induces x and fi so that, for example, 
U stands for U XjX+IJ/ and 9 stands for 9 X}X+fJi . 

Abelian dominance, found by Suzuki and collaborators fTTT] , is essen- 
tially the fact that the confining string tension can be extracted from the 
Abelian link variables alone. 

The matrix U remains abelian under U (1) gauge transformations of the 
form 

/ e i<** o \ f e w "'»+" \ / e iQ *+" 

(4.55) 
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We can proceed to make a further gauge fixing by choosing the angles a (x) 
so as to maximize the quantity 

4 

]T]T cos 2 (4-56) 

X fl=l 

This defines the maximal center gauge. Of course, this still leaves a remaining 
Z2 symmetry because 9 is only determined modulo ir. The, at each link 
(x, x + fi) we can define a value of Z XiX+li as follows: 

Zxp+n = sign (cos 9 X , X+Ii ) (4.57) 

so that Z x>x+ n takes the values +1 or —1. Center projection means the 
replacement of the full link variables by Abelian links A according to the 
rule 

U = a I + ia- a —> ZI = ^ ^ (4.58) 
in the computation of observables and Polyakov lines. 
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Chapter 5 

The confinement of SU (3) color 
charges 

If we wish to describe color confinement in terms of the Meissner effect of a 
dual superconductor, we need to adapt the Landau-Ginzburg model to the 
dynamics of quarks and gluons, so as to accommodate their color quantum 
numbers. The action ()3.1|) of the Landau-Ginzburg model describes a U (1) 
gauged self- interacting complex scalar field ip. Since the magnetic current 
of the dual superconductor is somehow related to the monopoles which are 
formed by topological defects in a given gauge, as described in Sect. I4.lt h 
might make sense to restrict the covariant derivative to the corresponding 
abelian gauge. An adaptation of the Landau-Ginzburg model to color SU (3), 
which respects Weyl symmetry, was proposed in the 1989 paper of Maedan 
and Suzuki [ID] . It was further developed in the 1993 paper of Kamizawa, 
Matsubara, Shiba and Suzuki [10^] and the 1999 papers of Chernodub and 
Komarov |lU7l H08j. We shall first present the model in the absence of quark 
charges. The latter will be introduced in Sect l5.2l 
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5.1 An abelian SU (3) Landau- Ginzburg model 



5.1.1 The model action and its abelian gauge invari- 
ance 



The model action, proposed by Maedan and Suzuki has the form: 



<V, {-l(dAB 3 ) 2 -±(dABs) 2 



3 r 



0=1 



~ \(d^ a + ig K ■ B,) ^ a )\ 2 - h (i; a r a - v 2 ) 2 



(5.1) 



The first term is the kinetic term of the abelian dual gauge fields B% and 
Bq which can be grouped together to form the vector B^ = (Bi^,B^). The 
model involves three complex scalar fields ip a =i,2,3- Each scalar field ifj a is 
gauged with the abelian covariant derivative _D M = <9 M + ig {w a ■ B^) where 
w a =i,2,3 are the three weight vectors ()C.14|) of the SU (3) group: 



w 1 



;i,o) 



w 2 



1 _v / 3 
2' 2 



W3 



■14 \ ("J 



2' 2 



Thus: 



D 3 



{w 2 ■ B> 



(5.3) 

The magnetic charges appearing in the covariant derivative are assumed to 
be proportional to the weight vectors w a , as discussed in Sect. 14.1.21 The 
explicit form of the second term of the action (|5.1|) is: 



1 3 



1 2 1 

- I (9^1 +^3^1! + - 



1 

+ 2 



<9^ 3 + ig 



1 a/3 
<9^ 2 + ig I --£3 - —fig ) ^ 2 



ifi 3 + ^j-B 8 ] ^3 



(5.4) 
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The last term is a potential by means of which the three scalar fields acquire 
non- vanishing values \ip a \ — v in the ground state. Note that the last term 

does not have the form — 16 (Xlo=i — v2 ) ■ 

In this and the following sections, repeated indices i, a, ... are not assumed 
summed unless it is explicitly stated. Repeated indices of the components 
/i, v, ... of Lorentz vectors and tensors are assumed to be explicitly summed. 

It is sometimes useful to use a polar representation of the scalar fields: 



in which case the action (J5.1j) acquires the form: 



(5.5) 



I(B 3 ,B B ,S,(p) 



d x 



0=1 



~(dAB 3 ) 2 -±(dAB s ) 2 



(5.6) 



The actions and ()5.6|) have a double gauge invariance. They are 

invariant with respect to the abelian gauge transformation: 



B 3 ^B 3 + (d(3 3 ) 
i) 2 -> ^ 2 e^^3 



Vl — >• ¥>1 - /?3 V?2 

as well as to the abelian gauge transformation: 

£? 8 ^ £? 8 + (0ft) 



^3 -> V^e 



¥?2 + ~P3 <^3 



¥?3 + ~/?3 



(5.7) 



^1 -> ^1 ^2 -> V^e 



^3 -> Vfce' 



V?2 + — P8 V?3 



<^3 - "yPs 



(5.8) 



y?i <pi y?2 

We can impose a constraint on the phases ip a of the fields ip a , namely: 

<Pi + ¥2 + <^3 = (5.9) 



The constraint ()5.9|) means that the degrees of freedom of the system consist 
of the two gauge fields B 3 and £?g, the three real fields S a= i t2 ,3 and two 
phases chosen to be, for example, (p 2 and (p 3 , in which case (pi = —tp 2 — V?3- 
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The gauge transformations (|5.7|) can then be expressed in the more sym- 
metric form: 

(W a ■ B) -> (W a ■ B) + (da a ) , 1p a -> ^tfi^^a, V?a ~ «a 

a a = K ■ /?) (ai + a2 + a 3 = 0) (5.10) 

where /3 = (/3 3 , /3 8 ). 

Such gauge transformations are compatible with the abelian gauge fixing, 
described in Sect J4. 11 which are only defined modulo residual U (1) transfor- 
mations, to which the gauge transformations (J5.1U)) belong. 

The model action (J5.1|) . as all model actions so far, is QCD inspired but 
not derived. Variants to the form ()5.6j) have been proposed and studied. For 
example, in reference |lU8j . the last term is replaced by — 16 Yla=i (^1 ~ v a) ■ 
In their 1990 paper, Maedan, Matsubara and Suzuki derive an effective 
Landau-Ginzburg model assuming the existence of magnetic monopoles. The 
latter are assumed to interact minimally with the gauge field B^ L to the dual 
field tensor F^ v and a additional phenomenological interaction between the 
monopoles is postulated [33]. The partition function of the dual supercon- 
ductor can then be obtained by summing the trajectories of the monopoles in 
a Feynman path integral, using a method developed by Bardakci and Samuel 
[109 . The summation is expressed in terms of an action involving a gauged 
complex scalar field. This derivation of the Landau-Ginzburg model is also 
considered in the 2003 paper of Chernodub, Ichiguro and Suzuki JTHj. 



5.2 The coupling of quarks to the gluon field 

The color-S't/ (3) quarks may be regarded as color-electric charges embedded 
in the dual superconductor, described by the action (|5.1|) . They may be 
coupled to the system with the help of Dirac strings, as done in Sect. 13.11 
We generalize the expression (|3.3|) by adding string sources (G 3 , G$) to the 
dual abelian field tensors (F 3 ,F 8 ): 

F 3 = dAB 3 + G 3 F 8 = dAB 8 + G 8 (5.11) 

The string sources G 3 and G 8 are antisymmetric tensors which satisfy the 
equations: 

daP?=jS d a G7=3% (5.12) 
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where {j 3 ,j 8 ) are the color-electric currents. If we assume that the quarks 
couple only to the abelian gluons A 3ll and A 8 ^, then they contribute a term 
to the lagrangian, of the form: 

q [ Ta , (id* - eT 3 A 3 ^ - eT 8 A 8fl ) + m) q (5.13) 

in which case the color-electric currents jtj and j 8 would be: 

Js = ~eqit,T 3 q jg = -eq^^q (5.14) 

which agrees, of course, with the color charges (jD.10|) of the quarks. 

In the presence of color-electric charge, the action ()5.1|) is replaced by: 

/ (B 3 , B 8 , V>, r) = Jd 4 x^~(dAB 3 + G 3 ) 2 -^(dAB 8 + G 8 ) 2 

3 Pi 1 l] 

+ E ^\{d^a + ig{w a -B,)ij a )\ 2 --b{^l-v 2 ) 2 ) (5.15) 



a=l 

Let us write: 



Bg=( X 8,B 



(5.16) 

and, in analogy to (|2.90)L let us express the string terms G 3:8 in terms of 
euclidean vectors E^ and Hf* 



3,8- 



Blt,3 



_(jr 3 - 2 ^jk n st 3 - -Lr 3 - --£ Lr 3Jk 



E lt,8 



8 — 2 Lr& d k 



Hlt,s = ~Gf = —/ Hlk G 8 , k (5.17) 
The action (|5.15|) can be expressed in terms of euclidean fields, as in (|3.22j) : 



d ] , { +1 {-d t B 3 - V X 3 + - 1 (-V x B 3 + Ef 2 



1 



+- i^-d t B 8 - Vxs + #f) 2 - ^ (-V x 5 8 + 4 st ' 2 



3 r 



0=1 



1 2 1 



V^ a — ig [w a ■ B) ip a 



2\2 



--6(^:-^ 2 ) 

The fields %3 and %8 act as constraints which we do not write down. 



(5.18) 
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5.3 The energy of three static (quark) charges 



Let us calculate the energy of three static color-electric charges, which, for 
the sake of argument, we shall call quark charges. The quark charges are 
listed in the table (|D.10J) . Consider the case where three quarks, red, blue 
and green, sit respectively at the points R, B and G. Such a configuration is 
described by static color-charge densities (p 3 ,p 8 ), with: 

p 3 (f) = -e5 (r — r) e5 (f—B 



Ps (0 



-.e6 



R) + 



e5[r-B 



e5 f-G 



(5.19) 



It is illustrated on Fig l5.ll In Sect 12.111 we showed that, in the presence 
of static charges, H st = 0. The fields are time-independent and the energy 
density is equal to minus the charge density. It is simple to check that, when 
if)if)* 7^ 0, the constraints imposed by X3 an d Xs are satisfied by X3 — Xs — 0. 
The energy obtained from ()5.18|) reduces to: 



£(B,iP,ij*) = / d 3 r 



-Vx5 8 + E, 



at, 8 



a=l L 



(5.20) 



For the color-electric charges (|5.19jl . the string terms and Eg can be 
written in the form: 



E 



st,3 [r) 



1 

2^ 



B 



R 



dZ 5 [r-Z 



Est,* (r) 




dZ 5 r 



(5.21) 

where, for example, f? dZ is a line integral along a path which stems from 
the point R and terminates at the point B. The expressions (15.21)1 satisfy 
the two equations (|5.12|) . Figure IB~T1 shows two examples of such strings. 

We may proceed as in Sect. 13.3.11 and use the Ball-Caticha trick 
which consists in expressing the string terms: 



E 



; sti3 (r) = El {f) + B° 3 (f) E 



st,s (r) 



E° s (r) 



B°8 (*0 



(5.22) 
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Figure 5.1: The flux tubes formed by three static color-electric sources de- 
noted R, B and G (red, blue and green). The two color-electric charges of 
each quark are denoted in parentheses. The full and dashed lines denote 
respectively the flux tubes formed by the color-electric fields E% and E s . The 
left figure denotes the A-shaped configuration and the right figure denotes 
the Mercedes configuration. The side of the triangle is L and M = 77= is the 
distance from a summit to the center of the triangle. 



in terms of the electric fields f E%, E$j and the gauge potentials yB^, B$ j 

which are calculated in the normal (non-superconducting) vacuum, in which 
ij} a — ip* — 0. The electric fields in the normal vacuum are the usual Coulomb 
fields produced by the color-electric charges (|5.19|) of the quarks: 



£3 (f) 




1 



f-B 



E°s (r) 



V 





- + - 






f- G 


f- R 




f- B 





The gauge potentials {B®,B®\ in the normal vacuum are given by expres- 



(5.23) 



sions analogous to (13.281) : 



B 

B2 (r) = — V r x / dZ 



Ait 



R 



f— z 
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When the forms (|5.22jl of the string terms are substituted into the energy 
fl5.20j) . an expression analogous to ()3.37j) is obtained: 



8tt 



R-B 



24tt 



R-G 



24vr 



B-G 



+ J d 3 r j~ (-V xB 3 + VxB 3 j 2 + | (-V xB 8 + VxB 



a=l L ' 



(5.25) 



5.4 Quantization of the electric and magnetic 
charges 

We can repeat the argument, given in Sect. 13.3.21 to show that the energy 
(I5.25J) does not depend on the shape of the strings ()5.24|) . A deformation of 
a 3-string (the string which defines .83) corresponds to the transformation 

Bl -f Bl + z^Vfi 3 of the field This adds a term x (vfi 3 ) to the 
energy (|5.25|) which can be compensated by a singular gauge transformation 
of type ([577)1 : 

B 3 - B 3 - ^Vft 3 

V>i -> Vie ie9 ^ ^2 -> foe'" 9 ^ ^3 -> ^e - ^ (5.26) 

thereby leaving the energy ()5.25|) unchanged. As explained in Sect. 13.3.21 the 
solid angle f2 3 (r) is a discontinuous function of f. We can, however, make 
the transformed fields ()5.26|) continuous and differentiable by imposing the 
condition: 

eg = Ann (5.27) 

Similarly, a deformation of an 8-string (the string which defines B 8 ) cor- 
responds to the transformation B$ — > B$ + e S/Qs of the field B®. The 
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— * 

added gradient V^s can be compensated by a gauge transformation of 
type (|EE1): 

B 8 - B 8 - ^Vfis 

thereby leaving the energy (|5.25jl unchanged. Again, the transformed fields 
(|5.28jl become continuous and differentiable if the charge quantization condi- 
tion ()5.27|) is satisfied. The factor of 2, which distinguishes the quantization 
condition (|5.27|) from the quantization condition ()2.8U|) . is due to the fact 
that electrons have electric charge e whereas the quarks have color-electric 
charge |e. 



5.5 Flux tubes formed by the electric and 
magnetic fields 

It is convenient to express the energy (|5.25jl in terms of the polar represen- 
tation ()5.5|) of the complex scalar fields: 

2 2 2 

£(B,<p,S) = - 



57T 



R-B 



24tt 



R-G 



247T 



B-G 



j d s r |i (v x B 3 + V x B°Y + l - (-V x B 8 + V x B^ 

+ E f^r «S ■ § ) ~ + \ (^ Sa ) 2 + \ ^ ~ v2 A I (5 - 29) 

„_i L - - J ) 



We can choose to work in the gauge defined by ot a = ip a , which is usually 
referred to as the unitary gauge. In this gauge, the energy ()5.29|) reduces to: 



£{B,<p,S) 



n R-B 24n R-G 24n B-G 



r/ : V ^(vxB 3 + Vx5 3 j 2 + i (-V xBs + V x b£ 
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(5.30) 

Z V / Z \ / Z ' I 

a=l L J J 

This expression of the energy, expressed in the unitary gauge, is not indepen- 
dent of the shape of the strings, nor should it be, as explained in Sect J3.3.41 
In the unitary gauge, the flux tubes develop around the Dirac strings (|5.24|) . 
which can be chosen to be straight lines joining the charges. This choice 
presumably minimizes the energy. 

The fields which minimize the energy (J5.30j) give rise to flux tubes formed 
by the color-electric fields E 3 and E 8 , as in the Landau- Ginzburg model 
described in Sect. 13.31 The electric fields are given by the analogue of the 
expression (|2.91|) : 

E a = -fxB a + E st>a = -VxB a + VxB° a + E° a (a = 3, 8) (5.31) 

and the magnetic currents by the analogue of ()2.94|) : 

j a = -V x E a (a = 3, 8) (5.32) 

Figure I5~TI shows two possible equilibrium shapes of the flux tubes. In the 
left figure, the flux tubes join the charges thereby forming a A-shaped pattern 
of electric fields. In the right figure, the E$ flux tubes converge first towards 
the center of the triangle, thereby forming a Mercedes shaped pattern. Let 
us estimate the energy contained in the flux tubes when the quarks are far 
apart. Each flux tube, emanating from a quark gives rise to a linear confining 
potential, the intensity of which is proportional to the squared charge of the 
quark it stems from. It therefore contributes an energy which is proportional 
to the square of the charge and to the length of the flux tube. For simplicity, 
assume that the quarks are at the summit of an equilateral triangle of length 
L, in which case the distance from a summit to the center of the triangle is 
M = In the A-configuration, the color-electric field E3 contributes an 

energy proportional to £3 (A) = L (|e) 2 and the color-electric field E§ an 

energy proportional to £$ (A) = 2L ^^^ej , so that the total energy of the 
A-shaped configuration is proportional to: 

1 \ 2 , / 1 \ 2 5 



<-(A)=^(A) + &(A)=L(|eJ +2L\^j=e) = j^Le 2 = 0.417 Le 2 

(5.33) 



105 



In the Mercedes configuration, the color-electric field E 3 still forms a flux 
tube joining the red and blue quarks and contributes the same energy as 
in the A-configuration. The color-electric field fig contributes an energy 

proportional to £$ (Y) = 2M (j^ge j + M (^-^ ej , so that the total energy 
of the Y-shaped configuration is proportional to: 

£ (Y) = £ 3 (A) + £ 8 (Y) = -Le 2 + \hle 2 = -Le 2 + \^=e 2 = 0.539 Le 2 

4 z 4 2 y3 

(5.34) 

Thus the A-shaped configuration is energetically favored in this estimate, by 
about 25% of the total energy. However, in a type I superconductor, flux 
tubes attract, so that this attraction could modify the preceding estimate. 
A numerical calculation would be required to check this. Lattice simulations 
appear to favor the A-shaped configuration as long as the distance L between 
the quarks is less than 0.7 fm |110j . 

The actual flux tubes formed by the color-electric fields E 3 and Eg, ob- 
tained by minimizing the energy ()5.25|) of the Mercedes configuration are 
computed and displayed in the 1993 paper of Kamizawa, Matsubara, Shiba 
and Suzuki [106J. 



5.6 A Weyl symmetric form of the action 

We are dealing with two gauge fields and three complex scalar fields and 
this introduces an apparent asymmetry in the model, which, in fact, has the 
virtue of respecting Weyl symmetry 1 . We can obtain a form of the action, 
in which Weyl symmetry is more explicit, by expressing the two abelian dual 
gauge fields B^ = (fig , fig) in terms of the three gauge fields b^ a = (w a ■ fi M ). 
Using the completeness relations (jC.16j) or (jC17|) . we can write the vector 
fi^ = (B%,Bg)as follows: 

B»= 2 -j^w a K • B") = ~ ^ wX (5.35) 

a=l a=l 

The three fields b afM are: 

b£ = (w a ■ B») b*t + b% + b 3 l = Q (5.36) 

lr The Weyl symmetry refers to the symmetry with respect to the exchange of the three 
colors defined in the fundamental representation of the SU (3) group. 
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They are not independent because they sum up to zero. They are, in fact, 
the following linear combinations of the two fields B% and B$: 

bl = (w x ■ B») = B$ 



% = ( W2 . S M) = J B % - ^-Bg b% = (w 3 ■ B») = + ^B 



(5.37) 



We have: 



(B, ■ B") = B 3 ,B» + B 8fl Bg = -J2 KK 



o=l 
. 3 



(d A B) 2 = (d A B 3 ) 2 + (d A B 8 ) 2 = - J2(d Ab a ) 



0=1 



d^ a + ig (w a ■ B*) t/j a = d^ a + igb^a (5.38) 

Similarly, the two string terms G = (G 3 , G$) can be expressed in terms 
of three string terms g a = (w a ■ G) as follows: 



G = - ^2 w a g a 



(5.39) 



o=l 



where: 

9a = K ■ G) gi + g 2 + gz = Q (5.40) 
The action (|5.15|) can then be written in the Weyl symmetric form: 

3 



J 0=1 



1 



(dAb a + g a ) + -\d^ a + igb a ^ip a \ - -b (tp a tp* a - v 2 ) 



= I d'x £ [-~ (d A b a + g a f + 1 [dS a f + {ba + (d Va )) 2 - \b {Si - v 2 ) 2 

a=l 

(5-41) 

When variations of the action are considered, the following constraints need 
to be taken into account: 



(5.42) 



0=1 



0=1 



0=1 
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The possibly offending factors | could be changed to \ by redefining the 
fields b a and the source terms g a but we do not feel compelled to do this. 
We can obtain a Weyl symmetric form of the action f)5.18j) by defining: 

K = (la, b a ) r ]l +V2+V3 = h + b 2 + b 3 = (5.43) 
We obtain: 

I, = j d'xJ2 [ + 1 {-dtl - V Va + ft) 2 -l(-Vxb a + et) 

+ l - \{d^ a + ig Va ij a )\ 2 - \ | (fVa - ^W>a) f " ^ - « 2 ) 2 } (5-44) 

where: 

e si , a = (w« • 4rt) ^i,a = (w a ■ H st j (5.45) 

When the color-charge densities (p 3 , ps)of the quarks are given by (|5.19|) . 
the corresponding color-charge densities p a = (w a ■ p) have the remarkable 
Weyl symmetric form: 

1 „ *\ 1 



p\ if) = w\ ■ p = -eS Rj — -eS B 

p 2 (r) = w 2 - P = ~^ e ^ (f- Rj + ^e<5 G ^ 
Ps (f)=w 3 -p = ^eS (r -Bj- \eb (f - (?) (5.46) 



The sources (I5.45j) can be computed from the sources (j5.21|) : 

e s t,i (r) = r- / dZ 6 (r- Z 
Ze Jr v 

e st ,2 (r) = J dZS(f- Zj e st<3 (f) = i- J dZ 5 (f- Zj (5.47) 

which is compatible with the constraint e st) \ + e stj2 + e stj3 = 0. In the expres- 
sions above, J? dZ denotes a line integral along a path which begins at the 
point R and ends at the point B. The sources g a are strings which satisfy 
the equations: 

d ■ g a = 3a (5.48) 
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where: 



J a 



(5.49) 



In the presence of static sources p a (r), we have h% = and the fields become 
time-independent. The energy, obtained fro; (|5.44|) . is: 



/"SB 



1 / * r - \ 2 1 
+-(-Vx6 n + e si , fl ) +- 



Vtp a - igb a ip a 



1. 



2\ 2 



+ 2 & - v 2 ) 



(5.50) 

The constraints imposed by the rj a , or, more rigorously, by the independent 
fields i] 3 and rjs are satisfied with r\ a = so that the energy is: 



r 3 r 

J a=l L 



+- ( -V x b a + e sU , 



+ 



\ 2 1 2 

Vip a - igb a tp a ) + -b (ip a ip* a - v 2 ) 



(5.51) 

When we attempt to minimize the energy (|5.51|) with respect to variations 
of the fields, we must remember the constraints: 

h (f) + b 2 (r) + h (f) = ip x (f) + <p 2 (f) + <f 2 (f) = (tff a = S a e i9 ^ ) 

(5.52) 

Figure (|5.2|) shows the flux tubes which are formed by the three electric 
fields: 



ei = E 3 e 2 = — -i? 3 — 



1 /3 
e 3 = "2^3 + ^-^8 (5.53) 



as well as the string terms ()5.47j) . The figure shows two possible equilibrium 
shapes of the flux tubes. In the A-shaped pattern, the flux tubes run in 
straight lines from one color source to the other, whereas in the Mercedes 
configuration, they converge first toward the center of the triangle. Flux 
tubes running in opposite directions attract each other and this may lower 
the energy of the Mercedes configuration. Flux tubes in the Weyl symmetric 
representation are discussed in the 1998 paper of Chernodub and Komarov 

una- 
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B 3 G b G 



Figure 5.2: The Weyl symmetric representation of the flux tubes formed by 
three static color-electric sources denoted R,B and G (red, blue and green). 
The left figure denotes the A-shaped configuration and the right figure de- 
notes the Mercedes configuration. 
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Appendix A 



Vectors, tensors and their 
duality transformations 

A.l Compact notation 

Scalar and vector products of two vectors A 11 and B^ are written as: 

A-B = A^ {A A B) = A^-AvB, (d A A) = (d^)-(d u A„) 

(A ' 1} 

If S^ u = —S Ufl is an antisymmetric tensor, its contractions with a vector A^ 
are written in the form: 

(A ■ Sf = A U S U » = -S» V A U = -(S- Af (A.2) 

If S and T are two antisymmetric tensors, we write: 

S-T = ST= \s^T» v S-S = S 2 = \s^ u (A.3) 

The dot "•" is not necessary when it is clear which symbols are vectors, which 
are antisymmetric tensors, etc. In such cases we can write, for example: 

(TA)^ = T^A" ASB = A^ V B V etc. (A.4) 
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A. 2 The metric g^ u and the antisymmetric 
tensor e^ va ^ 



In Minkowski space, we use the metric: 



( 1 











\ 





-1 
















-1 







\o 








-1 


/ 



det g — — 1 



(A.5) 



The antisymmetric tensor e^ va ^ is: 



,0123 



— e 



1023 



fiua/3 



(A.6) 



When indices of e^ ual 3 are contracted, it is useful to think in terms of scalar 
products of antisymmetric "states": 



(A.7) 



e^e m , p , Y = -(aP 1 \a'f3' 1 ') = -det I g\, g\ g\ 

\ 9l> g}> g}> 

\s^E» va , p , = -(a/3\ a'P') = - det ( 9 % 9 \ ) 



e a ^e aPlu = -^\v) = -g» 



3! 



-1 



(A.8) 



The rule is that upper indices appear in the bras, lower indices appear in the 
kets. 



A. 3 Vectors and their dual form 

The metric g^ v acts as the unit operator acting on a vector A 11 : 

gA = A g^A u = A» 



(A.9) 
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A. 3.1 Longitudinal and transverse components of vec- 
tors 

We define transverse and longitudinal projectors: 

They have the properties: 

T 2 = T L 2 = L TL = LT = T + L = g A = (T + L) A (A.ll) 

if: 

A T = TA A L = LA A = A T + A L (A.12) 

Then: 

d-A = d-A L dAA = dAA T d-A T = o dAA L = o (A.13) 
A. 3. 2 Identities involving vectors 

In the following, A^ (x) , B^ 1 (x) and C M (x) are vector fields and n M is a x- 
independent vector: 

A- (BAG) = (A- B)C- (A -C)B (A.14) 

n- (d A A) = (n- d) A — d(n- A) (A.15) 

d ■ (n A A) — (n ■ d) A — n(d ■ A) (A.16) 

(A ■ WACf = -e^AaBpC^ (A. 17) 

(n ■ FXA) " = -e^ria {dpA^) (A. 18) 

(0 ■ ITAAf = e^n a (dpAy) = -n • (d~AA) (A. 19) 

d a (n A (n ■ (d A A))) a » = (n ■ df A»-{n ■ d) (n • A)-ra" (n ■ d) (d ■ A)+d 2 n» (n ■ A) 

(A.20) 

d-(dAA)=d 2 TA d-dXA = (A.21) 
d-(TA) = d-(LA) = d-A dA(TA) = dAA dA (LA) = (A. 22) 



113 



In the following identities, surface terms are neglected: 

J d 4 x (n ■ d A A) ■ (n ■ TAB) = - J d 4 x (n ■ d A B) ■ (n ■ d~A~A) 

= J d 4 x e^A^ (n ■ d) d a B p 

= - J d 4 x A(n-d)[n- dT^B] = + J d 4 xB{n-d)[n- dl^A\ (A.23) 

- J d 4 x [n • (d A A)] [n ■ (d A A)] 

= J d 4 x [A (n ■ d) 2 A - A(n ■ d) d (n ■ A) - An (n ■ d)(d ■ A) + (A - n) d 2 (n ■ A)] 

(A.24) 

J d 4 x (d A A) (d A B) = - J d 4 xA [d • (d A B)] — — J d 4 xAd 2 TB 

J d 4 x (8 ■ A) (8 ■ B) = - J d 4 xAd {d-B) = - J d 4 xAd 2 LB (A.25) 

A. 3. 3 Identities involving vectors and antisymmetric 
tensors 

Other useful identities, valid for an antisymmetric tensor S and a vector A: 

A ■ {A ■ S) = (A.26) 
L (d ■ S) = T (d ■ S) = (d ■ S) (A.27) 

J d 4 x {d AA)S = - J d 4 xA-{d-S) J d 4 x (d7^A)S=- J d 4 xA-(d-S) 

(A.28) 

A. 4 Antisymmetric tensors and their dual form 

The operator: 

Gfiu,a^ = 9^a9u/3 ~ 9nP9tiP (A. 29) 

acts as the unit operator on an antisymmetric tensor S : 

GS = S \G^S aP = S, u (A.30) 
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A. 4.1 The dual of an antisymmetric tensor 

The dual of an antisymmetric tensor S is denoted as S: 

= ^e^Sap = (sSy iu (A.31) 

We can regard e as the operator = e^ va p so as to write: 

= (eS)^ = EpvapS 01 ? (Se)^ = S al3 e a f3^ u = S (A. 32) 

e 2 = -G (A.33) 

and the duality transformation of antisymmetric tensors is reversible with a 
minus sign: 

S = eS S = - (eeS) = -eS (AM) 
If S and T are two antisymmetric tensors, we have: 

TS = -TeeS = -TS TS = TeS = TS (A.35) 
A. 4. 2 The Zwanziger identities 

An antisymmetric tensor S^ u is entirely defined by the two vectors n ■ S and 
n ■ S, where n M is a given vector. Let S and T be two antisymmetric tensors. 
If n • S = n • T and n ■ S = n ■ T, then S = T. To understand why, let us 
choose n M space-like: 

n^ = (0,n) n 2 = -n-n (A.36) 

Consider, for example, the six components of the electromagnetic field tensor 
F^ u = ^E, Hj which can be written in terms of two Euclidean 3-dimensional 

vectors E and H as in (j2.1j) and (|2.2|) . The argument is valid, of course, for 
any other antisymmetric tensor. The components of the 4-vectors n ■ F and 
n ■ F are: 

(n -F)° = -n-E (n ■ Ff = - (n x H^j 
(n-F)° = -n-H (n-Fy=(nxE^ (A.37) 

We see that n ■ F represents the longitudinal part of E and the transverse 
part of H, whereas n ■ F represents the transverse part of E and the longi- 
tudinal part of H . The vectors n x E and n x H have only two components 
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because they are orthogonal to the vector n. Therefore the eight components 
of n • F and n ■ F represent, in fact, the six independent components of the 
antisymmetric tensor F^ v . 

The first Zwanziger identity j3H| consists in writing any antisymmetric 
tensor F^ v in the form: 



F = — [n A (n- F) -nA (n -F)j (A.38) 

which is correct because the left and right hand sides give the same value for 
n ■ F and n ■ F. 

The second Zwanziger identity reads: 

F 2 = -^((n-F) 2 -(n-F) 2 ) (A.39) 

From (| A. 37(1 . we see that:: 

(n-F) 2 = (n-E\ 2 - (n x H Y (n ■ F) 2 = (H-hY - (n x 

(A.40) 

so that: 

(n- Ff-{n- F) 2 = (n-E) 2 +(nx E\ 2 -(nx H\" ' -(n ■ H\ = n 2 (e 2 - H 2 ^j 

(A.41) 

Since n 2 = —n ■ n, the Zwanziger identity ()A.39j) is verified because F 2 = 
H 2 — E 2 . 

A. 4. 3 Longitudinal and transverse components of an- 
tisymmetric tensors 

We define the projectors K^ v>a p and E^^p which are the differential opera- 
tors: | 

E = eKe E^ a/3 = ^e^ ap K^ 5 e lSa( 3 = e^ ap ^ (g^d p d 5 ) e, l&af3 (A.42) 
The projectors K and E are related by the equations: 

K 2 = K E 2 = —E KE = K -E = G (A.43) 
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We have: 

KS = ^dA{d-S) 

ES = eKeS = eKS = e— (d A (d ■ S)) = —9 A (d • S) {AAA) 

so that {K — E) S — S = GS follows from (jA~38jl . 
If a and 6 commute with K and then: 

-AT+ \e ) (aK + bE) = K- E = G (A.45) 
a b J 

In the following, 5 and T are antisymmetric tensor fields, and we neglect 
surface terms in the integrals: 

d ■ {d ■ S) = (A.46) 

d ■ {ES) = d ■ {KS) = T {d ■ S) = {d ■ S) d-S = d- {KS) (A.47) 

J d 4 x {d ■ S) {d ■ T) = — J d 4 x Sd 2 KT (A.48) 

J d 4 x (d ■ S) (d ■ f) = - J d 4 x Sd 2 ET (A.49) 

K {d A A) = d A {TA) = d A A E {d A A) = (A.50) 

K (dTwT) =0 E (BT^A) = -d~A~A (A.51) 

The projectors K and E can also be defined in terms of a given vector n M . 
One simply replaces <9 M by n M . For example, the projector K can be defined 
thus: j 

^,a/3 = -j {g^n v n p - g^n^np + g v ^n a - g^n u n a ) 

and we still have: 

E = eKe K 2 = K E 2 = —E KE = K — E = G (A.52) 

The identity K — E = G is a statement of the Zwanziger identity (|A.38|) . 
Indeed, we have: 

KS = —n A (n • 5) 
n 2 

KS = \nA{n- S) ES = eKS = e— (n A (n • S)) = —n A (n-S) 
n 2 n 2 n 2 

(A.53) 

so that (K — E) S = S = GS is a statement of (jA.38|) . 
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A. 5 Antisymmetric and dual 3-forms 

Let T a/ g 7 = —Tp ai = ... be a tensor which is completely antisymmetric with 
respect to the exchange of its indices. We define: 

T 2 = l -T aPl T a ^ (A.54) 

there being 3! = 6 identical terms in the sum. The dual of the tensor T 01 ^ 1 
is the vector defined thus: 

T» = l -e^T aPl (A.55) 

and we have: 

f 2 = — T 2 (A.56) 

which, in explicit form, reads: 

-e^T^-e^jvrW = -~T a/37 T^ (A.57) 



This result can be checked using (|A.8|i . 

In the particular case where the tensor T is defined in terms of the deriva- 
tive of an antisymmetric tensor $ A " / : 

d a <$> Pl + dpQ^a + <9 7 $ Q/3 (A.58) 

we have: 

?" = -e^T aM = - («9 • <5) M (A. 59) 
6 

so that: 

br 2 = - ^f 2 = - ~ (d ■ $) 2 (A.60) 
Note also that, neglecting surface terms: 

J d ± x l T 2 = -jd A x^(d-$) 2 = ^J d A x (~ (d • $) 2 - U>d 2 § 

(A.61) 

We can also define the dual of a vector V 1 to be the antisymmetric tensor 
T a/37 defined as: 

yafr = e *fhny (A.62) 



118 



and we have: 

V 2 = -V 2 l -V a ^V aM = -VVp (A.63) 
Note also that, neglecting surface terms: 

J ±T 2 = -Jd A x^(d-$) 2 = ^J d A x (-± (d • $) 2 - Ud 2 ^j 



Note that, if: 
then: 

However, if: 



(A.64) 

<5> = dAA (A.65) 
T a/37 = f » = (A.66) 



<S> = dAA (A.67) 

then: 

= d a d A Apy + d d A A ia + 8^8 A A a/3 

= 8 a e Ml2 (PA 2 ) + dpe ial2 (d x A 2 ) + (d l A 2 ) (A.68) 
and the dual of T is: 

A. 6 Three-dimensional euclidean vectors 

Formulas are taken from [TF] . 

The cartesian components are 2 = (x, y, z) and the 3-dimensional anti- 
symmetric tensor is: 

£123 = -£213 = ... = 1 (A. 69) 

The cartesian component i = (x, y, z) of a three-dimensional vector a is 
denoted by aj. 

a ■ (b x cj — b- (c x a) — c- (a x b ^ = EijuaibjCk (A. 70) 

ax ^0* x cj — (a ■ c)b — (^a ■ bj c (A. 71) 

(axb\ ■ (cxdj =(a-c) (b ■ dj - (a-dj (fc • c) (A.72) 
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V x = (A.73) 
V-(vxo)=0 (A.74) 

V x (Vxfl)=v(faj-V 2 fl (A.75) 

V- (ipa) = a- V^ + ^V-a (A.76) 

V x (ipa) = xa + fVxo (A.77) 

V (a-tj = (a- v) b+ (b ■ v) a + ax (V xi)+6x|vx a) (A.78) 

V-(ax&)=6-(vxa)-a(vxft) (A.79) 

V x (a x fc) = a (v-&) -fc(V-a) + (fc-V)a- (a-V)fc (A.80) 
We define longitudinal and transverse projectors: 

^" = ^ T^8 i3 -^ L* = L T 2 =T LT = (A.81) 
A vector field may be decomposed into longitudinal and transverse parts: 

A = A L + A T = ^V(V-A) -ivx (Vxl) (A.82) 

A L = i,V(V-/) I T = i-iv(v-l)=-ivx(vxl) 

(A.83) 

y d 3 r (v -A} =- J d 3 r A L V 2 A L J d 3 r (v x A} = - J d 3 r A T V 2 A T 

(A. 84) 

For the unit vector e — we have: 

V-f=3 Vxf=0 (A.85) 

V-e = - Vxe = ( e = - ) (A.86) 



r \ r 



A- V)e = - U-eU-e)) = ^ (A.87) 
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J d 3 r (v x Aj-B = J d 3 r A - (V x (A.8£ 



J d 3 r(v x A) - (y x A) = J, 



V x V x A 



V x .4 ] = / d 3 r A 

— » / . , — » — * 

A ■ = / d r V • A (divergence theorem) 
Jv 

I \pds= [ d 3 r Vi) 
Js Jv 

[ dsx A = [ d 3 rV x A 
Js Jv 

A ■ dl (Stoke' s theorem) 

c 

J dsxV^J = / i)dl 



V x A) -ds 



(A.89) 

(A.90) 
(A.91) 
(A.92) 
(A.93) 
(A.94) 



is Jc 

Note that in four dimensions, A A B is a six-component antisymmetric 
tensor, whereas in three dimensions, a x b is three component vector. That 
is why the four-dimensional identity A ■ (B A C) = (A ■ B) C - (A ■ C) B 

plays the role of the three-dimensional identity ax (j) x = (a ■ c)b — 

(a ■ c. That is also why a four-dimensional vector A^ is represented by the 
identity A = 4j [n ■ (n A A) + (n ■ A) n], whereas a three dimensional vector 
A is represented by the identity A = (n ■ Aj n — n x (n x /)]. 

A. 6.1 Cartesian coordinates 

In cartesian coordinates, the position vector is f = (x,y,z). The three unit 
vectors e^ =x ^ y>z ) are: 



e'(i) = y e(i) • e(j) = <^ e(j) x e^-) = e ijk e^ k ) (A.95) 
Any vector A can be expressed in terms of its cartesian components A i=x ^ y y. 



— * 

A = A x c( x) + A y c( y ) + A z e^ z ) 



(A.96) 
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Then: 

_^ , _ dip _ dip _ dip 

W = e ( ^ + e (?/) - + e (2) - (A.97) 
- 7 dA x dA v dA z 

V-A=—^ + —^ + —^ A.98 

dx dy dz 

VxA-e ( 9Az - dA A e ( dA *_ dA A g fdA y _dA x 
^ \ dy dz J (y) \ dz dx J (z) \ dx dy 

(A.99) 

W 2 ip + + (A 100) 

<9x 2 dy 2 dz 2 

A. 6. 2 Cylindrical coordinates 

In cylindrical coordinates, the position vector is r = (p, 6, z) where the 
cylindrical coordinates are expressed in terms of the cartesian coordinates 
r— (x,y,z) as follows: 

x = pcos6 y = psinO z = z (A. 101) 

The three unit vectors e^ =Pt e iZ ) are defined by: 



//»oo /»oo /' 

d 3 r = / pdp / dz 
Jo J-oo Jo 



2tt 

de 



%) ' %) = ^3 %) x %) = £ yfce( fc ) (i, j = p, 0, z) (A. 102) 

Any vector A can be expressed in terms of its cylindrical components A i=Pt g iZ : 

A = A p e (p) + A e e {e) + A z e {z) (A. 103) 

Then: 

dip 1 dip dip 

V1, = ew J!L + em +t0l) JSL (a.104) 



v ' a =m; ( ^ )+ p1» (A - 105) 



= e ^ ( p ^ - ^7 J ~ W ) +e{z) ~p (dp- {P e) " ~de 

(A.106) 
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1 d ( dty\ 1 d 2 ip : d 2 ip 



^ = Wp\!%) + 7W + B? <A ' 107) 

poo p2n poo 

/ pdp / d9e^ cose f (p) = 2vr / pdpj (kp) f (p) (A.108) 
Jo Jo Jo 

A. 6. 3 Spherical coordinates 

In spherical coordinates, the position vector is r = (r, 9, </?) where the spher- 
ical coordinates are expressed in terms of the cartesian coordinates f = 
(x,y,z) as follows: 

x = r sin 9 cos ip y = r sin 9 sin if z = rcos9 (A. 109) 

The three unit vectors e^ =Pt g !ip ) are defined by: 

d 3 r = / r 2 dr / sin 9 d9 d(p 
Jo Jo Jo 

Hi) ■ = Sij e(i) x = e ijk e^ k ) = r, 9, ip) (A. 110) 

— * 

Any vector A can be expressed in terms of its spherical components A i=r fi jlf : 

A = A r e^ r ) + A e e {e) + A^e^ (A. Ill) 

Then: 

- _ dib _ 1 dtp _ I dip /a tt c\\ 

W = e <" * + e ^rW + ^T^edj (A - U2) 

V • ~ (r 2 A r ) + - J- A {sin 9 A e ) + (A.113) 

^7 - 1 / ' 9 , ■ n , ^ 9A e \ _ ( 1 0A r 19. .A _ 1 

(A.114) 

- 2 ( 1 9 / 1 d f . n dtj}\ 1 <9 2 V> 

V 2 ^ = -^tt: r 2 -^ + _ — — m*9-£ + — ; tt.~t~t A.115 



r 2 <9r \ dr / r 2 sin 9 d9 \ d e J r 2 sin 2 cfy? 2 
r 2 <9r V <9r / r <9r 2 
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Note the identity: 

e ^=- v n^^v (A - 117) 

y d 3 re i *' f> / (r) = 4tt r 2 drj (fcr) / (r) = 4tt r 2 dr S -^-f (r) 

(A.118) 
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Appendix B 



The relation between 
Minkowski and Euclidean 
actions 

The Minkowski action leads to canonical quantization and it is used to 
calculate matrix elements of the evolution operator e~ lHt . The Euclidean 
action is used to calculate the partition function tre~^ H . Lattice calcula- 
tions are formulated in terms of the Euclidean action. In Minkowski space 
9\iv — (1, — 1, — 1, — 1) and detg = —1, whereas in Euclidean space g^ = 
(1, 1, 1, 1) = 5^ and det g — +1. As a rule of the thumb, a Euclidean action 
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can be transformed into a Minkowski action by the following substitutions: 



Euclidean 




Minkowski 






5V — ctiag [l, — i, —i, —i) 


[t,r) 




(it, r) 


^t; Ui) — \yti V l) 




(—ipl, Ft.\ — (—nFl, —XJ\ 
\ lUt,Ui) — ( ICJt, V i) 


(A ,Ai) 





(iA Q ,A l ) 


1 7n. 7 1 




1 % in . 7 1 


A A 




At 1 A 


p,2 _ pi2 1 n2 _ Q2 , Q2 




Pl'2 a p)[L 

—a — — c M a^ 


f a a A\ Pi A Pi A 




n ( Pi a A\^ 1 „' ^AO /ti A* /t0\ 
— Z ^<7 A AJ = — l [O A — a A ) 


fa a a\ Pi A Pi A 

[O A A) tj — OiAj — OjAi 




(Ft A F)i A7 J_ Fti A* 

— [O A A) — — C /i J + u J A 


1 (Ft A 4^ 




1 (Ft A A\ (Ft A A\^ v 








c- _ c 0123 _ 1 
&0123 — e — 1 




^0123 _ _£: 0123 = l 


dAA 0t = ±eoi jk (d AA) jk 




-I e o»J*(9A^). fc = -9a/ 


d A Aij = e ij0k (d A A) ok 




z9 A A 3 = ie l ^ k (d A A) ok 


±8 A A^d A Aft,, 




-±dAA IJV dAA'~ 


J d 4 x = J dtdxdydz 




j d A x = j dtdxdydz 


{action) 




— (action) 


E 




-IE 


— * 

H 




—i 

H 



(B.l) 



For example, the Minkowski Landau-Ginzburg action ()3.8|) of a dual super- 
conductor is: 

h (B, S,<p) = J d'x (-\ (d A B + G) 2 + 9 -^f (B + d^f + i (OS) 2 - U (S 2 - t, 2 ) 2 ) 

(B.2) 

whereas the Euclidean action is: 

h (B, S,<p) = jd A x Q (dAB + G) 2 + 9 -^f(B- d^f + \ (dSf + \b (5 2 - , 2 ) 2 ) 

(B.3) 

The table (jB.ljl can be used to recover the Minkowski action ()B.2jl from the 
Euclidean action ()B.3j) . The change in sign of the action is chosen such that 
the partition function can be written in terms of a functional integral of the 
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Euclidean action, in the form: 



J D(B,S,<p)e 



(B.4) 



In general however, the functional integrals need to be adapted to the acting 
constraints.. 

We can choose to represent the Euclidean field tensor F^ v = in terms 
of Euclidean electric and magnetic fields E and H thus: 



( 

E T . 



-E x 




E„ H z 



-E y 




\ E z —H y H x 



-E z \ 

Hy 

-Hx 
1 



— 1 



2 ^ nuafiFafi 



( o 

Hy 



-H x 


-E z 

E„ 



-Hy 

E z 




(B.5) 

If we want to express the Euclidean field tensor as F = d A A then 
the relation between the Euclidean and Minkowski electric and magnetic 
fields is the one given at the end of table IB.ll The Euclidean electric and 
magnetic fields E and H are expressed in terms of the Euclidean gauge 



potential A, 



A I as follows: 



E = -d t A + V0 H = -V x A 



(B.6) 



In the Euclidean formulation, e 2 = G and the duality transformation of 
antisymmetric tensors is reversible without a change in sign: 



Su.v — ~£ u.unftSnft S„.v — —& iLvrrftSnft (B.7) 

The projectors K and E are defined by (jA.42j) with g^ u = and we have: 



K 2 = K = eKe E 2 = E KE = K + E = G (B.8) 



with: 



(B.9) 



-Ey 

E x 
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Appendix C 

The generators of the SU (2) 
and SU (3) groups 



C.l The 577(2) generators 

The three Pauli matrices are: 



°'={\\) -"-(r-ij ™ 

o { Oj = Sij + ie ijk ai (C.2) 
The three generators of the SU (2) group are T a = \o a . 

C.2 The SU (3) generators and root vectors 

The eight Gell-Mann matrices are: 



A 2 = \ i A, 







A 4 = A n 
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A 6 = | 1 | A 7 = | -i 

i 

2 11 

Kh = ^ab + {dabc + ifabc) A c T a T b = -5 ab + - (d abc + if abc ) T c (C.4) 

tr T a T b = -5 ab T a T b — T b T a = if abc T c T a T b + T b T a = -5 ab + d abc T c 

(C.5) 

where T a = |A a are the generators of the SU (3) group. There are two 
diagonal generators, namely T 3 and T 8 . They are said to form a Cartan 
subalgebra in SU (3). 

The f abc are antisymmetric in their indices /i 23 = — /213 = ••• and the non 
vanishing values are: 

1 

/l23 = 1, fl47 = /l65 = ^246 = ^257 = ^345 = ^376 = ^' / 4 58 = ^678 = 

(C.6) 

The d abc are symmetric in their indices d\ 2 z = d 2 \ 3 = ••• and the non vanishing 
values are: 



1 1 

^118 = ^228 — ^338 = ~^=, ^146 = ^157 = ^256 = ^344 = ^355 = ^ 



1 1 

^247 = ^366 = ^377 = ~Zi <^448 = ^558 = ^668 = ^778 = ^888 = ^118 



2 5 ""Via "'OOS DDo "'liO "-BOO 2^/3 

(C.7) 

The remaining six non-diagonal generators can be grouped together as 
follows: 

E± x = ±={T x ±iT 2 ) E ±2 = ^=(T 4 T^T 5 ) E ±3 = -±= (T 6 ± iT 7 ) (C.8) 
The commutators with the diagonal generators are: 

[T 3 , E ±1 ] = ±E ±1 [T 3 , E ±2 \ = t\e ±2 [T 3 , E ±3 ] = t\e ±3 

[T 8 ,E ±1 }=0 [T 8 ,E ±2 ] = T } yE ±2 [T 8 , E ±3 ] = ±^-E ±3 (C.9) 
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The gluon field can be expressed thus: 

3 

A» = A»T a = A%T 3 + A%T 8 + ]T C?E a + C£E_ a (CIO) 



0=1 

where the non-diagonal gluon fields are: 

CT = ^ K + C£ = i= K - iK) C% = ^ (A + Mtj 

C? = -L K - Cf = ^= K + <T - ^= K - itf) 

(C.ll) 

C.3 Root vectors of 5C7 (3) 

We can represent the two diagonal generators T 3 and T 8 by a two-dimensional 
vector H: 

H=(T 3 ,T 8 ) (C.12) 

The commutators (jC.9|) of if with the non-diagonal generators can be ex- 
pressed in the form: 

[H, E ±a ] = ±w a E ±a (a = 1, 2, 3) (C.13) 

The vectors w a are called the root vectors of the SU (3) group. Their com- 
ponents are: 

»i=(i,o) -»=(4-f) -*=(4t) < c ") 

The root vectors have unit length. They are neither orthogonal nor linearly 
independent, because they sum up to zero: 

(wi ■ wi) = 1 (w 2 - w 2 ) = 1 (w 3 ■ w 3 ) = 1 wi+w 2 +w 3 = (C.15) 

They form an over-complete set: 

X>X = |<% (C16) 

o=l 

130 



where w l a is the component % of the root vector w a . It can be useful to think 
in terms of bras and kets and to write w l a = (i\ w a ) = (w a \i) with = 5y. 
In this notation, the completeness relation (jC.16|) reads: 



2 3 



-J2M(w a \ = l (C.17) 



a=l 



The abelian projection of the gluon field is: 

= A 3 ^T 3 + A 8fl T 8 = (A„ -if) A^ = (A 3fl , A 8fl ) (C.18) 

In view of d03i, we can write: 



2 3 



i^-V^^.iJ (C.19) 

6 a=l 

and the abelian projection of the gluon field is then: 

2 3 2 3 

A M • if = - V (if ■ w a ) (w a -A (t ) = -J2 a ^ta (C.20) 

a=l a=l 

where the generators £ a are: 

t a = (H- w a ) h + t 2 + t 3 = 

*i = (if • = T 3 

t 2 = (F.^ 2 ) = _lT 3 -^T 8 t 3 = (H-w 3 ) = -^T 3 -^-T 8 (C.21) 
and where the fields a afl are: 

<V = K • A*) + fl 2M + = 

ai„ = | (™i • A„) = (C.22) 

2 If r \ 2 1/ 

fl2 M = g (w 2 • Afjt) = - y-A 3f , - V3A 8m J a 3fM = - (w 3 ■ A J = - y-A 3fl + 

(C.23) 
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Appendix D 

Color charges of quarks and 
gluons 

D.l SU (2) color charges 

In SU (2), the charge operator is: 

= ^4(o -e) < D1 ' 
The quarks form a doublet in the fundamental representation of SU (2): 

The red and blue quarks have respectively color charges and — |e. 
The gluon field has the form (|C.10|) : 

A» = A%T a = A%T 3 + C^Ex + CfE_i (D.3) 

and the color charges of the gluons are given by the commutator: 

[Q, A»] = e [T 3 , A»] = eC?Ex - eCf E- X (D.4) 

The color charge of A3 is zero and the fields C±* and Cf have respectively 
color charges e and — e. 



9b 


?b 








e 
2 


e 
2 





e 


— e 



(D.5) 
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D.2 SU (3) color charges 

In SU (3) there are two color charge operators, associated respectively to the 
diagonal matrices T 3 and T 8 : 



Q 3 = eT 3 = - -e Q 8 = eT 8 





2^ 

(D.6) 

The charge operator can be represented by the vector Q = (Q3,Q$) = 
e(T 3 ,T 8 ). 

The quarks form a triplet in the fundamental representation of SU (3): 

Qr \ ( Qr\ i ( Qr 





Qs | Qb I = 77 e ~QB Qs \ Qb \ = —j=e q B 

2 V / \Qo) 2VS \-2g G 

(D.7) 

The red, blue and green quarks have respectively color charges Q 3 equal to 
|e, — |e and and color charges Qs equal to ^^e, and — -^e. 
The gluon field has the form (|C.10|) : 

3 

A" = A»T a = A%T 3 + ( C a* E a + C aE-a) (D.8) 

a=l 

We calculate the commutators of the charge operator ()D.6jl written in the 
form (lUlOl) : 

3 

[Q, A»}=eJ2 «>« (CT J5« - C^_ a ) (D.9) 

0=1 

where w a are the root vectors (jC.13|) . The color charges of A 3 and A 8 are 
zero and the fields C%* and C^ 1 have respectively color charges ew a and — ew a . 
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cr 
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cr 




cr 
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(D.10) 



A model which confines only color charges will not confine the diagonal 
gluons which have zero charge. 



133 



Acknowledgments 

I wish to thank professor Wolfram Weise for inviting me to write up these 
lectures, for constructive criticism and for constant encouragement at times 
when it was needed most. I thank Nicolas Wschebor for illuminating dis- 
cussions. I also benefited from numerous exchanges with the physicists at 
ECT*, in particular from discussions with Pietro Faccioli and Luca Girlanda 
and Evgeni Kolomeitsev. 



134 



Bibliography 



[1] H. B. Nielsen and P. Olesen. Vortex-line models for dual strings. Nucl. 
Phys. B61, page 45, 1973. 

[2] Y.Nambu. Strings, monopoles, and gauge fields. Phys. Rev. D10, page 
4262, 1974. 

[3] M. Creutz. Phys. Rev. D10, page 2696, 1974. 

[4] G. T'Hooft. High Energy Physics. Editorice Compositori Bologna, 
1975. 

[5] G. Parisi. Phys. Rev. D10, page 970, 1975. 

[6] A. Jevicki and P. Senjanovic. String-like solution of the Higgs model 
with magnetic monopoles. Phys. Rev. Dll, page 860, 1975. 

[7] S.Mandelstam. Phys. Rep. C23, page 245, 1976. 

[8] V.L.Ginzburg and L.D.Landau. Zh.Eksp. Teor.Fiz. 20, page 1064, 1950. 
An English translation may be found in "Men of Physics: L.D.Landau" , 
vol.I,p.l38. 

[9] P.A.M. Dirac. Phys. Rev. 74, page 817, 1948. 

[10] Shinji Maedan and Tsuneo Suzuki. An infrared effective theory of quark 
confinement based on monopole condensation. Progr. Theor. Phys. 81, 
page 229, 1989. 

[11] Pietro Faccioli and Thomas A. DeGrand. Evidence for instanton in- 
duced dynamics from lattice QCD. 2003. hep-ph/0304219. 



135 



[12] Tsuneo Suzuki and Ichiro Yotsuyanagi. Possible evidence for abelian 
dominance in quark confinement. Phys.Rev. D42, page 425, 1990. 

[13] S.Kiura Y.Matsubara O.Miyamura S.Ohno S.Hioki, S.Kitihara and 
T.Suzuki. Abelian dominance in SU(2) color confinement. Phys.Lett. 
B272, page 326, 1991. 

[14] Hirochi Shiba and Tsuneo Suzuki. Monopoles and string tension in 
SU(2) QCD. Phys.Lett. B333, page 461, 1994. 

[15] Steven D.Neiman John D. Stack and Roy J.Wensley. String tension 
from monopoles in SU(2) lattice gauge theory. Phys.Rev. D50, page 
3399, 1994. hep-lat/9404014. 

[16] V.Bornyakov M.Mueller-Pruessker G, Bali and K. Schilling. Dual su- 
perconductor scenario of confinement: A systematic study of Gribov 
copy effects. Phys.Rev. D54, page 2863, 1966. hep-lat/9603012. 

[17] E.-M.Ilgenfritz-T. Suzuki Y.Koma, M.Koma and M.I.Polikarpov. Du- 
ality of gauge field singularities and the structure of the flux tube 
in abelian-projected SU(2) gauge theory and the dual abelian Higgs 
model. 2003. hep-lat/0302006. 

[18] Katsuya Ishiguro M.N. Chernodub and Tsuneo Suzuki. Determina- 
tion of monopole condensate from monopole action in quenched SU(2) 
QCD. 2003. hep-lat/0308023. 

[19] L. Montesi A. Di Giacomo, B. Lucini and G. Paffuti. Colour confine- 
ment and dual superconductivity of the vacuum -. Phys. Rev. D61, 
page 034503, 2000. hep-lat/9906024. 

[20] L. Montesi A. DiGiacomo, B. Lucini and G. Paffuti. Colour confine- 
ment and dual superconductivity of the vacuum - II. Phys. Rev. D61, 
page 034504, 2000. hep-lat/9906025. 

[21] A. DiGiacomo B. Lucini Anf G. Paffuti J.M. Carmona, M. D'Elia. 
Color confinement and dual superconductivity of the vacuum - III. 
Phys. Rev. D64, page 114507, 2001. hep-lat/0103005. 

[22] Paolo Cea and Leonardo Cosmai. A gauge invariant study of the 
monopole condensation in non abelian lattice gauge theories. 2000. 
hep-lat/0006007. 



136 



[23] Paolo Cea and Leonardo Cosmai. Abelian monopole and vortex con- 
densation in lattice gauge theories. JEEP 0111, page 064, 2001. 

[24] Paolo Cea and Leonardo Cosmai. Abelian chromomagnetic fields and 
confinement. JEEP 0302, page 031, 2003. hep-lat/0204023. 

[25] A. Di Giacomo and C. Paffuti. The abelian projection revisited. 2003. 
hep-lat/0309019. 

[26] M.N. Chernodub. Monopoles in abelian polyakov gauge and projection 
(in) dependence of the dual superconductor mechanism of confinement. 
2003. hep-lat/0308031. 

[27] L.Polley and U.-J. Wiese. Monopole condensate and monopole mass 
in U(l) lattice gauge theory. Nucl. Phys. B356, page 629, 1991. 

[28] L.Polley M.J.Polikarpov and U.-J. Wiese. The monopole constraint ef- 
fective potential in U(l) lattice gauge theory. Phys. Lett. B253, page 
212, 1991. 

[29] Alan H.Guth. Existence proof of a nonconfining phase in four dimen- 
sional U(l) lattice gauge theory. Phys. Rev. D21, page 2291, 1980. 

[30] D.A. Browne V. Singh and R.W. Haymaker. Structure of Abrikosov 
vortices in SU(2) lattice gauge theory. Phys. Lett. B306, page 115, 
1993. 

[31] G.S. Bali K. Schilling and C. Schlichter. A Ginzburg-Landau analysis 
of the colour electric flux tube. Prog. Theor. Phys. Suppl. 131, page 
645, 1998. hep-lat/9808039. 

[32] G.S. Bali. The mecanism of quark confinement. Newport News 1998, 
Quark confinement and the hadron spectrum, 1998. hep-ph/9809351. 

[33] Christoph Schlichter Gunnar S. Bali and Klaus Schilling. Probing 
the QCD vacuum with static sources in maximal abelian projection. 
Progr. Theor. Phys. Suppl. 131, page 645, 1998. hep-lat/9802005. 

[34] M.I. Polikarpov F.G. Gubarev, E.-M. Ilgenfritz and T. Suzuki. The 
lattice SU(2) confining string as an Abrikosov vortex. Phys. Lett. B468, 
page 134, 1999. hep-lat/9909099. 



137 



[35] Paolo Cea and Leonardo Cosmai. Dual superconductivity in the SU(2) 
pure gauge vacuum: a lattice study. Phys.Rev. D62, page 094510, 1995. 
hep-lat/9504008. 

[36] J. Giedt J. Greensite L. Del Debbio, M. Faber and S. Olejnik. Detection 
of center vortices in the lattice yang-mills vacuum. Phys. Rev. D58, 
page 094501, 1998. hep-lat/9801027. 

[37] J. Greensite. The confinement problem in lattice gauge theory. 
Progr.Part.Nucl.Phys. 51, page 1, 2003. hep-lat/0301023. 

[38] J. Giedt J.Ambjorn and J. Greensite. Vortex structure vs. monopole 
dominance in abelian projected gauge theory. JHEP 02, page 033, 
2000. hep-lat/9907021. 

[39] S.N.Syritsyn A.V. Kovalenko, M.I. Polikarpov and V.I. Zakharov. In- 
terplay of monopoles and P- vortices. 2003. hep-lat/0309032. 

[40] M.I. Polikaprov V.G. Bornyakov, A.V. Kovalenko and D.A. Sigaev. 2q 
and 2q2q systems in terms of p- vortices. 2003. hep-lat/0309049 (Talk 
presented by M.I. P. ant Lattice 2003, Tsukuba). 

[41] Attilio Cucchieri and Daniel Zwanziger. Static color-coulomb force. 
Phys. Rev. Lett. 78, page 3814, 1997. hep-th/9607224. 

[42] Daniel Zwanziger. Lattice coulomb hamiltonian and the static coulomb 
field. Nucl.Phys. P>485, page 185, 1997. hep-th/9603203. 

[43] Adam P. Szczepaniak and Eric S. Swanson. Coulomb gauge QCD, 
confinement and the constituent representation. Phys. Rev. D65, page 
025012, 2002. hep-ph/0107078. 

[44] Jeff Greensite and Charles B. Thorn. Gluon chain model of the confining 
force. JHEP 02, page 014, 2002. hep-ph/0112326. 

[45] Jeff Greensite and Stefan Olejnik. Coulomb energy, vortices and con- 
finement. 2003. hep/lat-0302018. 

[46] P.A.M. Dirac. Proc. Roy. Soc. (London) A133, page 60, 1931. 

[47] John David Jackson. Classical Electrodynamics. John Wiley and Sons, 
second edition, 1975. 



138 



[48] Daniel Zwanziger. Local-lagrangian quantum field theory of electric 
and magnetic charges. Phys. Rev. D3, page 880, 1970. 

[49] Alfred S. Goldhaber. Role of spin in the monopole problem. Phys. Rev. 
B140, page 1407, 1965. 

[50] R. Jackiw. Dirac's magnetic monopoles (again), hep-th/0212058, 2002. 

[51] W. Meissner and R. Ochsenfeld. Naturwiss. 21, page 787, 1933. 

[52] Steven Weinberg. The Quantum Theory of Fields. Cambridge Univer- 
sity Press, 1986. 

[53] A. Vilenkin and E.P.S. Shellard. Cosmic Strings and Other Topological 
Defects. Cambridge University Press, 1994. 

[54] James S. Ball and Ariel Caticha. Sperconductivity: A testing ground 
for models of confinement. Phys. Rev. D37, page 524, 1988. 

[55] Yoshimi Matsubara Shinji Maedan and Tsuneo Suzuki. Abelian con- 
finement mechanism and the QCD vacuum. Prog. Theor. Phys. 84, 
page 130, 1990. 

[56] H.G.Dosch M.Baker, N.Brambilla and A.Vairo. Field strength correla- 
tors and dual effective dynamics in QCD. Phys. Rev. D58, page 034010, 
1998. hep-ph/9802273. 

[57] A. A. Abrikosov. On the magnetic properties of superconductors of the 
second group. Sov. Phys. JETP 5, page 1174, 1957. 

[58] E.B. Bogomolnyi. Sov. J. Nucl. Phys. 24, page 449, 1976. 

[59] Shinji Ejiri Yoshimi Matsubara and Tsuneo Suzuki. The (dual) Meiss- 
ner effect in SU(2) and SU(3) QCD. Nucl. Phys. Proc. Suppl. 34, page 
176, 1994. hep-lat/9311061. 

[60] E;-M. Ilgenfritz Y. Koma, M. Koma and T. Suzuki. A detailed study of 
the abelian-projected SU(2) flux tube and its ginzburg-landau analysis. 
2003. hep-lat/0308008 vl. 

[61] H.Toki and H.Suganuma. Dual Ginzburg-Landau theory for confine- 
ment and chiral symmetry breaking. Progr. Part. Nucl. Phys. 45, page 
S397, 2000. 



139 



[62] L.N. Cooperl956. Phys. Rev. 104, page 1189, 1956. 

[63] L.N. Cooper J. Bardeen and J.R. Schrieffer. Phys. Rev. 108, page 1175, 
1957. 

[64] T.D.Lee. Particle Physics and Introduction to Field Theory. Harwood 
Academic Press (New- York), 1981. 

[65] H. Toki and H. Suganuma. Dual Ginzburg-Landau theory for confine- 
ment and chiral symmetry breaking. Progr. Part. Nucl. Phys. 45, page 
S397, 2000. 

[66] S.Sasaki H. Suganuma, H.Ichie and H.Toki. In Color Confinement and 
Hadrons, page 65. World Scientific, 1995. 

[67] A.Hart and M.Teper. Phys. Lett. B371, page 261, 1996. 

[68] K.Originos R.Brower and C.Tan. Phys. Rev. D55, page 6313, 1997. 

[69] H. Suganuma A. Tanaka- H.Toki M. Fukushima, S.Sazaki and 
D.Diakonov. Phys. Lett. B399, page 141, 1997. 

[70] J. Z inn- Just in. Quantum Field Theory and Critical Phenomena. Claren- 
don Press, Oxford, 1989. 

[71] J.P.Blaizot and G.Ripka. Quantum Theory of Finite Systems. MIT 
Press, 1986. 

[72] John W.Negele and Henri Orland. Quantum Many-Particle Systems. 
Perseus Books, Reading, Mass, USA, 1998. 

[73] A.Di Giacomo and H.Panagopoulos. Gauge invariant field correlators 
in QCD at finite temperature. Nucl. Phys. B>483, page 371, 1997. hep- 
lat/9603018. 

[74] A. Di Giacomo M. D'Elia and E.Meggiolaro. Field strength correlators 
in full QCD. Phys. Lett. B4O8, page 315, 1997. hep-lat/9705032. 

[75] Michael Kalb and PRamond. Classical direct interstring action. Phys. 
Rev. D9, page 2273, 1974. 



140 



[76] Peter Orland. Frustration and dual superconductivity in lattice gauge 
theories. Proceedings of the Argonne conference Gauge Theory on a 
Lattice 1981 

[77] Peter Orland. Extrinsic curvature dependence of Nielsen-Olesen 
strings. Nucl. Phys. b\28, page 221, 1994. hep-th/9404140. 

[78] Jiri Hosek and Georges Ripka. Macroscopic chromodynamics of a con- 
fining medium. Zeit. Phys. A354, page 177, 1996. hep-ph/9412285. 

[79] Ulrich Ellwanger and Nicolas Wschebor. Solvable confining gauge the- 
ories at large. Phys. Lett. B517, page 462, 2001. hep-th/0107093. 

[80] Ulrich Ellwanger and Nicolas Wschebor. Confinement with Kalb- 
Ramon fields, hep-th/0107196, 2001. 

[81] Shinichi Deguchi and Yousuke Kokubo. Quantization of massive 
abelian antisymmetric tensor field and linear potential. Mod. Phys. Lett. 
A17, page 503, 2002. hep-th/0204087. 

[82] D.Antonov and E.Ebert. Confining properties of abelian (-projected) 
theories. Eur. Phys. J. C12, page 349, 2000. hep-th/9812112. 

[83] D. Antonov and D. Ebert. String representation of field correlators in 
the dual abelian Higgs model. Eur. Phys. J. C8, page 343, 1999. 

[84] Y.Nambu. Magnetic and electric confinement of quarks. Phys. Rep. 
23, page 237, 1976. 

[85] Filippo Neri Richard A. Brandt and Daniel Zwanziger. Lorentz invari- 
ance from classical particle paths in quantum field theory of electric 
and magnetic charge. Phys. Rev. D19, page 1153, 1979. 

[86] M.I. Polikarpov F.V. Gubarev and V.I. Zakharov. Monopole- 
antimonopole interaction in the abelian Higgs model. Phys. Lett. E>438, 
page 147, 1998. hep-th/9805175. 

[87] P. Senjanovic. Path integral quantization of field theories with second 
class constraints. Annals of Physics 100, page 227, 1976. 



141 



[88] M. Blagojevic and P. Senjanovic. One-potential formulation of the 
quantum field theory of magnetic poles. Nucl. Phys. B161, page 112, 
1979. 

[89] A.M. Polyakov. JETP. Letters 20, page 194, 1974. 
[90] G. T' Hooft. Nucl. Phys. B79, page 276, 1974. 

[91] G. 't Hooft. Topology of the gauge condition and new confinement 
phases in non-abelian gauge theories. Nucl. Phys. B190, page 455, 
1981. 

[92] M.Lavelle A.Kovner and D.McMullan. On the nonexistence of magnetic 
charge in pure yang-mills theories. 2002. hep-lat/0211005. 

[93] F.V. Gubarev and V.I. Zakharov. Interpreting the lattice monopoles 
in continuum terms. 2002. hep-lat/02 11033. 

[94] F.V.Gubarev and V.I. Zakharov. Gauge invariant monopoles in SU(2) 
dynamics. 2003. hep-lat/0204017. 

[95] F.V.Gubarev. Gauge invariant monopoles in lattice SU(2) dynamics. 
2002. hep-lat/0204018. 

[96] M.I. Polikarpov M.N. Chernodub, F.V. Gubarev and V.I. Zakharov. 
Monopoles and confining strings in QCD. hep-lat/0103033, 2001. 

[97] M.I. Polikarpov M.N. Chernodub, F.V. Gubarev and V.I. Zakharov. 
Towards an abelian-like formulation of dual gluodynamics. Nucl. Phys. 
B600, page 163, 2001. hep-th/0010265. 

[98] A.J. Van der Sijs. Abelian projection without ambiguities. 
Prog. Theor.Phys.Suppl. 131, page 149, 1998. hep-lat/9803001. 

[99] G.Shierholz A.S.Kronfeld and U.J.Wiese. Topology and dynamics of 
the confinement mechanism. Nucl. Phys. B293, page 461, 1987. 

[100] M;N; Chernodub and M.I. Polikarpov. Abelian projection and 
monopoles. In Cambridge 1997, Confinement, Duality and Non- 
Perturbative Aspects of QCD., page 387, 1997. hep-th/9710205. 



142 



[101] Gunnar S. Bali. QCD forces and heavy quark bound states. Phys. Rep. 
343, page 1, 2001. hep-ph/0001312. 

[102] Hiroko Ichie and Hideo Suganuma. Dual higgs theory for color confine- 
ment in quantum chromodynamics. 1999. hep-lat/9906005. Based on 
the doctoral thesis of H. Ichie. 

[103] Hiroko Ichie and Hideo Suganuma. Monopoles and gluon fields in QCD 
in the maximal abelian gauge. Nucl. Phys. B574, page 70, 2000. hep- 
lat/9808054. 

[104] William Tucker John D. Stack and Roy J. Wensley. The maximal 
abelian gauge, monopoles and vortices in SU(3) lattice gauge theory. 
Nucl. Phys. B639, page 203, 2002. hep-lat/0205006. 

[105] J. Greensite L. Del Debbio, M. Faber and S. Olejnik. Center dominance 
and z(2) vortices in SU(2) lattice gauge theory. Phys. Rev. D55, page 
2298, 1997. hep-lat/9610005. 

[106] Hiroshi Shiba Satoshi Kamisawa, Yoshimi Matsubara and Tsuneo 
Suzuki. A static baryon in a dual abelian effective theory of QCD. 
Nucl. Phys. B389, page 563, 1993. 

[107] M.N. Chernodub and D.A. Komarov. String representation of SU(3) 
gluodynamics in the abelian projection. JETP Lett. 68, page 117, 1998. 
hep-th/9809183. 

[108] M.N. Chernodub. Classical string solutions in effective infrared the- 
ory of SU(3) gluodynamics. Phys. Lett. B4?4, page 73, 2000. hep- 
ph/9910290. 

[109] K. Bardakci and S. Samuel. Local field theory for solitons. Phys. Rev. 
D18, page 2849, 1978. 

[110] P. de Forcrand C. Alexandrou and A. Tsapalis. The static three quark 
SU(3) and four quark SU(4) potentials. Phys. Rev. D65, page 054503, 
2002. hep-lat/0107006. 



143 



